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Resume 



Le cadre general dans lequel se place cette these est I'etude des proprietes de la matiere 
produite lors de collisions nuclcaires a haute cnergie, ou I'on s'attend a crcer artificiellement 
les conditions requises pour la formation d'un plasma de quarks et de gluons. En particulier, 
nous etudions certains aspects relies a revolution intrinsequement hors d'equilibre de ce 
type de systeme. 

La premiere partie est consacree a I'etude de la formation possible de condensats chiraux 
desorientes (DCC) lors du passage rapide de la transition de phase chirale. Nous calculous 
d'abord la toute premiere estimation de la probabilite de formation d'un champ de pion 
classique potentiellement observable, lors de I'expansion spherique rapide d'une "bulle" de 
matiere chirale chaude. Ce calcul necessite la construction d'une methode d'echantillonnage 
des conditions initiales pour le champ chiral. Ensuite, en pratiquant une analyse dctaillee 
de la structure d'isospin du champ classique, nous montrons que le modele le plus simple 
utilise jusqu'a present ne permet pas d'expliquer le caractere collectif de la configuration 
DCC, ce qui contredit une idee tres largement admise. 

Dans la seconde partie, nous etudions la thermalisation des gluons produits dans les tout 
premiers instants de la collision. Nous modelisons I'effet des collisions elastiques par une 
approximation de temps de relaxation auto-coherente, et comparons differentes conditions 
initiales proposees dans la litterature. Nous arguons que les criteres utilises dans les travaux 
antecedents pour caracteriser I'equilibration ne sont pas satisfaisant et proposons plutot 
de mesurer le degre d'anisotropie de differentes observables. Nos conclusions contredisent 
celles obtenues precedemment, nous montrons en particulier que les collisions elastiques 
sont insufHsantes pour thermaliser le systeme aux energies de RHIC. 



Abstract 



The general framework of the present thesis is the study of the properties of the matter 
produced in high energy heavy ion colhsions, where one expects to reach the conditions 
under which nuclear matter turns into a quark-gluon plasma. In particular, we study some 
aspects related to the intrinsic non-equilibrium evolution of such systems. 

The first part is dedicated to the study of disoriented chiral condensate (DCC) forma- 
tion during the out of equihbrium chiral phase transition. We compute the first estimation 
of the probability of forming a classical pion field after the rapid expansion of a spherical 
droplet of chiral matter. This requires the construction of a method for sampling the initial 
configurations of the chiral field. Then, by performing a detailed analysis of the isospin 
structure of the classical field, we show that the simplest model, used up to now, cannot 
explain the collective nature of the DCC. This contradicts a widely admitted idea. 

In the second part, we study the thermalization of initially produced gluons. We con- 
sider only elastic scatterings and use a self-consistent relaxation time approximation. We 
compare different initial conditions proposed in the literature. We argue that the criteria 
used in previous works to characterize equilibration are not reliable and propose instead 
to test the isotropy of various observables. Our conclusions are in contradiction with those 
previously obtained, we show in particular that at RHIC energies, elastic collisions are not 
effective enough for the system to thermalize. 



Sommaire 



Remerciements 
Resume 
Abstract 
Sommaire 

Introduction generale 

1 LES CONDENSATE CHIRAUX DJi!S0RiENTE5| 



1 DCC : introduction 

1.1 La symetrie chirale) 

1.2 Un LASER a pions 

1.3 Un scenario microscopique . . . 

1.4 Expansion et conditions initiales 

1.5 Effets quantiques 



2 La probabilite de formation 
d'un champ de pion classique 

2.1 Le formalism^ 



2.1.1 Approximation de champ moyen 





2.1.2 Expansion splierique 


2.1.3 Particules physiques, champ interpolant 




2.2 Echantillonnage des conditions initiales 




2.2.1 Les conditions initiales 



2.2.2 Equilibre thermique loca^ 



2.3 Resultats et discussion 

2.3.1 La strategic : recapitulatif 

2.3.2 Coherence physique du modele et choix des parametres 

2.3.3 RcsiiTtat^ 

2.4 Conclusion et commentaire 



3 Description microscopique de la formation d'un DCC 

3.1 Le formalism^ 

3.1.1 Le modele de Rajagopal et Wilczek 

3.1.2 Observables 



3.2 Les resultats 



viii 



3.3 Discussion et conclusion 58 

|3.3.1 Resume 58 

|3.3.2 Comparaison avec des travaux antecedents 59 

|3.3.3 Speculation^ 60 

13.3.4 Conclusion 61 



A Creation de particules dans une| 



geometrie en expansion 



B Equilibre thermique et equilibre thermique local 




B.l Equilibre thermodynamique 






B.2 Equilibre thermique local 





63 

69 

69 



C Calcul d'integrales 



D Scenario du trempage : Petat initial 



E Modele gaussien 



75 
77 
83 



III VERS LA FORMATION D'UN PLASMA DE MAtIERE DE 



CONFINEE 



89 



4 Equilibration thermique des gluon^ 



91 

4.1 Le modele 97 

99 



4.2 



4.1.1 Lois de conservation 



4.1.2 Expansion .vs. Collisions] 100 



L' approximation du temps de relaxation] 101 

4.2.1 La methode des moments 102 



4.2.2 Procedure de resolution 103 



4.2.3 Le choix du moment 103 



4.3 Conditions initiales et resultats 105 



4.3.1 "Temps de vie" 105 



4.3.2 Le scenario de saturation 106 



4.3.3 Le scenario des minijets Ill 

4.3.4 Robustesse des resultats 112 



4.3.5 Conclusions et perspectives 113 



F L'integrale de collision de Landau 



G Le calcul des moments 



H Calcul de Ojto) dans le scenario des minijets 



Conclusion 



117 
123 
127 
129 



Reference^ 



136 



Introduction generate 



La matiere nucleaire, telle que nous la connaissons a Techelle sub-nucleaire, est remar- 
quablement bien decrite par la chromodynamique quantique (QCD), la theorie moderne 
des interactions fortes. Le bien-fonde de la description du monde hadronique en termes de 
quarks et de gluons, les degres de liberte fondamentaux du lagrangien de QCD, s'appuie 
sur des bases experiment ales et des fondements tlieoriques solides. Une des proprietes les 
plus importantes des interactions fortes est la liberte asymptotique : I'intensite effective de 
I'interaction entre les quarks et les gluons decroit avec I'echelle de distance a laquelle on la 
mesure. Ceci a entre autres pour consequence le fait qu'a tres haute temperature et/ou tres 
haute densite, les quarks et les gluons forment un plasma de matiere deconfinee faiblement 
couplee. L'etude des proprietes de la matiere nucleaire dans ces conditions "extremes", est 
un sujet de recherche qui progresse tres rapidement. Le diagramme de phase des interac- 
tions fortes commence a etre bien connu theoriquement, en particulier dans la region de 
potentiel chimique baryonique nulle. L'existence de nouvelles phases de la matiere est une 
prediction tout a fait remarquable de QCD. 

II existe certains systemes physiques ou ces conditions extremes sont naturellement 
realisees, comme par exemple dans le coeur des etoiles a neutrons, ou encore durant les 
premiers instant de I'univers primordial. De plus, on s'attend a pouvoir creer momentane- 
ment ces conditions "en laboratoire" lors de collisions entre noyaux lourds a tres haute 
energie. De nombreuses experiences ont eu lieu a I'Alternating Gradient Synchroton (AGS, 
BNL) et au Super Proton Synchrotron (SPS, CERN) durant ces dix dernieres annees, 
d'autres se deroulent actuellement au Relativistic Heavy Ion Collider (RHIC, BNL), et un 
programme est prevu a partir de 2005 au Large Hadron Collider (LHC, CERN), avec des 
energies de plus en plus eleveesQ. De nombreuses signatures de la matiere deconfinee ont 
ete proposees et sont activement etudiees tant experiment alement que theoriquement (voir 
par exemple |[l|). La dynamique d'une collision d'ions lourds est cependant extremement 
compliquee et il est tres difficile d'interpreter les donnees de maniere non-ambigiie. En 
particulier les donnees obtenues au SPS, si elles fournissent de fortes indications en faveur 
de la formation de matiere deconfinee, ne permettent pas d'affirmer quoi que ce soit de 
definitif. 

C'est dans ce cadre general que se place le sujet de cette these. Nous etudions cer- 
tains aspects de la physique des collisions d'ions lourds, relies en particulier au fait que 
la matiere produite n'est a priori pas a I'equilibre thermodynamique. Notre interet pour 
cette problematique est double. D'une part, certains phenomenes qui n'ont lieu que dans 
des situations hors d'equilibre pourraient donner lieu a des signatures uniques, fournissant 
ainsi une source d'informations supplementaire (ou complementaire). C'est le cas de la for- 
mation de condensats chiraux desorientes (DCC) lors du passage rapide de la transition de 
phase chirale. La detection de ce phenomene hypothetique pourrait par exemple etre utile 



^Collisions Au-Au avec — 200 GeV par nucleon a RHIC et Pb-Pb avec = 5.5 TeV par nucleon a 
LHC!! 
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pour localiser cette transition dans le diagramme des phases. Dans la premiere partie de 
cette these, nous etudions la possibilite qu'un tel phenomene se produise. Un DCC est une 
region de I'espace oil le parametre d'ordre de la symetrie chirale oscille dans une direction 
differente de celle qu'il a dans le vide. C'est une configuration collective classique du champ 
de pion. Nous calculous la toute premiere estimation de la probabilite pour qu'un champ 
classique potentiellement observable soit cree lors de I'expansion spherique rapide du sys- 
teme. Pour ce faire, nous proposons une methode originale d'echantillonnage des conditions 
initiales, habituellement choisies de maniere arbitraire dans la litterature. Nous obtenons 
une limite superieure de cette probabilite relativement faible, typiquement de I'ordre de 
10-3 (Chap. I). 

Nous nous interessons ensuite a I'aspect collectif du DCC : tous les modes du champ 
(classique) oscillent dans la meme direction de I'espace d'isospin. Cette propriete est a 
I'origine de la signature la plus remarquable du phenomene, qui est aussi la plus utilisee 
dans les strategies experiment ales de detection : la distribution anormalement large de la 
fraction neutre du nombre total de pions emis. Nous montrons que le champ produit dans le 
modele microscopique le plus simple, utilise jusqu'a present, n'exhibe pas ce comportement 
collectif : les differents modes sont comme autant de DCC dont les orientations dans I'espace 
d'isospin sont independantes les unes des autres. Ce resultat contredit une idee largement 
admise et remet en question la possibilite de former un DCC dans une collision d'ions 
lourds (Chap. ^). 

La deuxieme raison pour laquelle il est important d'etudier les aspects hors d'equilibre 
dans ces collisions est le fait que les calculs actuels concernant les signatures experiment ales 
reposent sur I'hypothese selon laquelle le systeme est en equilibre thermique local. II est 
important de savoir si cette hypothese est justifiee dans une collision reelle. Si ce n'est pas 
le cas, il est alors necessaire d'avoir une idee de la fagon dont le systeme evolue pour pouvoir 
interpreter les donnees. Dans la deuxieme partie de la these, nous nous interessons a I'equili- 
bration thermique du systeme de gluons initialement produits (Chap. Nous decrivons 
la dynamique des gluons a I'aide d'une equation de Boltzmann que nous modelisons par 
une simple approximation de temps de relaxation. Nous comparons differents scenarios 
proposes dans la litterature pour decrire I'etat initial : les scenarios des minijets et de 
saturation. En calculant le temps de relaxation de fagon auto-coherente nous reproduisons 
de maniere semi-quantitative certains resultats exacts, recemment obtenus numeriquement 
dans le scenario de saturation. Pour caracteriser I'ecart a I'equilibre nous mesurons le degre 
d'anisotropie du systeme. Nos conclusions contredisent celles de travaux precedents ou des 
criteres differents sont utilises pour mesurer le degre d'equilibration. En particulier nous 
montrons que les collisions elastiques ne sont pas suffisantes pour thermaliser le systeme 
aux energies de RHIC. 

Au cours de notre presentation, nous essayons de mettre en avant les idees physiques 
sous-jacentes aux problemes etudies, releguant les details techniques ou les derivations non 
essentielles mais instructives dans les annexes. Les differents chapitres sont relativements 
independants. Pour des raisons pratiques et esthetiques, nous avons choisi des notations 
legerement differentes d'un chapitre a I'autre, notamment en ce qui concerne les indices 
chiraux. Nous travaillons partout avec le systeme d'unites h = c = ks = 1 {ou ks est la 
constante de Boltzmann). 



Premiere partie 



LES CONDENSATS CHIRAUX 
DESORIENTES 



Chapitre 1 

DCC : introduction 



Dans ce chapitre nous passons en revue certaines idees pertinentes a la physique des 
condensats chiraux desorientes. Apres un bref rappel historique des motivations initiales, 
nous presentons les principaux developpements qu'a connu le sujet depuis sa naissance, au 
debut des annees 1990. Nous presentons divers concepts et images physiques qui constituent 
la toile de fond des travaux presentes dans la premiere partie de cette these. 

1.1 La symetrie chirale 

La theorie des interactions fortes, ou chromodynamique quantique (QCD) est approx- 
imativement invariante sous les operations du groupe SUl{2) x SUr{2) qui consistent en 
des transformations unitaires independantes des composantes chirales droite et gauche 
du doublet de quarks legers g = {u,d), c'est la symetrie chirale. Elle est realisee de 
fagon spontanement brisee, phenomene decrit par la valeur non nulle du condensat de 
quarksQ {0\qLqji\0) et se manifestant dans le spectre d'excitations des interactions fortes 
par la faible masse des pions 140 MeV) comparee a I'echelle de masse typique des 
hadrons (~ 1 GeV). Les proprietes de transformation du produit bilineaire qlQr sont 
analogues a celles d'un vecteur cf) = (a, n) sous les rotations quadri-dimensionelles du 
groupe 0(4). Dans ce language, les composantes, dites d'isospin, 7^1,^2,1^3, representent 
les modes collectifs ou bosons de Goldstone accompagnant le phenomene de brisure spon- 
tanee : les excitations elementaires (quanta) du champ vr sont les pions. La valeur moyenne 
dans le vide du champ chiral (p pointe dans la direction it : on a (O|0|O) = (/7r,0), ou 
/tt = 92.5 MeV est la constante de desintegration du pion. La symetrie residuelle 0(3) 
des rotations autour de cet axe est la symetrie d'isospin des interactions fortes. De fagon 
generale, la valeur moyenne du champ 0, qui mesure I'intensite de la brisure spontanee, 
est appelee parametre d'ordre de la symetrie chirale. 

1.2 Un LASER a pions 

Un condensat chiral desoriente (DCC) est un etat dans lequel le parametre d'ordre a une 
orientation differente de celle qu'il a dans le vide. Plus precisement, c'est une configuration 
classique ou le champ de pion oscille de maniere coherente dans une direction donnee de 
I'espace d'isospin. Pour preciser la nature du phenomene et en comprendre I'interet, il est 
utile de faire une petite retrospective historique. 

'11 existe d'autres parametres d'ordre associes a la brisure spontanee de la symetrie chirale. Nous ne les 
considererons pas ici. 
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Chapitre 1. DCC : introduction 



Dans le contexte des collisions hadroniques ou nucleaires a tres haute energie, la pro- 
duction multiple de pions dont les energies transverses ne depassent pas quelques centaines 
de MeV est chose courante. Certains auteurs ont propose d'interpreter ce phenomene comme 
resultant de la desintegration d'un etat coherent ||2[ ||, ^, |5|, I'idee (deja presente dans 
d'anciens articles de W. Heisenberg |^) etant de voir le processus de production multiple 
comme le rayonnement d'un champ classique. On parle aussi d'etat multipions ou encore 
de LASER a pions. Cette idee est longtemps restee marginale jusqu'au debut des annees 
1990 ou, sur des bases theoriques plus solides, elle a ete largement etudiee et developpee. 

En effet, le developpement de theories effectives de basse energie dans le milieu des 
annees 80 a permis d'etudier la dynamique des excitations de grande longueur d'onde des 
interactions fortes dans le cadre de modeles bases sur des fondements theoriques solides^. 
L'exemple le plus simple est le modele a non-lineaire, qui decrit bien les proprietes des 
pions de tres basse energie. En termes du champ (f) introduit plus haut, Taction de ce 
modele s'ecrit 

5 = 1 y" d^xd^cf>-df'cf> (1.1) 

avec la contrainte 

<l^' = a' + n' = f^ (1.2) 

La pertinence de ces modeles pour le probleme de I'existence de champs de pions clas- 
siques est double : I'approximation classique est justifiee par le fait qu'on s'interesse a des 
excitations de grande longueur d'onde, et I'etude des solutions classiques permet de savoir 
quelles configurations sont autorisees par la dynamique sous-jacente. C'est essentiellement 
ce deuxieme point qui manquait aux etudes precedentes. C'est dans ce cadre conceptuel, 
que differents auteurs etudient les solutions classiques du modele a non-lineaire au debut 
des annees 1990 |@, |9|. Ces etudes montrent en particulier qu'il est naturel de consid- 
erer des configurations classiques du champ de pion de grande longueur d'onde ayant des 
geometries non-triviales dans I'espace d'isospin. Illustrons ce point par un exemple. En 
utilisant des conditions aux limites idealisees modelisant la region centrale d'une collision 
tres energetique J. P. Blaizot et A. Krzywicki se ramenent a un probleme a une di- 

mension : le systeme est invariant sous les boost longitudinaux et le champ ne depend que 
du temps propre r = Vf^ — z"^, ou z est I'axe de la collision. Les equations de conservation 
des courants vectoriel = tt x d^n et axial = Trd^a — ad^n 

^f,y^' = o ■ a^A^ = o 

sont facilement integrees : 

Vq = TT X TT = — ; Aq = TTO" — (TTT = — 

T r 
ou les vecteurs a et b (a-b = 0) sont des constantes d'integration. Le champ de pion decrit 

^Les proprietes de symetrie chirale etant a I'origine de I'existence des excitations de grande longueur 
d'onde (les modes de Goldstone), il est clair qu'elles jouent un role important quant a leur dynamique. 
En exploitant les proprietes de symetrie chirale du lagrangien fondamental de QCD, on peut ecrire un 
lagrangien efFectif de basse energie qui se presente sous la forme d'un developpement en gradients (qui 
correspond a un developpement en puissances de Techelle d'energie typique a laquelle on s'interesse) et 
qu'il est necessaire de tronquer en pratique. Le modele a non-lineaire, dont il est question dans la suite, 
est obtenu en ne retenant que le premier terme, qui contient deux derivees. 



1.2. Un LASER a pions 
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done une trajectoire elliptique dans I'espace d'isospin (c = a x b, = 



TTa = 













K In 




TTft = 


— sin 1 








To 
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Kin 
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VTc = 


— cos 






K 




TO 



oH To est une constante delimitant I'hypersurfaee sur laquelle sont loealisees les sources du 
champ classique de pion ||9|. Cette solution appartient a la classe plus generale des solutions 
de type onde planes proposee par A. A. Anselm (voir aussi pl ). 

Parallelement, J.D. Bjorken suggere la possibilite de former une configuration partic- 



uliere du champ de pions qu'il denomme condensat chiral desoriente [T^, |T^]. II donne une 
image intuitive de la collision et de la formation du DCC qu'il est bon d'avoir en tete car 
elle est tres parlante, c'est le scenario "Baked- Alaska" : les produits primaires de la colli- 
sion s'eloignent de la zone d'impact avec une vitesse proche de celle de la lumiere, formant 
une "boule de feu" en expansion rapide et dont I'interieur est isole du vide environnant. 
Si la densite d'energie y est sufRsament faible, I'interieur ressemble de tres pres au vide. 
Cependant, ce vide n'etant pas en contact avec le vide physique, I'orientation du parametre 
d'ordre n'a pas de raisons d'y etre la meme, en particulier les composantes d'isospin de ce 
dernier peuvent etre non nulles ; c'est un vide desoriente. Apres un certain temps (de I'ordre 
de quelques fermi), la surface de la boule de feu se desagrege (hadronisation), I'interieur et 
I'exterieur ne sont plus separes et le vide desoriente relaxe vers le vrai vide en rayonnant 
ses modes collectifs, les pions. Ce scenario phenomenologique correspond a I'idealisation 
adoptee par Blaizot et Krzywicki dans le cas d'une expansion longitudinale. Un cas parti- 
culier de leur solution correspond a un champ oscillant dans une direction quelconque de 
I'espace d'isospin (c'est le cas oii I'intensite initiale du courant axial est grande comparee 
a celle du courant vectoriel : 6 » a. Le champ de pions oscille alors dans la direction b.). 
C'est le DCC. La prediction la plus frappante de ce phenomene hypothetique est la distri- 
bution evenement par evenement de la proportion de pions neutres emis par rayonnement. 
En effet, le nombre total de pions avec une composante d'isospin donnee est proportionnel 
(dans le cas du DCC uniquement, voir le Chap. ^) a I'integrale du carre du champ sur tout 
le volume de la bulle 

Jv 

Les interactions fortes etant invariantes sous les rotations d'isospin, toutes les directions 
sont equiprobables et la proportion de pions neutres 



est alors distribuee selon la loi^ 

dPjf) _ 1 
df 2^7' 



(1.3) 



^Cette distribution est donnee explicitement dans |^ ou les auteurs n'ont pourtant pas isole la config- 
uration DCC, ce qui n'est pas correct. Ceci est corrige dans |^. Notons que ce resultat est deja present 
dans 1^] ou I'auteur considere un etat coherent d'isospin total nul. Le DCC est un tel etat, en efi^et le champ 
oscille dans une direction donnee et est, par consequent parallele a sa vitesse, ce qui signifie que I'isospin 
total de cette configuration est nul : V = ^ x = 0, et I = /y V. 
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Chapitre 1. DCC : introduction 



a comparer avec la distribution binomiale tres piquee autour de la valeur moyenne / = 1/3 
prevue dans le cas de la production incoherente de pions. Ce phenomene pourrait expliquer 
les evenements Centauro et anti-Centauro enregistres dans le domaine des rayons cosmiques 
de tres haute energie, dans lesquels des amas de pions tous charges (Centauro) ou tons 
neutres (anti-Centauro) ont ete observes |15|. 

Si on comprend mieux les configurations classiques possibles, compatibles avec la dy- 
namique des interactions fortes, on ne salt cependant pas comment un tel champ pent etre 
cree lors de la collision. En effet, dans les considerations precedentes, on a suppose I'exis- 
tence d'un champ coherent, classique cree par des sources localisees sur le cone de lumiere 
(la surface de la boule de feu). La question de la formation du champ classique est cachee 
dans I'hypothese de I'existence de ces sources. 



1.3 Un scenario microscopique 

L'idee du DCC regoit une impulsion considerable quand, en 1993, K. Rajagopal et 
F. Wilczek proposent un scenario prevoyant I'emergence d'un champ fort de pions apres 



une collision de noyaux lourds ultra-relativistes |16|. Dans ce type de collisions, la densite 
d'energie par unite de volume deposee dans la region centrale pent atteindre des valeurs 
tres elevees (plusieurs GeV/fm^) et on s'attend a ce que la matiere soit dans un etat ou la 
symetrie chirale est restauree^ du fait des fortes fluctuations du champ chiral. L'expansion 
rapide du systeme engendre une chute brutale de la densite d'energie et par consequent une 
suppression soudaine de ces fluctuations. Le parametre d'ordre "se fige" dans une direction 
aleatoire de I'espace d'isospin. Ce phenomene est analogue a la formation de domaines 
d'aimantation lors du trempage d'un materiau feromagnetique. C'est l'idee sous-jacente au 
scenario de Rajagopal et Wilczek, aussi appelle scenario du trempage, et que nous allons 
maintenant decrire un peu plus en detail car il est a la base des travaux decrits dans 
les chapitres suivant. En fait, ce modele sera etudie dans le chapitre § ou nous aurons 
I'occasion d'en donner tous les details. Pour le moment, il est suffisant d'en presenter les 
grandes lignes et les resultats importants. 

Considerons la region centrale d'une collision d'ions lourds ultrarelativistes, c'est a dire 
la region autour de z = ou z est I'axe du faisceau dans le referentiel du centre de masse 
de la collison. A tres haute energie, les nucleons des noyaux incidents ne sont pratiquement 
pas ralentis et ont completement evacue la region centrale apres un temps tres bref. Cette 
derniere est alors formee de matiere non-baryonique ||l^ (la densite de baryons est egale a 
celle des anti-baryons). Nous nous interessons a la dynamique des modes de grande longueur 
d'onde et nous considererons le champ de pions dans cette region. Nous negligeons en 
particulier les effets possibles dus a la presence d'autres mesons ou baryons ainsi qu'a celle 
de matiere deconfinee. Nous supposerons de plus que le systeme a atteint un etat d'equilibre 
thermique local apres un certain temps. Si la densite d'energie dans la region centrale est 
suffisante, la symetrie chirale est restauree. Pour decrire revolution ulterieure du champ 
de pion, nous devons utiliser un modele capable de prendre en compte la possibilite d'un 
etat symetrique 0(4). II est clair, du fait de la contrainte (|1.2|) que ce n'est pas le cas du 
modele a non-lineaire. L'extension la plus simple consiste a relacher la contrainte : c'est le 



■*Supposons pour un instant que le systeme est decrit par un gaz de pions de masse nuUe a I'equilibre 
thermodynamique. La temperature de ce gaz de pions (relativistes) est reliee a la densite d'energie par la 
loi : e = 3(71^/30) T**, on g — 3 est la degenerescence du triplet de pions. A une densite d'energie e « 1 
GeV/fm"^, correspond une temperature T ~ 400 MeV. La temperature critique de la transition de phase 
chirale est Tc ~ 150 MeV. 
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modele cr-lineaire^ 

5 = Jd'^x jl^^^.a'^^-^ {((> ■ (!> - v^f + Ha^ . (1.4) 

Le terme Ha brise explicitement la symetrie et rend compte de la masse non-nulle des 
pions. Au niveau classique, on a 

H = Uml 
ml = \{fl-v^) 
ml = \{2,fl-v^) 

On retrouve le modele non-lineaire en prenant la limite chirale [m-^ — > 0) et en integrant 
le degre de liberte lourd a. 

Rajagopal et Wilczek ont etudie numeriquement la dynamique du parametre d'ordre 
dans le cadre du modele a lineaire classique en supposant un etat initial symetrique : les 
composantes du champ sont des variables aleatoire gaussiennes de meme variance et 
distribuees independemment sur les differents noeuds d'un reseau cubique. Bien qu'ayant 
en tete un systeme en expansion, Rajagopal et Wilczek travaillent dans une configura- 
tion statique. Le trempage des fluctuations initiales est modelise par la faible valeur de 
la variance dans I'etat initial. Durant revolution ulterieure du champ de pion, les modes 
de Fourier de grande longueur d'onde sont fortement amplifies en comparaison des modes 
de plus petite longueur d'onde, et oscillent dans le temps de maniere coherente avec une 
frequence 27r/ sj m^ + k"^ 2T:/m.,^. L'energie est ensuite repartie de maniere equivalente 
entre les modes (equipartition) par la dynamique non lineaire. Cependant dans une situa- 
tion en expansion, la dilution du systeme entraine le decouplage des modes (freeze-out). S'il 
a lieu suffisament tot, celui-ci pent empecher le systeme d'atteindre I'equipartition. Dans 
ce cas, la configuration finale du champ consiste en une superposition de modes de grande 
longueur d'onde et de grande amplitude. C'est le champ fort dont nous avions besoin. Le 



resultat de la Ref. [|l^] est reproduit sur la Fig. |Ll| ou Ton voit revolution temporelle du 
module au carre des composantes de Fourier de (j)'^, moyenne sur des bins de largeur 6k 
centres en \ \k\ \ = k pour diffferentes valeurs de k. 

Le phenomene decrit plus haut pent etre compris qualitativement a I'aide d'une simple 



approximation |16, 19, Les equations du mouvement s'ecrivent 

{d^ + X{ct>^ - v^)) (f} = Hn^ , (1.5) 

ou I'on a utilise la notation (jP' = (f) ■ (f) , et ou est un vecteur unitaire dans la direction 
a de I'espace chiral. En remplagant le terme non lineaire cfP' par sa valeur moyenne sur 
le volume (dans ce qui suit, nous denotons la valeur moyenne spatiale par des crochets : 
{O) = (1/1^) Jy d^xO{x)) et en prenant la transformee de Fourier de I'equation obtenue, 
on obtientQ 

-^ + k'' + mljf{t)W{k,t) = Q, (1.6) 



^Le modele a-lineaire a 6t6 propose par M. Gell-Mann et M. Levy en 1960 ||18| (en fait avant, par 
Schwinger, voir dans pour decrire les interactions fortes entre nucleons. Ces derniers etaient representes 
par le doublet d'isospin (p, n) et I'interaction pion-nucleon par un couplage de Yukawa. 

^Cette approximation est equivalente a la limite A'^ oo ou A'^ est le nombre de composante du champ 
(A'^ = 4), c'est une approximation de type champ moyen. Nous aurons I'occasion d'en parler d'avantage au 
chapitre B. 
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temps 

Fig. 1.1: Le module au carre moyen des composantes de Fourier de 03 en fonction du temps (en 
unites du pas du reseau a) dans le scenario du trempage de la Ref. ||l6|. On a calcule la moyenne 
du module au carre de tous les modes tels que k = \ \k\\ est dans le bin de largeur Sk — 0.057a~^ 
centre autour de (de haut en bas sur la figure) ka — 0.20, 0, 26, 0, 31, 0, 37, 0, 48, 0, 60, 0, 71, 0, 94, 
1,16 et 1,39. 



ou 



(1.7) 



est la masse effective (au carre) des excitations du champ a I'instant t. File represente 
la courbure instantanee du potentiel effectif vu par le champ. 



L'evolution temporelle de mgjj(t) dans le scenario de la Ref. |16| est represente sur la 



Fig. |L2[ Quand {4>^){t) < cette masse effective est imaginaire pure, ce qui se traduit 



par une croissance exponentielle des amplitudes des modes de Fourier tels que k < —m^j^ 
(cf.Fq. (p^). Ce mecanisme est bien connu dans divers domaines de la physique, comme 
la physique des solides ou encore la fragmentation nucleaire, sous le nom d'instabilite spin- 
odale. Mentionnons le fait que le mecanisme de I'instabilite spinodale n'explique pas toute 
I'amplification observee en resolvant numeriquement les equations (|1.5|). Fn effet, on pent 



voir sur la Fig. |L2| que le carre de la masse effective ( |1.7D n'est negatif que pour des temps 
^ ^ 10 (sn unite de pas du reseau) alors que I'amplification des modes de basse frequence 
se poursuit jusqu'a t ~ 50. On pent toutefois comprendre qualitativement I'amplification 
ulterieure a I'aide du mecanisme de resonance parametrique du aux oscillations regulieres 
de la valeur moyenne (cr) autour de sa valeur asymptotique |21, 22, 23 1. 



1.4 Expansion et conditions initiates 

Le scenario du trempage fournit done un scenario microscopique simple pour la forma- 
tion possible d'un champ fort de pion lors d'une collision entre noyaux tres energetiques. 
De ce fait, il permet d'obtenir des informations au moins qualitatives sur les proprietes du 
champ classique ainsi produit. Le resultat de Rajagopal et Wilczek a par la suite ete con- 
firme 25 1 et le modele du trempage a regu un interet croissant. Diverses ameliorations 
y ont ete apportees, les plus importantes etant I'inclusion de certains effets quantiques 
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Fig. 1.2: La masse effective au carre ^ tfti f) en fonction du temps (en unite du pas du reseau a) 
dans le scenario du trempage de la Ref. ||l6|. 



et I'abandon de I'hypothese concernant les conditions initiales (trempage "a la main" des 
fluctuations initiales) en faveur d'un scenario prenant explicitement en compte I'expansion. 

La modelisation de I'expansion consiste simplement a se placer dans le referentiel en 
co-mouvement, qui se deplace avec un element de volume du systeme. Par exemple dans 
le cas d'une expansion longitudinale a la vitesse de la lumiere, le temps mesure dans le 
repere en co-mouvement est le temps propre r = \/{^ — : au temps t, le systeme a 
une extension spatiale —t<z<t. Autrement dit, on specific les conditions aux bords, 
necessaires a la resolution des equations du mouvement, non pas sur un hyperplan t =cte, 
mais sur une hypersurface r =cte (voir le modele de Blaizot et Krzywicki decrit plus 
haut). En definissant la variable de rapidite r\ = (1/2) ln(t — z)/{t + z), on pent reecrire 
dans (|L5|) 9f — Of = 9^ + (1/r) dr — (1/t^) d"^. En plus du terme habituel d'acceleration 
of, on voit apparaitre un terme de friction oc qui traduit le fait que la densite d'energie 
dans un co- volume decroit du fait de I'expansion. Differents auteurs [|^, 27, 24 1 ont etudie 
numeriquement les solutions de ( |1.5| ) dans une geometric en expansion longitudinale, avec 
un etat initial stable (ou la masse effective (|l.7D est initialement positive). II ressort de ces 
etudes que, pour des conditions initiales appropriees, I'expansion entraine le systeme dans 
la region d'instabilite et le champ fort est cree. Cependant, le phenomene depend fortement 
des conditions initiales. 

Dans la Ref. |28|, J. Randrup propose une methode systematique d'echantillonnage des 
configurations initiales du champ chiral 0(x) a partir d'un ensemble thermique. La fonc- 
tion de partition du systeme en equilibre thermique a la temperature T dans le volume V 
est calculee en traitant les fluctuations thermiques du champ autour de sa valeur moyenne 
spatiale (</>) dans une approximation de champ moyen (I'approximation de Hartree). On 
pent alors calculer la distribution de probabilite des valeurs possibles de (</>) ainsi que celle 
des fluctuations 5cf){x) = cf){x) — (0). Dans un article ulterieur |^, Randrup etudie les tra- 
jectoires classiques generees par les equations (|l.5|) avec des conditions initiales aleatoires, 
echantillonnees selon cette methode. La masse effective (Eq. ( |1.7|) ) s'ecrit ((50) = 0) 



X[{ct>f + {5cl?)-v'') , 



ou {(p)'^ = {4>) • {4>). La region d'instabilite fn^jj < est representee sur la Fig. L3 dans 
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le plan ((0), (5^^)^/^). Randrup represente ses trajectoires dans ce meme plan et etudie 
I'incursion du systeme dans la region d'instabilite en fonction de la temperature initiale 
d'une part et du taux d'expansion d'autre part. En effet, supposons une expansion D- 
dimensionnelle, le temps mesure dans le repere en co-mouvement est alors : r = ^/t^ — 
ou i = 1...D, et le terme de friction dans les equations du mouvement devient {D/T)dr : 
I'expansion est d'autant plus efficace que D est grand {D < 3). 



<0(^ > 




m 3,<0 



Fig. 1.3: Representation shematique de deux trajectoires typiques dans le plan {{(j)), \J {5(p))- 
L'evenement (1) traverse la region spinodale {^'iff < 0) et subit une forte amplification. Dans 
la region d'instabilite, apres le refroidissement brutal des fiuctuations (trempage), la valeur des 
fiuctuations du champ est ~ t;/2, comme dans le modele ad hoc de 

Randrup observe que lorsque le systeme est prepare a une temperature T = 400 MeV, il 
n'entre dans la region d'instabilite que pour D > 1. Pour D = 3, I'instabilite apparait pour 
des temperatures allant de 200 MeV a 500 MeV (au moins). Au debut de revolution, I'ex- 
pansion rapide supprime fortement les fluctuations thermiques tandis que {(/)), initialement 
petit pour des temperatures superieures a la temperature critique Tc ~ 150 MeV, ne change 
que tres peu. Le systeme entre dans la region d'instabilite. Les modes de grande longueur 
d'onde, et en particulier ((/>), le mode zero, sont alors fortement amplifies et le systeme est 
"pousse" hors de la region d'instabilite. Si la temperature initiale est trop importante, les 
fluctuations n'ont pas le temps de decroitre suffisamment avant que {4>) n'augmente (vers 
sa valeur dans la phase brisee) et le systeme n'est jamais instable. A I'inverse, si la tem- 
perature est trop faible, on est trop proche de la temperature critique et {(p) est si grand 
que Ton passe a cote de la region interessante. II est interessant de remarquer que, dans le 
cas ou le systeme entre dans la region d'instabilite, la valeur typique des fluctuations apres 
la breve periode de refroidissement (trempage) est {dcjP'Y/'^ ~ t;/2, ce qui correspond a la 
configuration initiale ad hoc de Rajagopal et Wilczek 



29|. 



Nous comprenons done comment, a partir d'une configuration chaotique du champ, 
oil les differents modes sont initialement independants, la dynamique hors d'equilibre du 



1.5. EfTets quantiques 
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systeme genere moment anement une configuration formee par la superposition de modes 
de grande longueur d'onde et de grande amplitude. 



1.5 EfFets quantiques 

Le plienomene d'amplification des modes de grande longueur d'onde dans notre systeme 
en expansion est analogue^ a celui de la creation massive de particules durant la phase d 'in- 
flation de I'univers primordial. Cette analogic a conduit differents auteurs a appliquer au 
premier probleme des methodes developpees pour le second, notamment en ce qui concerne 



la prise en compte de certaines corrections quantiques a la dynamique. Dans la Ref. |31| 



les auteurs travaillent avec un systeme statique en modelisant le trempage "a la main" de 



la meme fagon que dans |16|. Les fluctuations thermiques et quantiques sont decrites par 



une matrice densite supposee gaussienne, ce qui conduit a des equations du mouvement 



dont la structure est similaire a celle des Eqs. (|1.6D et (|1.8D , oil (...) doit etre maintenant 
compris comme la moyenne sur ces fluctuations. Les auteurs des Refs. |32, 33, 34 1 pren- 
nent en compte les effets quantiques dans I'approximation de champ moyen (equivalente a 
I'ansatz gaussien de ||3l[]) dans un modele avec expansion. Les resultats de ces etudes sont 
similiaires a ceux obtenus dans le cas classique^ : on observe une augmentation importante 
du nombre de quanta de basse energie durant la periode d'instabilite. Le champ est alors 
essentiellement une superposition de modes classiques de grande longueur d'onde. De meme 
que dans le cas classique, I'ingredient crucial est la dynamique fortement hors d'equilibre 
due a I'expansion et le phenomene d'amplification depend fortement des valeurs initiales 
des moyennes du champ et de sa derivee temporelle. Ce dernier point fera I'objet de I'etude 
presentee au Chap. ^, on nous construirons une methode originale d'echantillonnage des 
conditions initiales, que nous utiliserons pour calculer la probabilite de former un champ 
classique potentiellement observable. 

Dans le language de la mecanique quantique, le champ classique ainsi forme n'est autre 
que la valeur moyenne du champ quantique dans I'etat final du systeme, au moment du 
decouplage. Autrement dit, I'etat final est un etat quasi-classique, c'est a dire un etat 
coherent De fagon generale, independante de toute dynamique sous-jacente, on pent 
construire differents types d'etats coherents multi-pions prenant en compte les contraintes 
imposees par les lois de conservation, par exemple de la charge electrique ou encore de 
I'isospin m, ^, |5|. La structure de I'etat forme dans un scenario realiste depend de la dy- 
namique microscopique sous-jacente. Par exemple, il est facile de voir que I'approximation 
de champ moyen correspond a une description du systeme en termes des etats dits com- 
presses0 ("squeezed state") qui sont une generalisation des etats coherents |^]. Le DCC 
est un etat particulier ou tons les modes du champ sont orientes dans la meme direction 
de I'espace d'isospin (voir par exemple ||3^) : c'est un etat d'isospin total nul, ce qui est a 
I'origine de la loi ( |1.3[ ) ||5|. Pour decrire la formation d'un tel etat collectif dans un scenario 
microscopique donne, il faut prendre en compte les non-linearites de la dynamique. Au 
niveau quantique, cela recquiert de considerer des corrections au dela de I'approximation 
de champ moyen. En fait, nous verrons au Chap. ^ que le probleme se pose deja au niveau 
classique. En effet, nous montrerons que dans le modele de Rajagopal et Wilczek, ou I'etat 
initial est suppose incoherent, I'etat final n'est pas un etat collectif, les modes amplifies sont 



'^Dans le cas du DCC il s'agit cependant d'un probleme en couplage fort 

*Dans le modele avec expansioiijl'inclusion des effets quantiques donne des resultats quantitativement 
similaires a ceux du cas classique Pq ] 

^Ceci est directement relie aux transformations de Bogoliubov connectant les operateurs de creation et 
d'anihilation de quasi-particule a differents instants (voir Chap. 
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independants les uns des autres. Contrairement a une idee largement admise, le caractere 
collectif suppose du DCC n'a pas encore trouve d'explication microscopique. 

Jusqu'ici nous avons considere le modele le plus simple possible en faisant des approxi- 
mations parfois tres grossiere, ce qui nous a permis de degager I'essentiel de la physique du 
phenomene. Nombre de corrections ont ete considerees dans la litterature (tres abondante) 
sur le sujet[^. Mentionnons par exemple les effets dus aux interactions du champ de pion 



avec son environement dans une collision d'ions lourds (3^, la modelisation d'effets 

quantiques au dela du champ moyen par une source stochastique dans les equations du 
mouvement du parametre d'ordre [EH], ou encore I'etude de modeles alternatifs au modele 



cT-lineaire, comme les modeles de Nambu-Jona-Lasinio |41| ou de Gross-Neveu |42|. Ces 
etudes n'alterent pas I'image de la formation d'un champ classique que nous avons decrite 
ici, qui s'avere done etre relativement generique. 

La detection eventuelle d'un DCC dans des experiences de collisions hadroniques ou 
nucleaires a haute energie pourrait permettre d'obtenir des informations precieuses sur 
la structure du vide des interactions fortes, ainsi que sur les proprietes de la transition de 
phase chirale. Jusqu'a maintenant deux groupes ont tente de detecter ce phenomene : T864 
(minimax), une experience de collisions proton-antiproton a Fermilab |43| et WA98 au SPS 



du CERN (collisions d'ions lourds) |44|, toutes deux reportent des resultats negatifs. En 
fait a ce jour les seuls "signaux positifs" sont ceux reportes dans le domaine des rayons 
cosmiques au debut des annees 1980. Une partie du programme experimental du detecteur 
STAR|^ a RHIC (BNL) est dediee a la recherche de grandes fluctuations dans la distribution 
evenement-par-evenement de la fraction de pions neutres. Diverses autres signatures ont 
ete proposees dans la litterature, mentionnons par exemple I'augmentation du nombre de 
photons de basse energie ||4^, ^ par desintegration des pions neutres provenant du DCC, 
ou encore, plus recemment, I'anomalie dans les abondances de baryons multi-etrange Q et 

n m. 



Pour conclure ce chapitre d'introduction, reprenons les premieres questions-reponses de 
la 'DCC trouble list" de J. D. Bjorken, qui resument bien I'idee sous-jacente aux etudes 
presentees dans les deux chapitres suivant. 

Existence of DCC : 

Must it exist ? NO 

Should it exist ? MAYBE 

Might it exist ? YES 

Does it exist ? IT'S WORTH HAVE A LOOK 



"Voir la 'DCC home page' http ://wa98.web.cern.ch/WA98/DCC 
"Voir la 'STAR home page' http ://www.star.bnl.gov/STAR 



Chapitre 2 



La probabilite de formation 
d'un champ de pion classique 



Un scenario plausible pour la formation d'un champ de pion classique lors d'une collision 
nucleaire a haute energie a ete identifie |jl6| : I'expansion rapide du systeme entraine la 
suppression brutale des fluctuations thermiques (trempage), generant une amplification 
importante des modes de grande longueur d'onde. Ce mecanisme microscopique permet 
de faire des predictions qualitatives utiles pour I'elaboration de strategies experiment ales, 
comme par exemple le fait que le signal soit a rechercher dans les pions de basse energie. 
Un parametre crucial pour la phenomenologie est la probabilite pour qu'un tel champ fort 
soit produit. En effet on choisira une methode experimentale plutot qu'une autre selon que 
le nombre typique de collisions donnant lieu au phenomene cherche est de une sur dix ou 
de une sur un million. 

Dans ce chapitre, nous nous proposons d'estimer I'ordre de grandeur de ce nombre. 
II est clair que les modeles actuels ne permettent pas de dire quoi que ce soit de fiable 
au niveau quantitatif, neanmoins il est legitime de chercher a avoir une estimation, meme 
grossiere, de cette probabilite. En fait, nous en obtiendrons une limite superieure. Le travail 
expose dans ce chapitre a ete realise en collaboration avec Andre Krzywicki et a fait I'objet 



d'une publication dans la revue Physics Letters B |48| 



Cadre theorique et image physique 

Le scenario du trempage est relativement robuste en ce qui concerne les details de 
la dynamique. Ici, nous utiliserons le modele cr-lineaire en tant qu'approximation de la 
dynamique des excitations de basse energie de QCD que sont les pions. Le point essentiel 
est la forte dependance du phenomene d'amplification vis-a-vis des conditions initiales. 
Notre but etant d'estimer I'ordre de grandeur de la probabilite d'avoir une amplification 
donnee, il nous faut donner un poids statistique a chacune des configurations initiales 
possibles. 

Pour ce faire, nous aurons recours a des approximations, parfois tres simplificatrices, 
mais qui "capturent" I'essentiel de la physique du phenomene. Ainsi nous traiterons les ef- 
fets quantiques a I'ordre dominant dans un developpement en ou N est le nombre de 
composantes du champ chiral (f) {N = 4). C'est une approximation de type champ moyen 
qui, nous I'avons vu dans le chapitre precedent (Eqs. (|1.6|) et (|1.7D ), decrit la physique de 
I'instabilite spinodale, responsable du phenomene d'amplification. De plus, nous consid- 
ererons le cas d'une expansion spherique, la plus efficace pour le trempage des fluctuations 
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initiales. De maniere generale, nous nous placerons toujours dans le cas le plus favorable 
pour la formation d'un champ fort de pion, de fagon a avoir une limite superieure a la prob- 
abilite cherchee. En effet on s'attend d'ores et deja a ce que celle-ci soit faible et, notre 
objectif etant d'avoir une estimation qualitative, une limite superieure est un parametre 
pertinent pour la phenomenologie. 

Nous focalisons notre attention sur une petite bulle de matiere chirale chaude, for- 
mee lors de la collision, qui subit une expansion tres rapide et qui, selon son etat initial, 
subit eventuellement une periode d'instabilite durant laquelle les modes de grande longueur 
d'onde sont fortement amplifies. Nous supposons que les fluctuations statistiques du champ 
de pion a I'interieur de la bulle sont decrites par un ensemble thermique local. Exploitant 
cette hypothese, nous proposons une methode originale d'echantillonnage des configura- 
tions initiales du champ, habituellement choisies de fagon arbitraire dans la litterature. 



En combinant cette methode avec le formalisme de la Ref. |Q (voir aussi (3J]), nous cal- 
culous la toute premiere estimation de la probabilite pour qu'un champ fort (classique), 
potentiellement observable, soit forme dans une collision d'ions lourds. 

En toute rigueur, la probabilite obtenue doit etre interpretee comme une probabilite 
conditionnelle, a multiplier par la probabilite pour que notre bulle initiale soit formee dans 
une collision reelle. Le calcul de cette derniere quantite recquiert cependant un modele pour 
I'ensemble de la collision, ce qui depasse de loin le cadre de I'etude presentee ici. Neanmoins, 
I'hypothese d'equilibre thermique local pent etre vue, non pas comme une approximation de 
I'etat reel du systeme, mais comme une parametrisation de la distribution des configurations 
initiales possibles. Dans ce sens, le modele decrit ci-dessus est suffisamment generique pour 
pouvoir, malgre sa simplicite, donner une estimation raisonnable de la probabilite cherchee. 
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Le chapitre est organise comme suit : tout d'abord nous ecrivons les equations dy- 
namiques a I'ordre dominant en 1/N dans une geometrie en expansion (nous reprenons 



pour I'essentiel le formalisme des Refs. |33, |3^). Nous introduisons ensuite la notion de 
champ interpolant qui nous permettra de relier les nombres de particules dans les etats 
initial et final. La methode d'echantillonage des conditions initiales est exposee en details. 
Enfin nous presentons les resultats de notre calcul numerique, suivis d'une discussion et de 
la conclusion. 



2.1 Le formalisme 

La densite de Lagrangien du modele cj-lineaire dans une geometrie quelconque, carac- 
terisee par sa metrique g^u, s'ecrit (cf. (|1.4| )) 

ou g est le determinant de la metrique. L 'action S = f d'^xC est, comme il se doit, un 
scalaire sous les transformations du systeme de coordonees. 

Nous cherchons a etudier revolution hors d'equilibre d'un systeme quantique avec un 
nombre infini de degre de liberte. Le formalisme adapte est la theorie quantique des champs 
hors d'equilibre. La theorie quantique des champs a ete initialement developpee pour decrire 
des problemes de diffusion : on prepare le systeme a t = — cx) et on fait une mesure 
a t = +00, ou ±oo signifient "tres longtemps avant/apres I'interaction". II s'agit done de 
calculer des amplitudes de transition entre etats asymptotiques. La problematique qui nous 
occupe ici est tres differente : ayant prepare le systeme dans un certain etat initial, decrit 
par la matrice densite p{to), on cherche a suivre son evolution au cours du temps. On veut 
calculer des objets du type {0){t) = Tr{p{t)0), revolution temporelle etant donnee par 
les equations du mouvement p = [H, p], ou H est le hamiltonien du systeme. II existe un 
formalisme tres elegant qui permet d'utiliser les methodes developpees pour les problemes 
de diffusion dans le cas des problemes hors d'equilibre, c'est le formalisme dit de "chemin 



ferme" (Closed Time Path) |^, 50, 51 1. En particulier, cette formulation assure la causalite 
des equations obtenues. Les particularites propres a ce formalisme, comme par exemple la 
forme matricielle des propagateurs, ne se manifestant pas dans I'approximation de champ 
moyen, nous n'en parlerons plus dans la suite. 



2.1.1 Approximation de champ moyen 

L' approximation dite "de champ moyen" consiste a remplacer les interactions entre 
quanta (les modes de Fourier du champ, en nombre infini) par une interaction moyenne 
effective. Les modes sont ejfectivement decouples les uns des autres en ce sens que chacun 
d'eux ne "voit" les autres que par I'intermediaire du champ moyen qu'ils creent : chaque 
mode evolue dans un potentiel effectif cree par tons les autres modes. Nous avons deja vu 
un exemple de ce genre d'approximation dans le chapitre precedent (Eqs. (|l.6|) et (|T?^)). 
Rappelons-le, il nous sera utile dans la suite. Les equations du mouvement pour le champ 
quantique s'ecrivent (nous travaillerons dans le systeme de coordonnees de Minkowski 
tout au long de cette section : g^^ = diag(l, —1, —1, —1)) 



(2.1) 
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En remplagant le terme cfP' de la partie non-lineaire par sa valeur moyenne sur I'etat con- 
sidere, et en decomposant le champ en la somme de sa valeur moyenne 4> = {cj>) et des 
fluctuations autour de celle-ci : = + Scj), on obtient {{Scf)) = 0) 

(a^ + A(02 + {S^^) - v^)) (0 + Scf}) = Hn^ . 

La valeur moyenne de cette equation nous donne une equation pour la valeur moyenne du 
champ. En soustrayant cette derniere a I'equation ci-dessus, on obtient I'equation pour les 
fluctuations : 

{d^ + xixMix) = Hn^, (2.2) 
(d^ + x{x))S(t>{x) = 0, (2.3) 

ou 

x{x)=X{cl>\x) + {5cp\x))-v^) . (2.4) 

On voit bien ce qui se passe sur ces equations : le condensat (ou parametre d'ordre) 4> 
evolue dans un potentiel effectif qui n'est autre que le potentiel classique, corrige par I'effet 
des fluctuations quantiques (le terme {ScfP')), la dynamique desquelles est gouvernee par une 
equation de type Klein-Gordon avec une masse qui depend du temps. Les fluctuations du 
champ sont done effectivement decouplees et evoluent dans un potentiel efl'ectif quadratique 
dont la courbure est determinee d'une part par la valeur du condensat, et d'autre part par 
I'effet moyen des fluctuations elles-meme. La masse effective xit) represente la courbure 
locale du potentiel effectif a I'instant t. 

II existe different types d'approximations de champ moyen. Dans ce chapitre, nous 
travaillerons avec I'approximation dite "grand N" qui consiste a ne retenir que I'ordre 
dominant dans un developpement en ou N est le nombre de composantes du champ 
(A^ = 4). Les equations ainsi obtenues sont precisement les Eqs. (^^)-( p^ ). Bien entendu la 
derivation presentee ci-dessus ne constitue pas une preuve de cette affirmation et est seule- 
ment un moyen rapide d'ecrire les equations. II existe plusieurs manieres d'obtenir I'ordre 
dominant en (voir par exemple |5^, 53, 54, 55, Q). Ces differentes approches, bien 



qu' equivalent es, sont toutes aussi instructives les unes que les autres car elles eclairent le 
probleme sous des angles differents et mettent ainsi en lumiere divers aspects de I'approx- 
imation grand N en particulier, mais aussi, et plus generalement, de I'approximation de 
champ moyen. Dans la derivation ci-dessus, nous avons sacrifie le probleme de la lumiere 
a celui du plus court chemin, et cela pour plus de clarte. Nous procederons suivant cette 
ligne tout au long de ce chapitre, de fagon a ne pas alourdir notre discours. Cependant, 
certains complements utiles a la comprehension sont detailles dans les Annexes, auxquelles 
le lecteur est renvoye le cas echeant. 

Ici, nous partons des Eqs. (^^)- (|2.4| ) et ne faisons que mentionner les aspects de I'ap- 
proximation de champ moyen pertinents pour notre etude. Tout d'abord, et il est facile de 
le voir a partir des equations ci-dessus, il s'agit d'une approximation non-perturbative, dans 
ce sens qu'elle correspond a la resommation d'une classe infinie de termes de la serie pertur- 
bativeQ. C'est de plus une approximation semi-classique, ce dont on pent se convaincre si 
on reflechit a ce qu'est le champ cree par I'ensemble des particules du milieu. De fagon plus 
formelle, plagons-nous un instant dans I'image de Schrodinger : c'est maintenant la matrice 

^Dans le cas present, il s'agit des diagrammes dits "tadpoles" ou encore "daisy" et "super-daisy". Un 
autre exemple est donne par la resommation des boucles thermiques dures (HTL) dans les theories de 
jauge a haute temperature, ou le champ moyen decrit les effets collectifs de grande longueur d'onde (voir 
par exemple (57)). 
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densite du systeme qui depend du temps. L 'approximation de champ moyen correspond 
a imposer que celle-ci soit une gaussienne a chaque instant (voir par exemple ||3l[])0. Par 
exemple, si le systeme est initialement dans un etat coherent (c'est a dire semi-classique), 
il reste dans un etat coherent a chaque instant : dans 1' approximation de champ moyen, 
la dynamique est essentiellement classique. Ces deux aspects des Eqs. (p.2D -(^) sont tout 
a fait pertinents pour le probleme qui nous interesse : la formation d'un champ classique 



dans une theorie de couplage fort. Remarquons enfin que dans I'Eq. (|2.3D , les fluctuations 
dans les differentes directions d'isospin sont traitees sur un pied d'egalite. Get artefact de 
r approximation nous simplifie la tache. En effet, dans ce chapitre, nous nous interessons 
a I'intensite totale du champ de pion, c'est a dire sommee sur les directions d'isospin^. Le 
champ dcj) doit etre vu comme decrivant la fluctuation moyenne dans I'espace d'isospin. 

2.1.2 Expansion spherique 

Pour tenir compte de I'expansion du systeme, nous nous plagons dans un petit volume 
en co-mouvement, c'est a dire qui se deplace avec le "fluide" en expansion. Dans la suite, 
nous considerons le cas d'une expansion tri-dimensionnelle a symetrie spherique a la vitesse 
de la lumiere. Le temps mesure par un observateur dans un co- volume est le temps propre 
T = — r^, oil r est la distance a I'origine. La symetrie du probleme est telle que la valeur 
du condensat ne depend que de r. Vue dans le referentiel ou notre bulle en expansion est 
au repos, I'hypersurface r = cte est une sphere dont le rayon croit tres vite : a I'instant t, 
celui-ci vaut r = \Jt'^ — t'^ . En passant d'une hypersurface de temps propre constant a une 
autre, on decrit les differentes couches spheriques a I'interieur de notre bulle, elles-meme 
en expansion. Imaginez un oignon dont les differentes couches s'etendent les unes derriere 



les autres, comme dans le shema de la Fig. ^ 



Commengons par introduire quelques notations et definitions. 

Coordonnees 

L'element de longueur infinitesimal s'ecrit 

ds^ = dt^ - dr^ + r2(de2 - sin^ dc^^) . 
Les coordonnees adaptees a la description de notre co- volume sont 

r = - r2 , r/ = ^ In ^^-^^^ , , , 

en termes desquelles l'element de longueur infinitesimal s'ecritQ 

ds^ = g^,^ dx^'dx'' = dr^ - {Kj dx'dx^) , (2.5) 

ou hij est la metrique tri-dimensionnelle decrivant I'hypersurface r = 1 (encore appellee 
pseudo-sphere de rayon unite) 

hij = diag (1 , sinh^ r/ , sinh^ r/ sin^ 9) , 

et {g et h sont les determinants de g^^ et hij respectivement) 

^J—g = \/h , = sinh^ r/ sin . (2-6) 



^Ceci est intimement lie au fait que les operateurs de creation et d'anihilation des excitations (ou quasi- 
particules) du systeme a un instant t sont relies aux operateurs a I'instant initial par une transformation 
de Bogoliubov (voir plus loin, voir aussi 1' Annexe^). 

^L'etude de I'intensite dans les directions d'isospin individuelles fera I'objet du Chap. ^. 

■*I1 s'agit d'une metrique de Robertson-Walker [pq, p9[. 
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t>t 



Fig. 2.1: Schema de la bulle de matiere chirale en expansion vue dans le referentiel ou elle est 
globalement au repos. A I'instant initial t = to, on a une buUc dc rayon Rq a I'equilibre thermique 
local. Celle-ci "explose" : aux instants ulterieurs t > t^, \a surface externe de la bulle s'etend a 
la Vitesse de la lumiere, sa surface interne decrit I'hyperboloide r = \/t'^ — = tq = Iq. Les 
hyperboloi'des successives t > tq decrivent I'interieur de la bulle : ce sont des couches successives, 
en expansion les unes derriere les autres. 
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Laplaciens 

Le Laplacien quadri-dimensionnel est defini par 



^6 



_9 " ' 



ou designe la derivee covariante et A^^^ est le Laplacien tri-dimensionnel sur la pseudo- 
sphere unite 

Afin d'exploiter la symetrie spherique du probleme, il est utile de decomposer le champ sur 
I'ensemble complet des fonctions propres de A*^^), definies par^ ||60|] 

A(3)3;,-.(f) = -(s2 + i)3;.(f), (2.7) 

ou X = (r/, 6, if), s = (s, l,m) ; s est une variable continue, sans dimension < s < +00, / et 
m sont des nombres entiers < / < +00 et —I < m < I. Les "harmoniques hyperboliques" 
sont bien connues |61], elles forment une base orthogonale de I'espace des fonctions sur 
la pseudo-sphere unite. Nous choisirons leur normalisation comme suit 

d^x ^y*g,{x) yg{x) = 5(3) (s-- s') = 6is - s') 6iv 5mm' ■ (2.8) 

Dans la suite nous utiliserons les proprietes suivantes ||3^, 



yz = i-iry_, , \ys{x)\' = ^ , 



Im 



(2.10) 



avec les notations 



d^t 



+00 



ds 



E 

Im 



{s,l, 



-m) 



Notre observateur en co-mouvement mesure des observables physiques, comme par ex- 
emple le nombre de particules, a des "instants de temps propre" r. Pour decrire son systeme, 
il lui faut done quantifier la theorie sur les hypersurfaces r = cte. Le formalisme adapte a 



ce probleme est celui de la theorie quantique des champs en espace courbe |58|. Nous allons 
suivre ici les Refs. |33, ^ (voir aussi |S1]). 



Quantification 

Les equations dynamiques a I'ordre dominant en 1/A^ dans la geometric en expansion 
avec cf){x) = 4>{t) s'ecrivent 

0(r) + ^0(r) + x(t) (/)(t) = Hrv,, (2.11) 

(AW + x(r)) 50(r,x) = 0, (2.12) 

®0n peut montrer que ys{x) — Uai{ri) Yl^{9, ip), ou Yi"^{6, ip) sont les harmoniques spheriques usuelles, 
et Uaiirj) forment un ensemble complet de fonctions reelles, regulieres a I'origine. 
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ou le point denote la derivee temporelle, et ou 

x(r) = A(02(r) + (<502)(^)_^2) _ (2.13) 

Les differentes composantes chirales 5<j)a de la fluctuation etant traitees sur un pied d'e- 
galite dans cette approximation, nous concentrons notre attention sur une d'entre elles et 
omettons les indices chiraux pour simplifier les notations. Nous mentionnerons au fur et a 
mesure les modifications (evidentes) de nos formules quand ces indices sont reintroduits. 
Tout d'abord, il est judicieux d'introduire les variables adimensionnees 

if{T,x) = T5(j){T,x) , 7r(r, x) = r (/9(r, x) . (2-14) 

En effet, dans la suite nous parlerons des representations de Schrodinger et de Heisenberg, et 
nous aurons besoin de definir I'operateur unitaire qui connecte ces deux representations. Or, 



on pent montrer (voir Annexe ^ qu'en termes des variables (2.14) le probleme s'exprime 
comme un ensemble (infini) d'oscillateurs harmoniques dont les frequences dependent du 
temps (voir Eq. ( |2.22| ) plus bas). En definissant la representation de Schrodinger a I'instant 
de reference tq, et en denotant par U{t,tq) I'operateur d'evolution, on a 

(^(r) = U{T,To)ip{To)U-\T,To), (2.15) 

7r(r) = U{t,to)tt{to)U-\t,to), (2.16) 
ou Ton a omis la dependance spatiale. De la on tire les relations 

<50(r) = ^U{t,to)6<P{to)U-\t,to), 

T 

Hir) = ^U{t,to)6(I){to)U'\t,to). 



On peut construire explicitement^ I'operateur d'evolution U pour un potentiel quadratique 
general |^. Cependant, pour notre propos, nous n'avons pas besoin de cette construction, 
seules les lois de transformation ci-dessus nous seront utiles. 
La decomposition du champ sur les modes s s'ecrit 

^(r,f) = I dh (^agMT)yg{x)+alrs{r)ym) ■ 
Definissons les projections (fg et ng du champ et de son moment conjugueQ 

ifgir) = j d'x^yt{x)^{T,x) = ^s{T)ag + rs{r){-lTal^, (2.17) 

ng{T) = [d^xVhy;{x)^'{T,x)=ij',{T)ag + rs'{r){-iral., (2.18) 



I I s'a git d'une construction formelle faisant intervenir les solutions de I'equation du mouvement (cf. 



Eq. ( 2.21 ). II est toutefois interessant de noter que Ton peut chois ir arb itrairement les conditions initiales. 



pourvu que celles-ci satisfassent la contrainte de Wronskien (Eq. ( 2.23 ). 

'^Les fonctions modes, ainsi que les operateurs de creation et d'anihilation portent un indice chiral : tps, 
et Ua.s- Les relations de commutation sont : 
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ou le prime designe la derivation par rapport k u = ln(m7rr) : ip' = Tip. Les operateurs de 
creation et d'anihilation, at et a? satisfont aux relations de commutation 



,5(3)(s--s') 



(2.19) 



lea autres commutateurs etant nuls. Les "fonctions mode" ips satisfont aux equations dif- 
ferentielles 



que I'on pent reecrire 



T 



iPUr)+u^Ar)Mr)=0., 



0, 



(2.20) 



(2.21) 



ou CJo 



to;, et 



^s(t) 



+ x(t). 



(2.22) 



On voit que le systeme se reduit a un ensemble d'oscillateurs Eq.( p.21| ), couples entre eux 
d'une part, et avec le condensat d'autre part, a travers le terme de masse xi^)- L'effet 
de ce couplage se traduit uniquement par la dependance en temps de cette masse. Les 
frequences physiques sont donnees par ( ^.22]) (voir Annexe |A|) . Le champ de fluctuation 
6(p et son moment conjugue 11 = ^/—g 5(p (voir Annexe ^ satisfont aux relations de 
commutation canoniques a temps propres egaux [5(p{T,x) ; n(r, x')] = i 5^'^\x — x'). Celles- 
ci sont equivalentes a ( ^.19|) si les fonctions mode satisfont la condition de Wronskien 

Ws = IPs (r) V:' (r) - (r) iP', {t) = i . (2.23) 

Pour completer la procedure de quantification, il faut specifier les conditions initiales 
pour les fonctions modes, ainsi que pour leur derivees. En effet, differents choix pour les 
Ips correspondent a differentes definitions des quasi-particules, ou, de fagon equivalente, a 
differentes definitions du vide. Nous suivons ici les auteurs de la Ref. |^, en choisissant le 
vide adiabatique d'ordre zero ||6^, ^ (voir Annexe c'est a dire 

V's(to) = gs{ro) , i^'si'To) = g'siro) , 



avec 



9s{t) 



exp 



\/2tDs(T, 

II est commode pour la suite d'introduire la notation 

1 



v/2cDs(r) ' 
"^s'(r) 



2Cbs{T) 

La condition initiale pour les fonctions mode s'ecrit alors 



(2.24) 

(2.25) 
(2.26) 

(2.27) 
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L 'approximation adiabatique peut etre vue comme une definition du vide physique. Le 
probleme de la definition d'un "bon" vide vient de la dependance en temps de la metrique, 
Eq. ( |2.5D . Cette derniere, couplee au champ quantique par I'intermediaire du terme \/—g 
dans Taction, agit comme un champ classique qui rayonne des particules. Un etat initale- 
ment vide, se remplit de quanta au fur et a mesure que le temps passe. Ce probleme est 
brievement discute dans I'Annexe Pour resumer, I'approximation adiabatique consiste 
a minimiser le taux de production de particules du a I'expansion. Le choix des conditions 
initiales ( |2.27| ) n'est cependant pas essentiel pour notre probleme, ce que nous avons verifie 
en essayant differents choix. Mentionnons toutefois que la condition adiabatique permet de 
renormaliser la theorie de maniere non-ambigiie, independante du temps. 



Renormalisation 



L'equation de gap ( p. 13 ) fait intervenir la valeur moyenne du produit de deux operateurs 
au meme point ((5(/>^(r, x)). Cette quantite diverge et il est necessaire de regulariser la 
theorie, par exemple en introduisant une coupure ultraviolette A. II est bien connu que ce 
type de modele devient trivial (la constante de couplage tend vers zero) quand A — > oo, on 
doit done garder A fini. Comme nous I'avons discute plus haut, le modele o"-lineaire est une 
approximation de la dynamique des excitations bosoniques de grande longueur d'onde de 
QCD. C'est un modele effectif de basse energie qui n'a de sens physique que s'il est defini 
avec une coupure finie A < 1 GeV, qui doit etre vue comme un parametre. Cependant, il est 
utile d'isoler les divergences et de les reabsorber dans la definition de quantites physiques 
de fagon a minimiser la dependance des resultats avec A. C'est dans ce sens qu'il faut 
comprendre la "renormalisation" decrite ci-dessous. 

Afin d'isoler et de traiter les divergences de l'equation de gap, il est necessaire de 
specifier I'etat du systeme sur lequel est prise la valeur moyenne. Anticipant sur la suite, 
nous ecrivons 

ou I'on a reintroduit les indices chiraux. Nous supposerons un etat d'equilibre thermique 
local a I'instant initial, de sorte que decroit exponentiellement vite pour les grandes 
valeurs de s. En utilisant les Eqs. ( ^.14 ) et ( p.l7| ), ainsi que les proprietes ( ^.10 ) et les 



relations de commutation ( p. 191) , on obtient 

^ - "+00 9, 



{H\t,x)) = Y,^{6^l{T,x)) = -^j^ {2ns + 1) \Mr)\ 



2 



La partie divergente de cette integrale vient de la contribution du vide ~ / s'^ds \'ips\'^- Pour 
^ 1, on ecrit, en ligne avec I'approximation adiabatique^, 

■0s (r) = Qsir) + reste. 



8t 1 



L'approximation adiabatique est exacte pour s +oo, en effet, dans ce cas on a (Ds — s, et la solution 



de I'Eq. (2.21) ayant pour condition initiale (2.27) s'ecrit 



ip3 ^ — ^ exp[-zsln(r/r())] = gs{r) . 
V2s 

Remarquons que si le systeme, au cours de son evolution temporelle, traverse une periode d'instabilite 
(x < 0), certaines frequences ont des valeurs negatives pour cette periode, et J^^ ljs n'est pas un nombre 
reel. Pour ces modes, I'ans atz a diabatique ne peut etre utilise en tant qu'approximation de la fonction 
mode car alors, la relation ( ^.23| ) n'est pas satisfaite. Cette objection ne s'applique cependant pas pour les 
modes de haute energie qui nous interessent ici, et dont les frequences sont toujours reelles. 
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ou Qs est donnee par I'Eq. (2.24), et ou le "reste" ne donne pas de contribution singuliere. 



Ceci nous permet d'isoler la partie divergent 

2 / s^ds\gs{T)\^ - 2 



A h2 



2^JWTW) 

Cette derniere integrale est calculee dans I'Annexe on a 

OU ... represente des termes dependant faiblement de A. Ceci est valable quel que soit le 
signe de x- La dependance quadratique dans la coupure est independante du temps. Elle 
est facilement eliminee en soustrayant a la masse effective au carre x sa valeur dans le 
vide, qui n'est autre que le carre de la masse physique du pion. En effet, dans le vide on a 
(^(r) = (/7r,0) et I'equation de gap s'ecrit 

avec 

,2..v f^k'^dk 1 



(0 1 (50^10) = N 







2^2 2^A;2 + ml 

En exprimant x fonction de m^, ce qui revient a absorber la divergence dans le parametre 
V = f (A), on s'affranchit du terme quadratique en A : 

i(x(r) - ml) = <P\t) -fl + { 5(t>' )(r) - ( | 5<t>'' | 0) . 

La dependance logarithmique restante est eliminee en definissant la constante de couplage 
renormalisee \r, telle que (voir Annexe ^) ||5^, ^ 

1 I N k'^dk 

-r + 



Ar a 87r2 (fc2 + ml )3/2 

1 N ^ k 

Comme promis, la procedure de renormalisation est independante du temps. L'equation 
de gap est entierement exprimee en terme de quantites physiques mesurables, et done 
independantes de la coupure A, qui doivent etre tirees de I'experience. Nous utiliserons les 
valeurs proposees dans les Refs. |3|, HQ de fagon a pouvoir verifier nos calculs en les 
comparant aux leurs 

m^ = 139.5 MeV , 

= 92.5 MeV, 

\r = 7.3. 

^La coupure A s'interprete comme I'impulsioii maximale permise. Ici I'impulsioii est s/r (cf. Eq. ( 2.22 )), 
Tintegrale sur s est done coupee a Smax = Ar. 

^"Dans ces Refs. (voir aussi [^), la constante de couplage renormalisee \r est extraite de I'amplitude 
de diffusion tttt -ktt dans I'onde s (moment angulaire 1 = 0), dans la voie isoscalaire (isospin total 7 = 0). 
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Nous concluons cette partie par une remarque concernant les valeurs possibles de la 
coupure. La constante de couplage renormalisee depend de I'echelle d'energie fi a laquelle 
on la mesure, ici nous avons choisi /i = m^. Reecrivons la relation entre la constante de 
couplage nue A(A) et la constante de couplage renormalisee Xnifi) sous la forme 

Am = ^^(^) 

ou /3(A, /x) = l3(A//x) est I'integrale definie plus haut (voir aussi Annexe]^). On voit que 
pour Xr > il existe une valeur critique finie Ac de la coupure pour laquelle la constante 
de couplage nue diverge. En supposant que Ac/m.,^ » 1, on obtient I3 ~ ln(2Ac/m7r) — 1, 
c'est a dire {N = 4, XR{mTr) = 7.3) 

^ ~ 1 gi+87rVAfAfl _ 20 . 
2 

C'est le pole de Landau de la theorie : pour A > Ac, la constante de couplage nue 
A(A) est negative, la theorie est instable. On doit done prendre une valeur de la coupure 
inferieure au pole de Landau. La valeur numerique de Ac est de I'ordre de 3 GeV, ce qui 
est bien au dela de I'echelle typique d'energie que nous cherchons a decrire ici (< 1 GeV). 



Toujours en suivant les Ref. [|^, nous prendrons 

A = 800 MeV . 



2.1.3 Particules physiques, champ interpolant 

Notre but est de calculer le nombre de pions produits dans I'etat final. Pour ce faire, il 
nous faut definir precisement ce qu'est cette quantite. C'est le sujet de cette partie. 

Nous avons choisi, comme instant de reference, ou les representations de Schrodinger et 
de Heisenberg coincident, I'instant initial r = tq. Dans la representation de Schrodinger, 
on a la decomposition suivante 

ipg = ips{ro)=Mro)ag + iPt{To){-iraK, (2.28) 
vr.- = Mro)=<{To)ag+i^:'{To){-iral^. (2.29) 

Avec la condition adiabatique ( p.27| ), at (ag) s'interprete comme la representation de 
Schrodinger de I'operateur qui cree (anihile) un quanta de frequence uJs{to). Ce sont les 
quantas appropries a la description de I'etat initial ou le systeme est un plasma de pi- 
ons interagissant fortement les uns avec les autres. Dans I'approximation de champ moyen 
ce systeme en interaction est remplace par un ensemble d'excitations effectives libres de 
masse ^/xo, en d'autres termes on "diagonalise" le probleme en identifiant les excitations 
ou modes physiques du plasma. L 'information concernant I'interaction entre les pions est 
entierement contenue dans la masse effective de ces quasi-particules. II est bien clair que, 
du fait de I'expansion la definition des quasi-particules change avec le temps : les excita- 
tions physiques du systeme a I'instant tq ne sont plus les excitations physiques (c'est a dire 
qu'elles ne diagonalisent plus le hamiltonien) a I'instant tq + 5t. En fait les operateurs de 
creation et d'anihilation correspondant aux quasi-particules a I'instant r sont relies a ceux 
correspondant aux quasi-particules a I'instant tq par le biais d'une transformation unitaire 
dite de Bogoliubov (cf. Eq. ( |A.13D dans 1' Annexe]^). 



"^Si A_B < 0, A < et la theorie est instable 
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II existe une infinite de maniere de decomposer le champ sous la forme ( 2.28| )-( 2.29 ), 



chacune correspondant a un choix particulier d'operateurs de creation et d'anihiliation ou, 
de fagon equivalente, a un choix de la fonction mode et de sa derivee (la seule contrainte 
sur ces fonctions est qu'elle doivent satisfaire la relation de Wronskien ( p. 23 )). Les quasi- 



particules physiques du systeme a I'instant final r/ sont des quanta de frequence uJs{Tf). 
En denotant par ftt et hg les operateur de creation et d'anihilation correspondant, dans la 
representation de Schrodinger, on a 

= 6(ro)&s+C(To)(-ir&L-, 
= aro)hs + es{ro){-irb\ 



-s 



Les fonctions modes pour des quanta de frequence iOs{Tf) sont donnees par la condition 
adiabatique (|2^ )- (|2.26D 



6(to) = ^f(r;), 
e(To) = 4')(r/). 



Dans la representation de Heisenberg, les operateurs de champ ( 2.17] ) et ( |2.18| ) peuvent se 
reecrire (cf. (|l|)-(|2l|)) 

^s{t) = Uro)h,[T) + C{ro){-irhl.{T) , (2.30) 
vrs(r) = arMT)+e:{r,){-irhl^{T) , (2.31) 

ou hgij) = U{t,to) bgU~^{T,To) est la representation de Heisenberg de bg. En utilisant 
les Eqs. (|2.17| )- (|2.18D , ( 2.30 )-( ^3l|) , et (|2.23| ), on trouve aisement la transformation de 



Bogoliubov connectant les operateurs ag, at et bg{T), fe^(r) 

bg{T) = a,(r) ag + /3,(r) (-l)^*" al^, 

avec 

ia*(r) = ^s{to)'4'*s' (t) - (,'^{To)i;*{T) , 
mtir) = 6(To)V;(r)-C;(ro)^.(T). 

On verifie, a I'aide de ( |2.23| ), que cette transformation est unitaire a chaque instant]^ : 

K(r)p-|/3,(r)p = l. 

Bien sur, bien que les formules ci-dessus soient valables pour tout r, elles n'ont de significa- 
tion physique que pour r = t/. En effet, les quanta crees par I'operateur 6^ ne correspondent 
aux excitations physiques du systeme qu'a cet instant (ce sont des excitations de frequence 



^^La definition des quasi-particules que nous avons adoptee ici differe de celle de la Ref. [S^], ou la base 
dite "adiabatique" est utilisee de maniere abusive. En effet, cette base n'a de sens que si les frequences de 
tons les modes sont toujours reelles j63[, ce qui n'est pas le cas ici. Quand x < 0, la transformation reliant 
les operateurs de la base adiabatique a differents instants n'est pas unitaire pour les modes < r^Xj 
puisque le facteur exponentiel entrant dans la definition de la base adiabatique n'est pas une phase pure. 
Nous comprenons que, dans leur calcul, le facteur d'amplification a ete obtenu a partir de la formule derivee 
pour des frequences reelles (le facteur de phase n'apparait pas dans la formule du facteur d'amplification). 
Avec la derivatio n pre sentee ici, nous obtenons la meme formule pour le facteur d'amplification dans le cas 



general (cf. Eq. (2.37) 
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uJsiTf)). Avec les notations introduites ici, on pent ecrire une formule valable pour toute 
valeur de rj : 

bs{Tf) = a,{Tf)ag + (3s{Tf)i-iral., (2.32) 



avec 



ialiTf) = 4'\v)i;:'{Tf)-4'Hrf)rs{rf), (2.33) 
i(3:{Tf) = i;i^\rf)i^iiTf)-4''\Tf)Mrf). (2.34) 

Le nombre de quasi-particules physiques a I'instant Tf est facilement obtenu. Nous 
verrons, dans la prochaine partie, que dans I'etat initial 

{ai^s'^b,s') = K^ab^^^Hs- s') , {aa,sab,s') = 0, 

ou Ton a retabli les indices chiraux. On obtient alors 

{blgiy)bb,s'{rf))=n'',{Tf)5ab5'''\s-s') , (M^/) = , (2.35) 

et le nombre moyen de quasi-particules de composante chirale a dans I'etat final est donne 
par 

Kiy) + ^ = My) (n^s + l) , (2.36) 

ou Ton a defini le facteur d'amplification 

As{Tf) = 2\/3s{Tf)f + l (2.37) 

Nous verrons dans la suite que pour tj > 10 fm, le systeme a atteint un regime stationnaire 
ou xi^f) — f^'i '■ les excitations physiques sont des pions libres, dont nous calculerons la 
multiplicite a I'aide de la formule ( 2.36| ) ci-dessus. 



2.2 Echantillonnage des conditions initiales 

Cette partie constitue le coeur de notre travail. Nous considerons une petite bulle 
spherique de matiere chirale chaude, de rayon initial Rq^ que nous supposons avoir ete 
formee a un instant Iq dans la region centrale d'une collision nucleaire a haute energie, 
et qui subit une expansion spherique rapide. Nous supposons de plus que dans ce volume 
initial Vq, le champ de pion fluctue dans I'ensemble thermique local. Les effets combines 
des fluctuations initiales et de I'expansion rapide du systeme peuvent generer un "acciden- 
t" : une periode d'instabilite resultant dans I'amplification importante des modes de basse 
frequence. Notre but est d'estimer la frequence avec laquelle cet accident arrive, et plus 
precisement la probabilite d'avoir une amplification donnee. Le fait que cet accident ait 
lieu ou non depend de la configuration initiale du champ de pion a I'interieur de la bulle. 
II nous faut done echantillonner les configurations initiales possibles. 

Le modele presente dans les parties precedentes decrit la dynamique du champ de pion 
a I'interieur de la bulle, c'est a dire dans le quadri- volume r > tq = to. Notre modelisation 
de I'expansion spherique suppose qu'a chaque instant t, le parametre d'ordre est constant 
sur toute la surface de la sphere. Ceci est assure par le fait que la valeur du parametre 
d'ordre sur cette surface est la meme que ce qu'elle etait dans le petit volume initial, quand 
les differents points de la surface etaient en contact causal. Les conditions "initiales" sur 
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1' hyper surf ace r = tq, dont nous avons besoin pour resoudre les equations du mouvement, 
sont done donnees par la configuration du champ de pion dans le volume Vq au temps Iq. 

Remarquons que I'idealisation adoptee ici revient a supposer que la surface de la bulle 
est indefiniment couplee a un bain thermique a la temperature T. Nous discuterons plus 
loin les implications du fait que cette source d'energie exterieure ne dure qu'un temps fini 
dans une situation plus realiste. 



2.2.1 Les conditions initiales 

Comment est caracterisee une configuration initiale du champ de pion ? II est facile de 



se convaincre (cf. Eq. ( |2.12D , voir aussi la Ref. ^^), que 1' approximation de champ moyen 
correspond a une matrice densite p(ro) = po gaussienne et qui plus est, diagonale dans 
I'espace des impulsions. L'etat du systeme est done completement specifie par la donnee 
des valeurs moyennes et "largeurs" de cette matrice densite : 

et 

{al,a,^,,) = N^6ab6^'Hs- s') , {-ir (a„,_.-a,,,-.,) = p- 6ab d^^\s - s') , 

oil a,b = 1, N [N = 4). En principe les valeurs de A''^ et fluctuent d'un evenement a 
I'autre et doivent etre tirees aleatoirement dans I'ensemble thermique local. Cependant, en 
ce qui concerne revolution temporelle du systeme, ces quantites ne jouent un role que dans 
la formule de la masse effective ( 2.13D , ou elles apparaissent sous une integrale sur tous les 



modes. Les fluctuations statistiques sont alors moyennees a zero et on pent remplacer ces 
quantites par leurs valeurs moyennes 

ou Us est la distribution de Bose- Einstein pour les quanta de frequence iOs{To) et est in- 
dependant de la direction d'isospin dans I'approximation grand : 

1 

no 



.i^s{T0)/T _ I 



Les quantites et P^ interviennent aussi dans le calcul). Ici, nous sommes interesses a la 
multiplicite moyenne et nous remplacerons partout A'^^ et P° par leurs valeurs moyennes. 

Les quanta de l'etat initial sont definis par la donnee des fonctions modes ipsiTo) et 
de leur derivees ^l^s{To), qui sont a leur tour completement determinees par les valeurs de 
la masse effective initiale x{to) = Xo et de sa derivee xi'^'o) = Xo- On obtient ces deux 
quantites a partir de I'equation de gap, en utilisant I'Eq. ( ^.241) (ci-dessous, on fait le 
changement de variable A; = s/r) : xo est la solution de 

Xo = A [cPl-v^ + NAiixo)] , (2.38) 
et Xo satisfait I'equation algebrique 



Xo 



J + -jA.ixo) 



2 ^0 • ^0 - 77— 



2^i(xo) + Xo^3(xo) 



A3 coth[^A(xo)/2T] 

27r2 ^a(xo) 

(2.39) 



avec cj>o = 4>{tq) et c/jq = 0(to), et oil Ton a introduit les notations suivantes : 



(2.40) 
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A u2 



1 



n r„ / COth . 

27r2 [Ek{p?)Y V 



2T 



(2.41) 



II est facile de verifier que I'Eq. (|2.39D , exprimee en termes de la constante de couplage 
renormalisee Ar, ne contient aucune dependance forte dans la coupure A, comme il se doit. 

En conclusion, on voit que la configuration initiale du systeme est completement spe- 
cifiee par la donnee des huit nombres reels {(fiQ , 4>o}- Un exemple en est donne dans la 
Fig. ^.2| , qui represente la variation temporelle de la masse effective au carre x{'t) pour un 
des choix des valeurs initiales du parametre d'ordre et de sa derivee faits dans |Q. Dans 
cet exemple x prend des valeurs negatives et on assiste a une forte amplification des modes 



de basse frequence. La Fig. 2^ montre le facteur d'amplification As{Tf) correspondant (cf. 
Eq. ( p. 37 )) en fonction de s pour r/ = 5, 10 et 15 fm. 



1 




-0.5 

5 10 15 20 

T(fm) 



Fig. 2.2: Evolution temporelle de xi^) (en fm~^) pour la consition initiale (cf. Q) ^(tq) = 
(0.3,0) (en fm~^), </)(to) = (-1,0) (en fm~^), a tq = 1 fm. Les parametres sont A = 800 MeV, 
= 139.5 MeV, /tt — 92.5 MeV, \b. — 7.3. La ligne en tirets montre la valeur de rn^ — 0.5 fm^^. 

Dans la litterature, les valeurs des conditions initiales sont generalement choisies de 
maniere arbitraire. Dans la partie suivante, nous allons montrer que I'liypotliese, largement 
utilisee, d'un etat initial en equilibre thermique local est suffisante pour determiner la 
probabilite pour que ces nombres prennent des valeurs donnees. 

2.2.2 Equilibre thermique local 

Comme nous I'avons vu plus haut, les conditions sur I'hypersurface r = tq sont donnees 
par la configuration du champ a I'interieur de la bulle de rayon Rq, en equilibre thermique 
local a I'instant tQ. Dans cette partie, nous focalisons notre attention sur cette derniere. 
Nous utiliserons des notations conventionnelles, a ne pas confondre avec les notations util- 
isees dans le reste du chapitre : x designe la position dans le systeme de coordonnees 
cartesiennes, x = {x,y,z), d^x = dxdydz. Nous omettrons d'ecrire la variable temporelle, 
partout egale a to- 

Pour exposer notre argument de fagon claire, nous negligeons dans un premier temps les 
complications dues a la mecanique quantique et considerons I'exemple d'un champ scalaire 
classique, sans structure d'isospin, en equilibre thermique local a la temperature T dans 
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Fig. 2.3: Le facteur d' amplification As{Tf) correspondant a I'evenement de la Fig. |2j, pour 
differentes valeurs de t/. 



une petite boule de volume Vq. Cela signifie que tout se passe comme si cette boule etait 
une sous-partie d'un grand systeme de volume V ^ Vq en equilibre thermodynamique a la 
meme temperature T. Dans ce grand systeme fictif, en I'absence de correlations a longue 
distance (ce qui est le cas pour des interactions a courte portee) la variance des fluctuations 
statistiques de la moyenne spatiale du champ sur le volume total V est tres petite, d'ordre 
1/V (elle est nulle dans la limite thermodynamique V ^ oo, cf. Eq. ( 2.51| )). La variance 
correspondante pour la moyenne spatiale du champ sur le volume Vq est plus grande d'un 
facteur V/Vq (pourvu que le rayon de la boule soit au moins de I'ordre de la longueur 
de correlation thermique^ Xt = l/y/x^- Le meme raisonnement s' applique a la derivee 
temporelle du champ moyennee sur le volume Vq. Le point essentiel est que, dans le systeme 
en equilibre local, les valeurs du champ et de sa derivee, moyennes sur le volume Vq de 
la bulle, fluctuent dans I'ensemble canonique et done de maniere parfaitement predictible, 
pourvu que le rayon initial Rq ne soit pas trop petit. 

Appliquons ces idees a notre probleme. Dans la theorie quantique, on definit les opera- 



^^Considerons le systeme de volume V, et divisons le en = V/Vo sous-parties identiques de volume Vq. 
La moyenne spatiale du champ dans chacune de ces sous-parties est 



= — I d?x ( 
Vo 



(x) , i = l,2, 



Choisissons le rayon i?o des cellules au moins egal a la longueur de correlation du champ a I'equilibre 
thermodynamique a la temperature T, Ro^ \t- Ainsi, par construction, les <f>i sont des variables aleatoires 
independantes, ayant toutes la meme distribution de probabilite. La somme de ces variables aleatoires, 
S — X^iLi <^ij 6st une variable aleatoire dont la moyenne statistique est egale a la somme des moyennes 
statistiques des 4>i, et dont la variance est egale a la somme des variances des </)i : 'E[S] = fcE [(/>], et 
Var[S'] = feVar[(^]. La somme S n'etant autre que la moyenne spatiale du champ sur le volume total V, a 
un facteur k pres : 4>v ~ S/k, on en deduit 



Varfi 



V_ 

Vo 



Var[( 
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teurs (a designe les composantes chirales) 



1 

Vo 
1 



Vo 



Vo 



Ceux-ci ne regoivent pas de contribution de la part des modes de petite longueur d'onde, 
moyennes a zero. Un observateur vivant dans la boule de volume Vq, et mesurant ces 
observables, identifiera done le resultat de sa mesure avec les valeurs du parametre d'ordre 
et de sa derivee respectivement. La distribution des valeurs possibles (j)a et (p^ des resultats 
de cette mesure est donnee par la distribution de Wignei0 (voir Annexe caracterisee 
par les moyennes 



E 



les variancesf^ 



et par la covariance 



Gov 



Var [(pa 
Var <pa 

1 



I T 



(2.42) 
(2.43) 

(2.44) 
(2.45) 



(2.46) 



ou {0)j, = Tr{pO) designe la moyenne sur I'ensemble thermique (p oc exp(— iJ/T), ou 
H est le hamiltonien du systeme). Le calcul des quantites ci-dessus est detaille dans 1' An- 
nexe^. L'etat d'equilibre thermodynamique est homogene et invariant par translation dans 
le temps. En denotant par 4>^ la valeur du parametre d'ordre dans cet etat et Xt ceWe de 
la masse effective au carre, on obtient, a I'aide de I'Eq. ( 2.11|) , 



H 



d'ou on deduit 



E 
E 



rio- , et cf)^ = , 



0. 



L'equation de gap thermique, dont Xt ^^t la solution, s'ecrit (cf. Eq. ( 2.41] )) 

A Ul - + N A^ix^)] , 



(2.47) 

(2.48) 
(2.49) 

(2.50) 



ou (j)^, est donne par I'Eq. ( 2.47| ). Dans I'approximation de champ moyen, le systeme en 
interaction est remplace par un ensemble d'excitations libres de masse effective ^/x^, c'est a 

^*La distribution de probabilite de la valeur mesuree d'une observable O est determinee par la fonction 
caracteristique (e*'''^). 

^^Tous les termes non-diagonnaux dans les indices chiraux sont nuls dans rapproximation grand A'^ 
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dire que le hamiltonien effectif est quadratique et la matrice densite de I'ensemble thermique 
p oc exjp{—H/T) est gaussienne. La distribution de probabilite des valeurs possibles de et 
(p est done completement determinee par les parametres ( 2.42D -( ^^ ), la covariance etant 
nulle. 

Quand le rayon Rq est grand devant la longueur de correlation dans I'ensemble ther- 
mique a temperature T, At = les fluctuations a I'interieur de la bulle sont indepen- 
dante de ce qui se passe a I'exterieur, et on pent donner une expression analytique pour 
les variances : 



Var [(pa 



1 



coth 



2T 



Var 



2Vo 



coth 



2T 



(2.51) 
(2.52) 



La dispersion de la variable (pa, calculee numeriquement (voir Annexe est plus petite 
d'un facteur 2 (1.5) que ce que donne la formule ci-dessus pour Rq = At (Rq = 2At), pour 
des temperatures allant de T = 200 a 400 MeV. Pour la variable i^^, I'ecart entre le resultat 
analytique ci-dessus et le calcul exact est de 20% (8%). Get ecart augmente rapidement 
pour Rq < At- Les fluctuations sont approximativement independantes de ce qui se passe 
hors de la bulle, pourvu que Rq ^ ^t- 

En resume, la configuration initiale du systeme est completement determinee par la 
donnee des quatres composantes chirales du parametre d'ordre et de sa derivee a I'instant 
To : 4>{to) et 0(to). Ces huit nombres sont des variables gaussiennes independantes dont 



les moyennes sont donnees par les Eqs. (|2.48|) et (|2.49|) , et dont les variances sont approx- 
imativement donnees par les Eqs. (|]5l|) et (|]5|). 



2.3 Resultats et discussion 

Dans cette partie nous appliquons la methode d'echantillonnage des conditions initiales 
decrite ci-dessus au calcul de la probabilite d'avoir une amplification donnee dans I'etat 
final. Nous resumons tout d'abord la strategic du calcul et examinons la question du choix 
des parametres. Ceux-ci sont contraints par la coherence physique du modele et des ap- 
proximations utilisees. De plus, ils doivent etre choisis en correspondance avec la question 
posee. Enfin, nous presentons et discutons les resultats du calcul, puis nous definissons les 
valeurs de I'amplification correspondant a un phenomene obsevable dont nous calculous la 
probabilite d'occurence. 

2.3.1 La strategic : recapitulatif 

Les parametres de notre modele, outre ceux deja fixes : m-n-, fn, et A, sont I'instant 
^0 = ■^0, le rayon Rq et la temperature T initials. Nous discuterons plus loin le choix des 
valeurs de ces parametres. Commengons par enumerer les differentes etapes du calcul : 

1. Nous calculous d'abord la longueur de correlation At du champ dans I'ensemble 
thermodynamique a la temperature T. Celle-ci est donnee par I'inverse de la masse 
effective des excitations (quasi-particules) du systeme, dont le carre est solution de 
I'equation d'auto-coherence (Eqs. ( |2.47D et ( |2.50D ) a I'equilibre thermique : 

A \XtJ Jo 2E,{x^) 
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Puis nous calculous les moyennes et dispersions (cti et £72) des distributions de proba- 
bilite des valeurs possibles du parametre d'ordre 0(to) et de sa derivee </>(to) a I'aide 
des formules ( pT7|) , (pl8|) , (^]2|) et (^^251) de 1' Annexe g : 



T/2 ^2 



T/2 ^2 
^0 0^2 



ARl / £ coth ) ^(A:/2o) 



4i2^ / dkEk coth 



2r 



oii F{y) = (sin 2/ - y cosyf'/y^ et E'fc = 
2. Nous tirons aleatoirement les valeurs initiales des composantes chirales du parametre 
d'ordre et de sa derivee : |0(ro); </>(to)|, avec lesquelles nous initialisons complete- 
ment le systeme. Nous calculous tout d'abord la masse effective des excitations 
physiques du milieu a I'aide de I'equation d'auto-colierence a I'instant initial (equili- 
bre thermique local) ( 2.38 ) : 

Xo = A - + N Aiixo)] , 

et dont la derivee temporelle a I'instant initial est donnee par I'Eq. O). Avec ces 
deux nombres, nous initialisons I'ensemble des fonctions modes ipsiTo) (Eqs. ( p.26 )- 



2(2's(ro) 



+ iu)s{ro) 



avec 



Ws(to) = touJs{to) = v^s^Trfxo et a}^(To) = (2to Xo + Tq xo)/2uJs{to). 
3. On pent alors resoudre les equations du mouvement ( |2.11| ), ( |2.20 ) et ( 2.13| ) 



4>{t) + -(j){T) + xir) cf>{T) 

T 

i^s{T) + -^s{T)+Ujl{T)i>s{T) 



et 



X(r) 
A 



.2, ^ 2 Z"^" 



s^ds 
'2^ 



0, 

:2n, + l)|V'.(r)|^ 



4. On obtient ainsi les valeurs des fonctions modes et de leur derivees a I'instant final : 
V's(rj) et ips{Tf), ainsi que celles de la masse des excitations physiques et sa derivee : 
x(r/) = Xf et x{Tf ) = Xf, a I'aide desquelles on calcule la multiplicite finale ( p. 36 ) 



<(r/) + ^=A(r/) + 



ou le facteur d'amplification est donne par 

A(r/) = 2|/3,(T/)|2 + 1 

avec (|2]3|) 



mrf)\' = \4'Hrf)^'s{rf)-4'Hv)Mrf)\'. 



2.3. Resultats et discussion 



35 



et (|2]25|)-(|]26|) 



V2a;s(7"/) 



ou Ws(r/) = ^s'^ + t'I Xf et ^^(t/) = (2t/ Xf + r] Xf)/2ujs{Tf). 

5. On repete les pas |2| a ^ autant de fois que possible, de fagon a avoir une bonne 
statistique. 

6. Enfin, on recommence avec d'autres valeurs des parametres. 

2.3.2 Coherence physique du modele et choix des parametres 

II est clair que le modele effectif de basse energie que nous utilisons ici n'a de sens 
que pour des temperatures T <C A. De fagon plus fine, le modele cr-lineaire n'a pas de sens 
pour des temperatures trop elevees devant la temperature critique de la transition de phase 
chirale Tc ^ 160 MeV. En effet, a ces temperatures, on doit tenir compte des degres de 
libertes colores que sont les quarks et les gluons, la temperature critique de la transition de 
deconfinement etant du meme ordre de grandeur. Pour la physique qui nous interesse, on 
doit cependant demarrer dans la phase symetrique (T > Tc), et on neglige le phenomene 
du deconfinement, sans toutefois aller trop haut en temperature. 

Voyons d'abord quelles sont les limites imposees par la coherence physique du modele 
sur les parametres Rq et tq = to- En ce qui concerne la taille de la bulle initiale, nos 
hypotheses ne sont valables que pour Rq > At- En effet, I'hypothese d'equilibre local n'a 
de sens que si les degres de liberte du systeme fluctuent plus ou moins independamment de 
I'environnement de la bulle. Dit d'une autre fagon, nous avons vu que I'approximation de 
champ moyen consiste a voir le systeme comme un ensemble d'excitations independantes de 
masse effective ^Jxx) ~ \/X^- Cette image n'a de sens que si le volume dans lequel vivent ces 
excitations est au moins superieur a leur extension spatiale ~ 1/ ^fXr- Nous n'excluons bien 
evidemment pas la possibilite que I'instabilite que nous cherchons a decrire se developpe 
dans une bulle plus petite, mais ce scenario est en dehors du domaine d'applicabilite de 
notre modele. L'instant initial to est le temps de formation de notre bulle lors de la collision. 
Celui-ci doit etre superieur au temps de formation d'une quasi-particule qui est au moins 
de I'ordre du temps caracteristique des interactions fortes ~ 1 fm. De meme, il ne fait 
aucun sens de parler d'equilibre pour des temps inferieurs a cette echelle. Enfin, on a 
to ^ ^0 ^ ~ 1 fni pour les temperatures considerees. 

Revenons maintenant a la question de depart : nous voulons estimer I'ordre de grandeur 
de la limite superieure]^ de la probabilite d'avoir une amplification notable. Nous nous 
plagons done dans les conditions les plus favorables pour qu'une instabilite se developpe. 
Nous avons deja vu que le trempage des fluctuations initiales est d'autant plus efficace que 
le terme de friction dans les equations du mouvement est grand. Celui-ci etant inversement 
proportionnel au temps, nous commencerons I'expansion le plus tot possible, c'est a dire 



Notre description repose implicitement sur I'hypothese que cette probabilite est rare. En effet, dans le 
cas contraire il faudrait s'inquieter de la possibilite que plusieurs bulles developpant des instabilites soient 
formfies lors de la collision, et prendre en compte les interactions entre des bulles voisines. En particulier, 
I'hypothese selon laquelle la valeur du parametre d'ordre est constant sur la surface de la bulle n'aurait 
alors aucun sens. 



36 



Chapitre 2. La probabilite de formation 
d'un champ de pion classique 



que nous prendronsQ to = Rq = Xt- Par souci de comparaison, nous etudierons aussi le 
scenario to = Rq = 2Xt- 

II pourrait sembler qu'une temperature initiale elevee (qui correspond a une agitation 
thermique plus grande, et done a une longueur de correlation plus faible) soit preferable. 
Cependant de trop fortes fluctuations initiales sont plus difficiles a supprimer par trempage. 
On peut le voir sur I'equation de gap (|2.13|) : il est d'autant plus facile d'entrainer le carre de 
la masse effective dans la region des valeurs negatives que la contribution des fluctuations 
) est petite (voir aussi la Fig. L3 dans le Chap. |l]). Dans la suite, nous considererons 
deux valeurs de la temperature initiale : T = 200 MeV et T = 400 MeV. 

Le dernier point concerne le temps rj auquel on mesure la multiplicite finale. Pour 
I'ensemble des evenements consideres, le systeme atteint un regime stationnaire ou x{t) ~ 

pour T > 10 fm. On a alors un gaz de pions libres en expansion. Nous prendrons 
Tf = 10 fm. 



2.3.3 Resultats 

Nous resolvons numeriquement les equations du mouvement en utilisant un algorithme 



de Runge-Kutta d'ordre quatre ||6j], les integrales sur les modes etant calculees par la 
methode de Simpson, corrigee par une methode de trapezes lorsque le nombre de points sous 
I'integrale est pair (rappelons que la borne d'integration depend du temps). Les equations 
de gap a I'equilibre thermodynamique et a I'instant initial sont resolues par une simple 
methode de bissection. Nous reproduisons les resultats des Refs. |33, |3^F^ (voir Fig. 2.2). 



Quand le systeme entre dans la region d'instabilite spinodale x < 0, les modes dont la 
frequence tOs devient imaginaire pure sont tres fortement amplifies. II est clair que I'am- 
plification est d'autant plus importante que s est petit : le facteur d'amplification est une 
fonction monotone decroissante de s. Dans la suite nous focalisons notre attention sur le 
mode s = O"*", le plus amplifie. Nous appelons P{A) la probabilite pour que le facteur 
d'amplification correspondant, calcule a tj = W fm, ait la valeur = ^- La Fig. 
represente les histogrammes des probabilites P{A) pour les quatres scenarios correspon- 
dant aux valeurs T = 200 et 400 MeV et Rq/Xt = 1 et 2. 

Comme on s'y attend (voir la discussion plus haut), le cas Ro = 2At est le plus defa- 
vorable en ce qui concerne la possibilite de fortes amplifications. Dans tous les cas, on voit 
clairement que seule une petite fraction des evenements correspond a une amplification 
importante. Les histogrammes correspondant au cas Ro = At sont assez bien representes 
par 

p(A) = 

(1 + 0.27^)3-21' 

pour T = 200 MeV, et 

PiA) 



(1 + 0.032^)6-959 ' 

pour T = 400 MeV. II est impossible de reproduire la queue des distributions par une loi 
exponentielle du type P{A) cx exp(— ayl''). 

II est interessant de chercher a caracteriser les etats initials destines a etre fortement 
amplifies. II semble qu'une condition soit la petitesse relative de I'intensite du courant iso- 
vectoriel classique, c'est a dire calcule avec le parametre d'ordre (p, a I'instant initial. Les 



^'^Remarquons de plus que pour que I'hypothese d'equilibre local a I'instant initial ait un sens, nous 
devons choisir -Ro aussi grand que possible, c'est a dire -Ro = to- 

^^Nous remercions M. A. Lampert, de nous avoir communique les resultats de ses calculs 
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Fig. 2.4: Le facteur d'amplification du "mode zero" dans les quatrcs scenarios suivant : T = 200 
MeV et Ro^ Xt = 1.17 fm (3.2 x 10* evenements) ; T = 400 MeV etRo = At = 0.68 fm (10* evts) ; 
T = 200 MeV et Ro = 2At (2.5 x 10* evts) ; T = 400 MeV et Ro = 2At (10* evts). 
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courants iso-vectoriel et iso-axial classiques sont definis par 

= TT X d^TV , 

ou Ton a ecrit cj) = (7r,cj). Les intensites initiales • et • de ces courants sont 
representees en fonction de 1' amplification finale A sur la Fig. |2]^ (dans le cas T = 200 MeV, 
Rq = At). II est interessant de remarquer que dans le modele classique en 1+1 dimensions 
(expansion longitudinale) des Refs. |l^, la petitesse du rapport des valeurs initiales 
des intensites des courants vectoriel et axial est la condition pour qu'aux temps longs 
le parametre d'ordre oscille dans une direction bien determinee de I'espace d'isospin (voir 
aussi le Chap. |l]). Autrement dit, c'est la condition pour avoir une configuration DCC. Dans 
le cas present cependant, le champ classique forme dans I'etat final est une superposition 
du parametre d'ordre (le mode zero) et des modes de grande longueur d'onde. II n'est pas 
clair que tous oscillent collectivement dans la meme direction de I'espace d'isospin. Cette 
question sera examinee dans le Chap. |3|. 

Nous voulons maintenant faire le lien avec la phenomenologie et determiner quelles 
sont les valeurs de I'amplification qui donnent lieu a phenomene (potentiellement) observ- 
able. Estimons la densite Edn/d^p de pions produits par unite invariante d'espace des 
phases. Ces pions sont produits par les differents elements de volume en co-mouvement, 
c'est a dire qu'ils sont emis sur I'hypersurface r = t/. Les pions emis par I'element de 
volume r| sinh^r/ sm.9 dr] dO d(p = r^Vhd^x, centre sur le point de coordonnees hyper- 
bolique {ri,6,ip) ont, dans le referentiel ou celui-ci est au repos, une energie ujs{Tf). Dans 
un referentiel quelconque, cet element de volume est en mouvement avec la quadri-vitesse 



u 



At 



dx^/dr, et la quadri-impulsion des pions emis est telle que p^Ufj_ = ujs{Tf). 
Cette relation nous dit quelles sont les valeurs de I'impulsion k = s/tj qui contribuent 
a la quadri-impulsion = {E,p) {p^ = Xf), compte tenu du mouvement de I'element 
de fluide. Par exemple, dans le referentiel ou la bulle est globalement au repos, on a 
p'^UfM = Ecoshi] — pcosO sinhr] = ujs{Tf). La contribution de cet element de volume a la 
densite invariante Edn/d^p est proportionnelle a ns{Tf) p^u^ d'^x \/—g 5{t — r/), ou Ton a 
ecrit I'element de volume sous une forme manifestement invariante. Au total, on a done ||65| 



^^(p) = j d'xV^6iT-Tf)fix,p)p^n^, (2.53) 

ouQ 



N 



L 'integrand dans ( 2.53| ) ne dependant que de la quadri-impulsion p^, Edn/d^p est une fonc- 



tion de I'invariant p^ = x = cte et est, par consequent, constante : on a un spectre plat^. 
Ceci est dii a I'idealisation que nous avons adopte, qui consiste a considerer que le systeme 
est couple a un bain thermique pendant une duree infinie ou, en d'autres termes, que I'ex- 
pansion dure un temps infini. Par exemple, dans le referentiel oil la bulle est globalement 
au repos, aussi grand que soit le module p de I'impulsion consideree, il existe toujours un 
ensemble d'elements de volume qui se deplacent suffisament vite (avec une grande rapidite 



^^La formule de Planck pour le rayonnement du corp noir est obtenue, a partir de I'Eq. ( 2.53 ), en 
remplagant la contrainte t — Tf par t = cte, de sorte que p^u^ = E, et en prenant J{x,p) = 
2(27r)~^[exp(i?/r) — le facteur 2 etant le nombre d'etats de polarisation possibles du photon. Cet 

exemple permet de fixer les normalisations. 

^"Le spectre decroissant obtenu dans les Refs. js^, est incorrect. Ceci a ete confirme en prive par 
F. Cooper. 
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Fig. 2.5: Les valeurs initiales des intensites des courants iso-vectoriel et iso-axial classiques cor- 
respondant aux amplifications observees dans I'etat final. On a calcule la valeur moyenne des 
intensites sur I'ensemble des cvcncmcnts dont I'amplification est comprise entre et 50, 50 et 
100, et 100 et 150. Les barres d'erreurs representent les incertitudes statistiques. En moyenne, les 
evenements ayant subit une forte amplification se distinguent des autres par la petitesse relative 
de I'intensite initiale du courant vectoriel. 
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radiale) et qui contribuent a la densite invariante. On peut modeliser grossierement I'effet 
de la duree finie du regime d'expansion en limitant I'integrale sur la variable de rapidite 



dans I'Eq. ( 2.53 ) aux valeurs ij < rjsup, oil rjsup reflete I'influence de I'environnement de la 
bulle. Dans ce cas, chaque element de volume en co-mouvement emettant essentiellement 
des pions de basse energie (dans son referentiel propre), le spectre (|2.53|) est coupe pour 
des valeur de I'impulsion psup ~ Xf smhr]sup- Une discussion de ces aspects necessite la 
modelisation de la collision dans son ensemble, ce qui depasse le cadre de notre etude. La 
multiplicite totale (irrealistement infinie dans notre cas) depend fortement de cette coupure, 
et ne peut etre estimee de fagon fiable dans ce modele. Cependant, la formule (|2.53| ) (voir 



aussi ( 2.54 ) plus loin) est probablement une estimation raisonnable de la densite invariante 
de pions de basse energie. On ecrit 



E 



N 

0) = — 

Xf 



^^f s^ds 



■ n 



:irf) 



Pour les evenements ayant subi une forte amplification, la multiplicite est approximative- 
ment proportionnelle au facteur d'amplification du mode le plus amplifie. Sur la Fig. p.6| 
on a represente la quantite c definie par : 



E 



dn 
d^ 



cA. 



(2.54) 



en fonction de A, pour un ensemble d'evenements tels que A > 20, dans le cas T = 200 MeV 
et Rq = At- Le nombre c fluctue relativement peu d'un evenement a I'autre pour A > 30. 
On observe la meme chose pour T = 400 MeV. On a c ~ 5 GeV~^ dans les deux cas. 
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Fig. 2.6: La constante c (Eq. ( 2.54 )) en fm^ pour quelques (268) evenements dans le scenario 
T — 200 MeV, Rq ~ Xt- A gauche, on a represente c en fonction de I'amplification A pour A > 20, 
a droite, I'histogramme correspondant, restreint aux evenements tels que A > 30. Pour ces derniers 
on observe < c >= 0.19 fm^, avec Ac — {< > — < c >2)i/2 = 8 x 10~^ fm^. On obtient un 
resultat similaire dans le cas T — 400 MeV. 

Cette estimation est a comparer avec la densite de pions mous produits par unite 
d'espace des phases lors d'une collision d'ions lourds. II est instructif de considerer I'exemple 
simple suivant : dans une collision typique, le spectre de pions emis est relativement plat 
dans la region centrale de rapidite y ~ 0, ou y = (1/2) ln[{E + pz) / (E —pz)] est la rapidite 
longitudinale (z est I'axe de la collision). Denotons par h la hauteur de ce plateau. Dans 
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la region centrale, on pent parametriser le spectre inclusif par 

, dn 



d^p 



dn 

ind dy d^pt 



Pt 

ind ^(P?) ''''^ V {pT, 



Dans les collisions centrales Pb-Pb au SPS du CERN, on observe ||66|| un nombre de I'ordre 
de 200 TT~ par unite de rapidite, c'est a dire, pour tons les pions, h = dn/dy 600. Done, 
la densite de pions de tres faible impulsion transverse ~ 0, vaut approximativement 
1900 GeV~^, ou on a pris {pf) =0.1 GeV^. Les fluctuations de cette densite sont ~ 
Vh/n{pf) w 75 GeV-2. Pour que notre signal soit observable, il doit se detacher de ces 
fluctuations statistiques. D'apres I'estimation ( |2.54| ), on voit que le signal est superieur a 
trois fois la fluctuation si A > 45. Dans le cas i?o = At, la probabilite correspondante est 
approximativement 

Proba(A > 45) 4 X 10"^ , (2.55) 

aussi bien pour T = 200 MeV que pour T = 400 MeV. Cette estimation nous donne I'ordre 
de grandeur de la limite superieure de la probabilite d'avoir une amplification observable. II 
s'agit, bien stir, d'une probabilite conditionnelle, puisque notre calcul repose sur Thypothese 
selon laquelle notre bulle est initialement dans un etat d'equilibre thermique local (ou, ce 
qui a un sens plus faible, que les fluctuations initiales sont correctement decrites par un 
ensemble thermique local). 



2.4 Conclusion et commentaire 

Resumons-nous : dans le cadre le plus courament utilise dans la litterature sur les DCC 
(modele cr-lineaire, approximation de champ moyen, equilibre local a I'instant initial), nous 
avons construit une methode originale d'echantillonnage des valeurs initiales du parametre 
d'ordre et de sa derivee temporelle, qui specifient completement I'etat initial du systeme. 
Nous plagant dans le cadre le plus favorable possible pour generer une forte amplifica- 
tion des modes de grande longueur d'onde (expansion spherique, taille initiale de la bulle 



aussi petite que possible), et combinant le formalisme des Refs. |33, 34, [G^] avec notre 
methode, nous avons calcule la distribution des valeurs possibles du facteur d'amplifica- 
tion (Eq. ( |2.37| )) du mode le plus amplifie. Une analyse phenomenologique simple nous 
permet alors d'estimer la probabilite d'avoir une amplification donnant lieu a un signal 
detectable (Eq. (^)). Celle -ci s'avere assez faible, grossierement de I'ordre de 10 pour 
une collision centrale Pb-Pb au SPS du CERN. Ce resultat doit encore etre multiplie par la 
probabilite pour que nos hypotheses soient verifiees dans une collision reelle. Le calcul de 
cette derniere necessite une description de I'ensemble de la collision, ce qui va bien au-dela 
de I'etude presentee ici. 

II est raisonnable de penser que notre estimation, meme si elle est tres approximative, 
donne I'ordre de grandeur auquel il faut s'attendre et doit done etre prise en compte aussi 
bien par les experiment at eurs chasseurs de DCC que par les theoriciens. En particulier, 
il ressort de notre etude que la probabilite pour que plusieurs bulles subissant une forte 
amplification soient formees lors d'une collision est faible^. De plus, du point de vue 
experimental, les analyses globales sont a proscrire, au benefice d'analyses "evenement-par- 
evenement". 
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Si plusieurs DCC sont formes durant une collision, leurs ori enta tions respectives etant independantes 



la distribution en charge du nombre total de pions emis (cf. Eq. (L3)) devient de plus en plus piquee autour 
de la valeur moyenne 1/3 a mesure que le nombre d'emetteurs augmente [jsTt - 
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Nous terminerons ce chapitre par un commentaire concernant I'attitude generale que 
nous avons adopte quant a sa presentation : nous avons calcule la probabihte d'avoir une 
amphfication donnee, ayant en tete le fait que ce phenomene correspond a la formation 
d'un DCC Cependant, une etude realisee ulterieurement (cf. Chap. |3|) indique que la 
situation n'est pas aussi claire. En effet, nous verrons que si le mecanisme du trempage est 
responsable de la formation d'un champ classique de pion - ce qui est une caracteristique 
essentielle de la configuration DCC - il ne suffit pas a generer la structure d'isospin du DCC. 
Le calcul presente dans ce chapitre, initialement congu comme un calcul de la probabilite de 
formation d'un DCC [^], correspond done plutot au calcul de la probabilite de formation 
d'un champ de pion classique. Ceci n'affecte cependant en rien la pertinence de notre 
resultat. 



Chapitre 3 



Description microscopique de la 
formation d'un DCC 



Nous avons vu dans les chapitres precedents que le scenario du trempage prevoit la 
formation d'un champ de pion classique lors du passage rapide de la transition de phase 
chirale dans une collision d'ions lourds ultra-relativistes. Ce scenario, initialement propose 



par Rajagopal et Wilczek en 1993 |16|, a ete largement admis comme une description 
microscopique de la formation d'un condensat chiral desoriente, un etat hypothetique du 
champ de pion propose au debut des annees 1990 ^, La correspondance entre 
la configuration du champ issue du trempage d'un etat thermique et celle proposee au 
debut des annees 1990 n'est cependant pas claire. Par exemple, les fortes fluctuations de 
la fraction de pions neutres prevues dans I'hypothese du DCC n'ont jamais ete observees 



dans le scenario du trempage |25, 3C 



La configuration DCC est caracterisee par trois points essentiels |14, 35, |6^ : 

1. c'est une excitation coherente du champ, une configuration essentiellement classique, 

2. chaque composante de Fourier du champ oscille dans une direction bien determinee 
de I'espace d'isospin, 

3. les orientations des differents modes sont alignees dans une seule et meme direction. 

L'objet de ce chapitre est I'analyse statistique detaillee de la configuration du champ de 
pion produit par trempage, le but etant de determiner si la structure decrite ci-dessus est 
realisee de fagon generique dans ce scenario. 



Pour ce faire, nous reprenons le modele original de ||16| dans lequel le systeme, ini- 
tialement dans un etat symetrique evolue selon les equations du mouvement du modele 
(T-lineaire classique. Les motivations de ce choix sont multiples. D'abord, les effets quan- 
tiques ne sont pas de premiere importance, le systeme etant essentiellement classique ||35|| . 
De plus, nous ne savons traiter ceux-ci de fagon simple que dans I'approximation de champ 
moyen oti les differents modes sont independants les uns des autres. Ce traitement n'est 
done pas adapte a I'etude des correlations entre les modes. Nous avons vu dans le chapitre 
precedent que dans un scenario plus realiste, ou le trempage est du a I'expansion rapide, 
le systeme n'entre que tres rarement dans la region d'instabilite, ce qui nous oblige a avoir 
une statistique importante pour I'etude des configurations interessantes. Dans le modele 



statique de |16| ou le trempage est realise de maniere artificielle, le systeme demarre dans 
la region spinodale et tous les evenements sont amplifies. Enfin, I'etude que nous nous 
proposons de realiser remet en question une idee largement admise. II est done preferable 
de raisonner dans le cadre le plus simple de fagon a ne pas obscusrcir notre argumentation. 
Nous commengons par rappeler en details le modele de Rajagopal et Wilczek, puis nous 
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construisons les outils nous permettant d'etudier la structure d'isospin du champ. Nous 
presentons ensuite les resultats de I'analyse statistique de I'etat final, que nous comparons 
a des travaux antecedents. Enfin nous discutons les implications de notre etude tant du 
point de vue plienomenologique qu'en ce qui concerne la comprehension microscopique de 
la formation d'un DCC. Ce travail a ete publie dans la revue Phys. Rev. D |7C]. 



3.1. Le formalisme 
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3.1 Le formalisme 

3.1.1 Le modele de Rajagopal et Wilczek 

Nous avons deja introduit ce modele dans le Chap. [l]. Nous le decrivons ici en detail. 
Le parametre d'ordre de la symetrie chirale est represente par un champ vectoriel a quatre 
composantes : (j) = (tTjCt), et son evolution temporelle est gouvernee par les equations du 
mouvement du modele ci-lineaire classique (cf. Eq. (|1.4|) ) 



- At;2 + A(/.2(x, t)) (t){x, t) = Hria , (3.1) 

ou les parametres v, X et H sont determines a partir des quantites physiques^ = 
135 MeV, = 92.5 MeV et m„ = 600 MeV a travers les relations 

H = Uml , 

Nous resolvons numeriquement ces equations sur un reseau cubique de maille a et de volume 
V = L^, ou L = Na. Pour ce faire, nous devons specifier les conditions aux bords, ainsi 
que I'ensemble des valeurs du champ cj> et de sa derivee temporelle dtcf) = cf) en chacun des 
noeuds du reseau a I'instant initial. 

Suivant |[l^ , la configuration initiale est echantillonnee a partir d'une distribution 
symetrique refletant un etat thermique a une temperature T > Tc : les valeurs de cj> 
et aux differents noeuds du reseau sont des variables aleatoires gaussiennes indepen- 
dantes de moyennes = = et de variances respectives^ {(pj) = f^/16 et ((^|) = v'^/A 
[j = 1, 2, 3, 4), en unites du pas du reseau. L'independance statistique des valeurs du champ 
et de sa derivee en differents noeuds du reseau signifie que le pas a a un sens physique : c'est 
la longueur de correlation du champ (degeneree en tt et cj) a la temperature T. Nous pren- 
drons, avec Rajagopal et Wilczek, a = (200 MeV)"^ = 1 fm. Nous avons deja mentionne, 
dans le Chap, [l], le fait que bien que les valeurs des variances du champ et de sa derivee 



soient choisies arbitrairement dans |16], ce choix a regu une justification avec les travaux 
ulterieurs de Randrup |29| : la valeur donnee ci-dessus est approximativement egale a la 
valeur typique du carre des fluctuations, apres que I'expansion ait emmene le systeme au 
coeur de la region d'instabilite spinodale (cf. Fig. |lj^). Ce trempage artificiel permet de 
nous limiter aux configurations qui subissent une forte amplification, les seules qui nous 
interessent ici. 

Nous denotons par ^{k,t) et <p{k,t) les composantes de Fourier du champ et de sa 
derivee temporelle a I'instant t et choisissons des conditions aux bords de Neumann, de 
sorte que ces composantes soient reelles. Ce n'est pas la un point essentiel, mais ce choix 
s'avere plus judicieux^que des conditions periodiques pour ce qui est de discuter le probleme 



^Nous suivons ici le choix de |0 pour les valeurs des quantites physiques, de fagon a tester nos resultats 
en les comparant aux leurs. 

^Dans EB], les variances sont donnees pour les longueurs des vecteurs (p et 4>, par exemple {<p^) — 

^Avec des conditions aux bords periodiques, les composantes de Fourier sont des nombres complexes 
dont les parties reelle et imaginaire jouent le role de deux degres de liberte distincts. C'est un artefact des 
conditions periodiques que ces deux degres de liberte soient reunis pour decrire un seul et meme mode 
de Fourier (a ce propos, voir 1' Annexe Pour discuter les questions des trajectoires des ondes k dans 
I'espace d'isospin, ainsi que des correlations entre celles-ci, il est preferable que les composantes de Fourier 
soient des nombres reels. 
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de la polarisation^ des ondes (^(/c, t) dans I'espace d'isospin. On a (cf. Annexe ^) 

(f{k,t) = (^J^ J d'^x cos{kxx) cos{kyy) cos{kzz) cf){x,t) , (3.2) 
ip{k,t) = \TJ J d'^x cos{kxx) cos{kyy) cos{kzz) 4>{x,t) ■ (3.3) 

L'espace de Fourier est un reseau discret de pas Ak = ir/Na et, les fonctions modes etant 
periodiques et paires dans chacune des directions, on se limite au cube < /cj < vr/a. 

L'evolution hors d'equilibre du systeme genere une amplification importante de I'ampli- 
tude des oscillations des ondes ip{k, t) correspondant au modes de grande longueur d'onde 



(cf. Fig. |3.lD , et ce dans chacune des directions d'isospin. Ce phenomene explique la for- 
mation possible d'un champ de pion fort (classique), c'est la premiere condition pour avoir 
une configuration DCC. Pour que les deux autres conditions enumerees dans I'introduction 
soient satisfaites, les amplifications de differentes composantes d'isospin d'un meme mode 
doivent etre telles que les ondes <^(A;, t) soient lineairement polarisees dans I'espace d'isospin 
(c'est a dire que le vecteur oscille le long d'un axe dans cet espace), et les amplifications 
de differents modes doivent etre correlees entre elles, de fagon a ce que les modes oscillent 
tous le long du meme axe. Nous reviendrons sur ces aspects de maniere plus precise dans 
la suite de notre expose, le point important ici est le fait que, dans une configuration DCC, 
ou plus exactement, dans un ensemble statistique dont une configuration generique est de 
type DCC, les differents modes sont fortement correles. 

Le point de depart de notre etude est la remarque suivante : les valeurs de (f) et (p en 
differents noeuds du reseau sont supposees etre des nombres aleatoires gaussiens indepen- 
dants. II est facile de voir (cf. Annexe |^) que cela implique que les valeurs des composantes 
de Fourier (f et cp sur les nceuds du reseau reciproque sont aussi des nombres gaussiens 
independants. En d'autres termes, il n'existe pas de correlations entre les modes dans I'etat 
initial. Pour creer une configuration DCC, le mecanisme du trempage doit, non seulement 
etre efhcace en ce qui concerne I'amplification des modes (et il Test), mais il doit aussi etre 
capable de generer des correlations entre les modes amplifies, ce qui est loin d'etre une 



evidence (cf. Eq. (|3.1j)) 



3.1.2 Observables 

Notre but est de tester si la configuration generique issue du scenario du trempage est 
de type DCC. En d'autres termes, nous voulons analyser certains details de I'ensemble 
statistique decrivant I'etat final, en particulier les polarisations des modes et leurs corre- 
lations. Pour ce faire, nous allons construire certains outils permettant de formaliser et de 
quantifier les notions dont nous avons besoin, en particulier la polarisation des ondes ip. 
Nous prenons le parti de construire le modele le plus simple possible, de fagon a ne pas 
obscurcir notre argument. 

Dans ce modele simple, bien que la taille de la boite soit fixee, nous avons en tete 
un systeme en expansion rapide (I'hypothese du trempage est une idealisation de I'effet 
de I'expansion). Dans cette image, apres un certain temps, le systeme est si dilue que les 
modes de Fourier decouplent les uns des autres et evoluent librement. La modelisation la 
plus simple^ de cet effet consiste a arreter "a la main" revolution couplee (regie par les 

^Dans la suite, nous appelerons (iso-)polarisation de I'onde k, la trajectoire dans I'espace d'isospin 
decrite au cours du temps par le vecteur ip{k,t). 

'Nous negligeons la duree finie de la periode de decouplage, pour laquelle une modelisation plus precise 
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equations du mouvement (3.1)) a un temps de decouplage tj (freeze-out). La configuration 
du champ a I'instant tf joue le role d'une source classique rayonnant des pions, on pent la 
voir comme la condition initiale de revolution ulterieure libre, les fonctions ^j{k,t > tf) 
representant les ondes associees a la propagation des pions rayonnes par cette source. Les 
proprietes de ceux-ci sont done entierement determinees par la configuration du champ a 
I'instant de decouplage tf. 

Donnons un aspect plus formel a cette image. Les quanta emis par une source classique 
sont dans un etat coherent (dans la suite, nous nous limitons au secteur des pions) 

\a,tf) = exp d^k a{k,tf) ■ a<{k)^ |0) , (3.4) 



avec 

3 

a • = Oj-aj , 
i=i 

ou aj{k) est I'operateur de creation d'un pion libre, de composante d'isospin j et d'impulsion 

k, et aj{k,tf) est la valeur propre de I'operateur d'anihilation correspondant {tf est un 
parametre), correspondant a I'etat propre (^^ ). C'est aussi la composante de Fourier de 
la source classique qui rayonne les pions |p5[ , elle est determinee par la configuration du 
champ classique a I'instant tf k travers la relation|^ 

a(fc;i,) = ^^^^lM±^I^SiM, (3.5) 



ou LOk = \/^ + rn^. Nous avons deja mentionne le fait que I'etat (|3.4| )- (|3.5D doit etre vu 
comme la condition initiale pour revolution ulterieure (t > tf), supposee libre et, de ce 
fait, contient toute I'information utile. En fait, I'idee de revolution ulterieure libre n'est 
pas necessaire, I'etat du systeme a I'instant tf etant tout ce qui nous interesse. II est 
cependant instructif d'introduire cette idee qui nous permettra de nous faire une image 
physique claire et simple des quantites pertinentes a mesurer. Pour alleger les notations, 
nous omettrons dorenavant d'indiquer explicitement la dependance en tf. Nous ecrirons 
par exemple, cx{k) pour a{k,tf), et a{k,t) pour cx{k,t > tf). De plus, tant que nous 
focaliserons notre attention sur un mode k donne, nous omettrons I'indice k, lequel sera 
reintroduit quand ce sera necessaire. L'evolution ulterieure s'ecrit done 

a{t) = ae-''^^^~^fK (3.6) 

Le nombre moyen de quanta d'isospin j associe a cette onde est independant du temps et 
vaut 

fij = {a\a^jaj\a) = \aj{t)\'^ = \aj\'^ . (3.7) 
Enfin, la fraction / de pions neutres dans un mode k donne est (0 < / < 1) 

/ = ^ • (3.8) 

ni+n2+ 71-3 



du phenomene est necessaire. Nous supposons que tous les modes decouplent au meme instant t / , qui peut 
etre vu comme la fin d'une periode commune de decouplage, laquelle est implicitement supposee etre tres 
courte devant I'Sclielle de temps caracteristique ~ 1 fm. 

^On peut voir la configuration classique du cha mp a I'instant tf comme la valeur moyenne du champ 



quantique correspondant dans I'etat coherent (3.4) 
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Pour caracteriser la structure d'isospin de la configuration du champ a I'instant tf, nous 
introduisons le concept de polarisation des ondes sortantes dans I'espace d'isospin : c'est la 
trajectoire decrite par I'extremite du vecteur 93°"* (t) = ip{t > tf) dans I'espace d'isospin. 
On a 



Rea(t) 

UJ 

2uj Ima(t) . 



(3.9) 
(3.10) 



Tout d'abord, il est facile de voir que ce mouvement est planaire : en effet, le vecteur 
I = <^°"* X (^°"*, qui est I'analogue du moment angulaire en mecanique du point]], est 
independant du temps. La trajectoire de 93°"* est done une ellipse^ dans le plan perpendic- 
ulaire a I. Appellons u et L (v et /) la direction et la longueur du grand (petit) demi-axe 
de cette ellipse, comme indique sur le schema de la Fig 3^1 (u^ = = 1, u.v = 0, 
I = ojILvl X v). On a alors 



cx{t) 



(3.11) 



ou 77 est un facteur de phase a determiner a partir de I'Eq. (|3.6|), ainsi d'ailleurs que les 
longueurs L et I, et les directions u et v. Avec cette parametrisation, il vient, pour le 

















L 



Fig. 3.1: Trajectoire elUptique du vecteur <p°"* dans le plan perpendiculaire a I = <^°"* x (p° 



nombre d'occupation moyen (3.7) 



et la fraction de pions neutres 



/ 



L2+/2 



(3.12) 



'^Ig est la composante k du generateur des rotations dans I'espace d'isospin J dPx (t>{x, t) x (pix, t) oc 
J dPkip{k,t) X ip{k,t), lequel est conserve. Pour t > tf les modes sont decouples et chaque composante Ig 
est conservee. 

^C'est la trajectoire periodique caracterisee par une echelle de temps unique (ici l/to) la plus generale 



3.1. Le formalisme 



49 



Une mesure directe de la polarisation de I'onde sortante (p'"** est donnee par I'excentricite^ 
de I'ellipse ( ^.11 ), definie comme le rapport des longueurs des petit et grand demi-axes 



e = //L (0 < e < 1). Le calcul de cette quantite a partir de la configuration a I'instant 
tf est simple. Le nombre total de pions dans le mode k, n = Yl^=i "-jj ^st egal a I'energie 
dans le mode k divisee par I'energie d'un pion : 



n 



2uj 



ou Ton a note (/)^ = • et de meme pour ip^ . En utilisant la parametrisation ( ^.11|) a 
t = tf, on obtient les relations 



f + L^ = ln, IL=- 

u> u> 



ou I est la longueur du vecteur d'isospin I = (p{tf) x ^(tj) defini plus haut. On en deduit 
les longueurs 



et I'excentricite 



2 n — ^/n^ — P 2 — P 
L — , Ij = 



2 n - Vn^ - P 
n + vn^ — 



Avec ces definitions en main, considerons le cas de la polarisation rectiligne (ou lineaire). 
L'onde k est polarisee lineairement si le vecteur (^|"* oscille le long d'un axe donne : yj^"* et 
0^"* sont collineaires. La polarisation rectiligne est done caracterisee par le fait que = 0, 
ou encore : = 0, = 0. On a done 

OLUnaireik^t) = a(fc) 6"*"'= , 

fUnaire{k) = COS^ 6*^ , (3.14) 

OU Oj: est Tangle entre la direction d'oscillation et I'axe vra de I'espace d'isospin. Aussi bien 
la dynamique (Eq. ^H] ) que I'ensemble statistique decrivant I'etat initial sont invariants sous 
le groupe 0(3) des rotations dans I'espace d'isospin : il n'existe aucune direction privilegiee. 
Done, si une telle onde, polarisee lineairement, est produite de fagon generique dans I'etat 
final, la distribution des valeurs de la fraction neutre f{k) sera donne^^par la loi en l/vT 
(cf. Eq. (|l.3|) ). Ceci reste approximativement vrai dans le cas, plus realiste, ou les ondes 
sortantes ont une polarisation quasi-rectiligne (/^ <C -L^). II est important de remarquer 
que, s'il est toujours possible d'ecrire la fraction neutre f{k) (Eq. (|3.8|)) comme le carre du 



®Pour se faire une image du sens physique de I'excentricite, il suffit de penser a une experience de 
polarisation d'ondes lumineuses (dispositif de type polariseur-analyseur). Dans le cas des "ondes de pions", 
on decompose la polarisation sur la base des etats de polarisation lineaire Ittq) et j7r+), les roles des 

analyseurs etant joues par des detecteurs sensibles a ces etats de charges respectifs. 

^"Dans un espace a trois dimensions, la direction aleatoire u est reperee par les deux angles {9,(fi). Si 
toutes les directions sont equiprobables, la probabilite pour que la direction u soit dans Tangle solide 
= sinOdOdif centre autour de la direction (dj'p) est (ffl/i-jr. La probabilite pour que coa^ 6 soit 
compris entre / et / + df s'ecrit P{f)df, avec 



P{f)= f sine de 5{f -cos^e)= [ dxSU-x^) 

Jo Jo 
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cosinus d'un certain angle x^) celui-ci n'a pas, en general, la signification d'une orientation 
aleatoire uniforme dans I'espace d'isospin. La loi en l/\/7 n'est valable que dans le cas de 
la polarisation lineaire (voir par exemple 1' Annexe |^. 

Considerons maintenant la configuration DCC "ideale", telle qu'elle a ete proposee au 
debut des annees 1990 : tons les sont lineairement polarises, et oscillent dans la mime 
directiorJP] u : 

aDCc{k,t) = a{k)e-''^^ u. 
Definissant la fraction du nombre total de pions neutres 



avec Nj = J d^knj{k), 



on a 



ou 6 est Tangle entre la direction u et I'axe tt^. Pour un DCC ideal, la fraction neutre totale 
est distribuee selon la loi en l/vT- Dans une situation plus realiste, on s'attend a ce que 
seuls les modes amplifies contribuent a la structure du DCC. II faudrait done ne consideder 
que les modes de grande longueur d'onde dans ( ^.15 ). De plus, il est peu probable que les 



ondes aient une polarisation strictement rectiligne, ni meme qu'elles soient strictement 
alignees les unes avec les autres, et on s'attend a des deviations a la loi ideale. Le point 
clef de la configuration DCC est le fort degre de correlation entre les modes. Introduisons 
la fonction de correlation 



\/{0!(k)),{02(t' 



ou {AB)c = {AB) — {A){B), (...) designe la moyenne statistique, et O represente une 
observable quelconque. Les observables pertinentes pour I'etude des correlations entre po- 
larisations, sont la fraction neutre / (Eq. ( |3.8| )) et I'excentricite e (Eq. ( p.l3D ). 

3.2 Les resultats 

Nous travaillons avec un reseau cubique de taille N = 64. Les equations du mouvement 



(3.1) sont resolues numeriquement a I'aide de I'algorithme dit "straggered leap-frog" 
avec un increment de temps At = 0.04a. Les composantes de Fourier ( |3.2P -(^) sont 
calculees numeriquement par une methode de type "Fast Fourier Transform", adaptee aux 
conditions aux bords utilisees ici (voir Annexe 0). Avec une configuration initiale 
donnee, on pent suivre revolution temporelle des (f{k,t). Nous reproduisons completement 
le resultat de la Ref. ||l6|, que nous rappellons brievement ici (voir aussi le Chap. ||). 
Nous calculous a chaque instant la moyenne de ip'j{k,t) sur le volume de I'espace des 

phases delimite par les spheres de rayons k ± 6k/2, ou k = \\k\\, avec 6k = 0.057a~^. 
L'evolution temporelle de cette quantite est tracee, pour differentes valeurs de k, sur la 



Fig. ^ dans le Chap. 1^. Dans les directions d'isospin {j = 1,2,3), les modes de basse 



^^Le DCC est un etat d'isospin total nul : / d^k = 0. Ceci est clairement specifie dans la ref. |35|, et est 
une hypothes e, p lus ou moins explicite, dans tous les articles originaux, ou I'idee du DCC a ete proposee. 

^^Les Figs. |l]l|et|L^ ont ete obtenues avec des conditions aux bords periodiques, pour des raisons tech- 
niques : le processus de moyennage de la Ref. decrit dans le texte, est alors plus facile h implementer. 
La quantite a moyenner est alors \ipj\^. De plus, la normalisation des coefficients de Fourier a ete choisie 
identique a celle de pour faciliter la comparaison. 
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frequence sont fortement amplifies et oscillent en phase avec une periode vr/wfc — ixlra-,^. 
Ni amplification, ni oscillations collectives ne sont observees dans la direction a {j = 4). 
Ce comportement moyen s'explique, de maniere qualitative, a I'aide des arguments de 
champ moyen ||T^ exposes dans le Chap. |^. En resume : pour les temps courts (t < 10a), 
la courbure du potentiel effectif, c'est a dire la masse effective au carre, est negative (cf. 
Fig. p^ ), et les modes de grande longueur d'onde, dont la frequence est imaginaire pure, 
sont amplifies. C'est I'instabilite spinodale, que nous avons eu maintes fois I'occasion de 
discuter. L' amplification ulterieure (lOo ^ t < 50a) est due aux oscillations regulieres de 
la valeur moyenne du champ autour de sa valeur asymptotique /jr. C'est la resonnance 



parametrique |21, 23]. Le phenomene d'amplification est evidemment transitoire et, 
pour des temps suffisament longs {t > 100a), I'energie initiale est equi-repartie entre tons 
les modes, le systeme est dans un etat d'equilibre stationnaire. 

Ici, notre but est I'analyse statistique detaillee de la structure d'isospin de I'etat fi- 
nal, oil le champ a ete fortement amplifie, la question etant de savoir si des correlations 
entre les polarisations des differents modes ont ete generees avec I'amplification. Au vu 
de la discussion precedente, il est interessant de discerner entre les differents mecanismes 
responsables de I'amplification. Dans la suite, nous presentons done des resultats pour 
deux valeurs du temps de decouplage '■ tj = 10a (avec une statistique de 21 x 10^ evene- 
ments), qui correspond a la fin de la periode d'instabilite spinodale, et tf = 56a (10.9 x 10^ 
evenements), temps auquel I'amplification moyenne (cf. Eq. (|3.17| ) plus bas) est maxi- 
male. Par souci de clarte, nous ne considererons que les modes]^ k = {k = nA/c,0, 0), 
ou Ak = n/Na ^ 10 MeV. De plus, le modele a-lineaire etant une theorie effective aux 
echelles < 100 MeV, nous ne considerons que la fenetre < n < 15. 

Definissons le facteur d'amplification du mode k a I'instant tf : 

oil 

3 

P{k,t) = LOk "^fijik^t) 

i=i 

est la densite moyenne d'energie dans le mode k a I'instant t. Les nombres d'occupation 
moyens nj{k, t) correspondant sont calcules a partir de la configuration du champ a I'instant 
t, au moyen des Eqs. (^) et (|3.7|). La dependance avec k de la moyenne statistique du 



facteur d'amplification {A{k,tf)) est presentee Fig. |3.2| . Le comportement observe est bien 
celui attendu : les amplitudes des modes de grande longueur d'onde sont amplifies par 
rapport aux autres. Pour tf = 10a on note un accord semi-quantitatif avec les arguments de 
champ moyen (voir le Chap. I, Eqs. ([r^)-(|r7D) qui prevoient une amplification monotone 
des modes tels que k < amplification d'autant plus importante que k est petit. La fenetre 
des modes amplifies est considerablement diminuee et I'amplification moyenne augmentee 
dans le cas tf = 56a. 

Le comportement qualitatif de I'amplification moyenne en tant que fonction de I'impul- 
sion k, en particulier la decroissance monotone, est done assez bien decrit par I'approxima- 
tion de champ moyen. II est interessant d'aller plus loin et de caracteriser plus precisement 
la distribution des valeurs du facteur d'amplification. Les histogrammes de cette distribu- 



tion pour le mode zero, le plus amplifie en moyenne, sont representes Fig. |3.3| , aux instants 



^^En pratique nous calculous les modes (n, 0, 0), (0, n, 0) et (0, 0, n) et prenons la moyenne des trois 
contributions, exploitant I'invariance du probleme sous les rotations spatiales et augmentant ainsi la statis- 
tique. Le raisonnement ne s'appliquant pas pour n = 0, la statistique pour ce mode est trois fois moindre 
que pour les autres. 
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Fig. 3.2: Le facteur d'amplification moyen {A{k,tf)) en fonction de rimpulsion k (en MeV) pour 
tf = 10a (a gauche) ettf = 56a (a droite). Les barres d'erreur representent I'incertitude statistique 
et les lignes en traits tires servent de guides pour I'oeil. L 'amplification moyenne est une fonction 
monotone decroissante de I'impulsion, comme le prevoit I'approximation de champ moyen. 



tf = 10a et tf = 56a. Dans les deux cas on observe de tres fortes fluctuations du facteur 
d'amplification autour de sa valeur moyenne. Supposons que, pour un evenement donne, 
I'amplification A{k) dans le mode k soit superieure a sa valeur moyenne {A{k)). L'am- 
plification A{k + 6k) dans un mode voisin est-elle aussi superieure a sa valeur moyenne ? 
Autrement dit, le facteur d'amplification est-il, evenement par evenement, une fonction 
decroissante de I'impulsion? La reponse est non, comme le montre la Fig. 3.4 qui represente 
la fonction de correlation C_4 definie par I'Eq. (|3l^ ) avec O = A.La decroissance monotone 
du facteur d'amplification avec I'impulsion n'est qu'une propriete moyenne. Cette absence 
de correlations entre les amplifications dans differents modes0 n'a pas de consequences tres 
importantes du point de vue phenomenologique. En effet, en pratique on mesure le nombre 
total (charges et neutres) de pions moyenne sur des bin dans I'espace des phases. Les am- 
plitudes totale des modes etant independantes, ce moyennage sur des bins est equivalent a 
la moyenne sur I'ensemble statistique pour des bins suffisament grands. La quantite perti- 
nente, du point de vue de I'amplitude totale, est done I'amplification moyenne. On retrouve 
le fait (Fig. p.2| ) que le scenario du trempage rempli tres efficacement le premier point du 
cahier des charges, decrit dans I'introduction de ce chapitre, concernant la formation d'un 
champ de pion fort (classique). 

Venons-en maintenant a notre etude proprement dite, c'est a dire a I'analyse de la 
structure d'isospin des configurations dans I'etat final et commengons par le second point 
de notre cahier des charges : I'etat de polarisation des modes individuels. Les distributions 
des valeurs de la fraction neutre f{k) sont representees pour differentes valeurs de k sur 
les Figs. 3^ et 3^, correspondant respectivement a = 10a et tf = 56a. Bien que toutes 
presentent de fortes fluctuations autour de la valeur moyenne (/) = 1/3, on remarque un 
changement de profil lorsqu'on passe des modes amplifies a ceux qui ne le sont pas. En 
fait pour tons les modes en dehors de la fenetre d'amplification (n < 15 pour tf = 10a, 
n < 4 pour tf = 56a), les distribution des valeurs de / presentent le meme profil lineaire. 
II est interessant de remarquer que cette loi lineaire est precisemment celle que Ton obtient 
dans le cas d'un ensemble thermique de pions (voir I'Eq. ( E.22|) ) de I'Annexe Cette 



observation conforte I'idee selon laquelle ces modes sont deja thermalises [16, |T^. En ce 



^*Lsl longueur de correlation est inferieure a la coupure infrarouge AA:. 
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Fig. 3.3: Histogrammes des valeurs de 1' amplification dans le mode fc = 0, dent I'amplification 
moyenne est la plus grande, pour tf — 10a (a gauche) et tf = 56a (a droite). Les fluctuations de 
I'ampliflcation autour de sa valeur moyenne sont tres grandes pour les modes les plus amplifles en 
moyenne. Ces histogrammes sont les analogues de ceux de la Fig. 2A du Chap. Dans le cas present 
cependant, on selectionne les evenements interessant (amplifies) en trempant artificiellement le 
systeme. 
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Fig. 3.4: La fonction de correlation reduite Cyi(0, fc) (Eq. ( ^.16 )) en fonction de k (en MeV) pour 
tf = 10a (a gauche) et tf = 56a (a droite). Les fonctions de correlation du type CA{ko,k) et 
C^(fco,p), ou fco et fc sont des impulsions dans la direction k^ tandis que p est une impulsion dans 
une direction differente, montrent des profils similaires : les amplifications des differents modes sont 
statistiquement independantes. 
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qui concerne les modes amplifies, qui sont ceux qui nous interessent, la distribution en / 
est sensiblement la meme que dans I'etat initial (cf. Eq. ([E.23|) de 1' Annexe g). 



P(f) 



n=0 



0.2 0.4 0.6 0.8 1 
f 



P(f) 



n=15 



0.2 0.4 0.6 0.8 1 
f 



Fig. 3.5: Histogrammes des valeurs des fractions neutres dans les modes k — nAk {Ak « 10 MeV) 
pour n — {a droite) et n = 15 (a gauche), a I'instant tf = 10a. Les histogrammes correspondant 
aux modes n < 14 sont tous compatibles avec la distribution correspondante dans I'etat initial 
(indiquee en traits tires pour n = 0), tandis que pour n > 15 on observe une distribution lineaire, a 
compare r avec la loi 2(1 — /) (indiquee en traits tires) obtenue dans le cas d'un ensemble thermique 
(cf. Eq. (|E.22D de 1' Annexe |). 
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Fig. 3.6: Meme chose que sur la Fig. 3.5 pour tf = 56a. Les histogrammes correspondant aux 
modes n < 3 sont compatibles avec la distribution initiale. Tous les modes n > 4 presentent une 
distribution lineaire. 

L'analyse des excentricites des polarisations elliptiques permet d'etudier plus precise- 
ment ce dernier point. La Fig. p.7| represente I'excentricite moyenne {e{k,tf)) en fonction 
de k a differents instants. Dans les deux cas tf = lOo et tf = 56a, les modes pour lesquels 
I'amplification moyenne est importante ont une excentricite moyenne tres proche, bien que 
legerement inferieure, de la valeur correspondante dans I'etat initial. Dans le cas t/ = 56a, 
on voit que I'excentricite moyenne des modes qui ne font pas, ou plus, partie de la fenetre des 
modes amplifies {k > 50 MeV) est independante de k et tres superieure a la valeur initiale 
correspondante. Cette observation s'avere vraie, non seulement en moyenne, mais evene- 
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ment par evenement, au niveau des distributions des excentricites des differents modes, 
representees aux instants initial et final, tf = 10a et tj = 56a, sur les Figs. IJ, 3^ et 3.10| 
respectivement. Pour les modes amplifies, les distributions finales sont legerement decallees 
vers les petites valeurs de e, a I'inverse, les modes tliermalises ont une polarisation plus 
circulaire (grandes valeurs de I'excentricite) que dans I'etat initial. Pour ces derniers, la 
distribution est independante de k. Bien que les modes de grande longueur d'onde aient 
une polarisation de plus en plus lineaire a mesure qu'ils sont amplifies, il s'agit d'un effet 
tres faible. On voit par exemple sur les Figs. |3]^, 3^ et S.lOj , que la proportion d'evene- 
ments pour lesquels I'excentricite du mode zero est inferieure a 0.1 est de : 13% dans I'etat 
initial, 16% a i/ = 10a, et 18% a = 56a. 
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Fig. 3.7: L'excentricite moyenne {e{k,tf)) en fonction de k (en MeV) dans I'etat initial (carres), 
ainsi que pour tf — 10a (losanges) et tf = 56a (triangles). Les lignes sont des guides pour I'oeil. 

Les modes qui, a un instant donne, sont dans la fenetre d'amplification ne voient pas leur 
polarisation notablement modifiee, tandis que ceux qui ne ressentent plus I'amplification 
et thermalisent, se retrouvent avec une polarisation "thermique", d'autant plus differente 
de leur polarisation initiale que leur frequence est grande. Le mecanisme responsable de 
I'amplification ne modifie done pas I'etat de polarisation des modes amplifies, la distribution 
de la fraction de pions neutres dans ces modes est essentiellement donnee par la distribution 



correspondante dans I'etat initial (cf. Figs. |3^ et |3^ 



II est cependant interessant de noter que, bien que les modes soient loin d'etre polarises 
strictement lineairement, cela n'entraine que de petits ecarts a la loi en l/\/7. En effet ici, 
de meme que dans le cas ideal, la fraction de pions neutres dans un mode k donne fluctue 
tres fortement d'un evenement a I'autre, si bien que du point de vue d'un experiment at eur 
la difference n'est pas tres importante. 

Le dernier point de notre etude concerne la question de I'alignement entre les direc- 
tions d'oscillations des differents modes dans I'espace d'isospin. La Fig. ^.11 represente la 



fonction de correlation C/(0,/c) (Eq. (3.16)) en fonction de k. Pour les deux valeurs de tf 



etudiees, les fractions neutres des differents modes sont completement independantes les 
unes des autres dans I'etat final. II en va de meme pour les excentricites. Bien qu'en moyenne 
les amplitudes ipj{k,t) des modes amplifies oscillent en phase dans chacune des directions 
d'isospin, les vecteurs (p{k,t) oscillent dans des directions completement independantes : 
dans un certain sens, les differents modes sont comme autant de DCC independants. La 
consequence phenomenologique immediate est I'attenuation, voire la disparition des flue- 



56 



Chapitre 3. Description microscopique de la formation d'un DCC 



n=0 




P(e) 



2 



1 - 



n=15 



0.2 0.4 0.6 0.8 1 



Fig. 3.8: Histogrammes des valeurs de I'excentricite (Eq. (3.13)) dans I'etat initial, pour les modes 
n = (a droite) et n = 15 (a gauche). 
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Fig. 3.9: Meme chose que sur la Fig. 3J, cette fois dans I'etat final a t/ = 10a. Les distributions 
correspondant aux modes tels que n < 11 sont deplacees vers les petites valeurs de e par rapport a 
I'etat initial. Ce decalage est le plus accentue pour n = 0, et s'estompe a mesure que n augmente 
pour s'inverser a partir de rt = 11 : au dela les distributions sont decalees vers les grandes valeurs 
de e. 
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Fig. 3.10: Meme chose que sur la Fig. 3^, pour tf — 56a. Le meme phenomene de decalage de 
la distribution est observe : decalage vers la gauche pour n < 2, et vers la droite au dela. Ce 
decalage vers les grande valeurs de e pour les modes de grande frequence se stabilise cependant : 
les distributions correspondant aux modes n > 7 sont toutes identiques (voir la distribution de 
n = 15 ci-dessus). 
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Fig. 3.11: La fonction de correlation reduite C/(0,fc) (Eq. ( ^.l(j )) en fonction de k (en MeV) 
pour tf — 10a (a gauche) et tf = 56a (a droite). De meme que pour I'amplification (Fig. |3.4| ), 
on a calcule les fonctions de correlation entre differents modes, pour la fraction neutre ainsi que 
pour I'excentricite. II en ressort que les polarisations des differents modes sont statistiquement 
independantes. 
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tuations de la fraction neutre du nombre total de pions dans un bin d'espace des phases. 
Ceci est illustre sur la Fig. 3.12| qui represente la distribution de la fraction neutre lorsque 
les contribution de seulement quelques modes equivalents (c'est a dire ayant la meme dis- 
tribution individuelle) sont pris en compte]^. 
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Fig. 3.12: Histogrammes des valeurs des fractions neutres provenant de la somme des contributions 
des modes n = 0, ...,4 (a droite) et n = 0, ...,9 (a gauche), a I'instant tf = 10a. Dans ce cas 
tous les modes consideres ont des distributions individuelles de / essentiellement identiques. On 
voit que le signal est rapidement detruit par I'interference entre les contributions statistiquement 
independantes. 



3.3 Discussion et conclusion 
3.3.1 Resume 

Partant d'un ensemble symetrique pour le champ chiral 4> = {''^i'^)-, nous resolvons 
numeriquement les equations du mouvement du modele cr-lineaire, qui pent etre vu comme 
une theorie effective decrivant la dynamique des excitations bosoniques de basse energie 
de QCD. Le trempage des fluctuations de I'etat initial est un mecanisme tres efficace pour 



generer des oscillations coherentes de grande longueur d'onde et de grande amplitude p6 |. 

L'analyse statistique detaillee de I'etat final montre que la structure d'isospin du champ 
de pion y est sensiblement la meme que dans I'etat initial. Bien que differentes de la loi en 
l/v7i propre a une polarisation rectiligne, les distributions de la fraction neutre dans les 
modes individuels sont tres larges, ce qui est le point important pour la phenomenologie. 
Cependant, les polarisations des differents modes sont completement independantes les unes 
des autres. La consequence phenomenologique importante est la suppression des grandes 
fluctuations de la fraction de pions neutres lorsque les contributions de plusieurs modes 
sont prises en compte, par exemple dans un bin de I'espace des impulsions. 

Le point clef est de realiser que I'hypothese d'un etat initial spatialement desordonne 
implique que les modes de Fourier du champ sont initialement independants les uns des 
autres. Notre resultat pent done etre enonce comme suit : les non-linearites de la dynamique 
ne generent pas de correlations entre les modes. Le chaos que Ton a postule a I'instant 
initial se retrouve dans I'etat final. Cela contredit I'idee largement repandue, selon laquelle 



'II s'arit du rapport r-^ r— , ou A'^, = n, (A:). 

A'^i + iV2 + ^ J w 
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la configuration generee dans le modele le plus simple utilise jusqu'a present, est un DCC, 
tel qu'il a ete propose a I'origine. Ce modele permet d'expliquer I'aspect classique du DCC, 
pas sa polarisation hypothetique. 



3.3.2 Comparaison avec des travaux antecedents 

De nombreuses etudes concernant les implications plienomenologiques du scenario du 



trempage ont ete realisees depuis Parmi les auteurs de ces travaux, certains |25, pO|, |68||, 
s'interessant a la distribution de la fraction de pions neutres dans I'etat final, ont obtenus 
des resultats analogues a celui de la Fig. ^.12 . 



Etudiant I'etat final du modele de Rajagopal et Wilczek en termes de domaines d'ori- 
entation donnee (par analogie avec les domaines d'aimantation formes lors du trempage 
d'un feromagnetique) , S. Gavin, A. Gocksch et R. D. Pisarski |25] arrivent a la conclusion 



suivante : le systeme est forme d'un grand nombre de petits domaines orientes aleatoire- 
ment et dont la taille est de I'ordre de la longueur de correlation, ce qui se traduit par 
une distribution binomiale de la fraction neutre du nombre total de pionsf^ Par la suite, 
K. Rajagopal a souligne le fait que le systeme en question, etant hors d'equilibre, ne pent 



etre caracterise par une echelle de longueur unique ||30[, de sorte que I'image de la Ref. ||2^ 
n'est pas appropriee. Dans ce chapitre, nous avons montre que les differents modes du 
champ se comportent comme des DCC independants : ils sont amplifies (ce qui justifie 
I'approximation classique), ont une orientation differente de celle du vide physique, mais 
leurs polarisations dans I'espace d'isospin sont statistiquement independantes. Bien sur, en 
prenant la transformee de Fourier inverse, on arrive a la meme conclusion dans I'espace des 
configurations : les orientations du champ sont independantes d'un site a I'autre. II n'est 
cependant pas correct de dire que Ton a produit une plethore de petits DCC (domaines), 
chacun de la taille du pas du reseau. En effet, les modes de grande longueur d'onde, qui ont 
ete fortement amplifies et qui sont, par consequent les degres de liberte interessants, occu- 
pent tout le volume de la collision ou le champ a ete trempe. De plus, la representation dans 
I'espace des impulsions donne une image plus proche de ce que les experiment at eurs voient. 
Quoi qu'il en soit, les "domaines" sont ephemeres, les detecteurs enregistrent finalement les 
impulsions des particules produites. 



Dans la Ref. |30|, Rajagopal etudie les implications phenomenologiques de son mod- 
ele, notamment quant a la distribution de la fraction de pions neutres de basse energie, 
en decoupant I'espace des phases comme dans une experience reelle. II calcule, pour un 
evenement donne, la distribution sur I'ensemble des bins, de la fraction neutre des pions 
dont I'impulsion ||fc|| < 300 MeV. II obtient une distribution relativement piquee autour 
de / = 1/3, ce qu'il interprete comme : "an admixture of a l/\/J distribution." II a en 
tete I'image d'un condensat chiral desoriente, forme par la superposition coherente des 
modes de grande longueur d'onde, plonge dans le bain thermique incoherent forme par les 
autres modes. Cette image dans la representation en impulsion, est plus satisfaisante. De 
plus, I'observation du fait que des modes de differente nature (amplifies ou thermalises) 
contribuent de fagon incoherente, supprimant ainsi le signal, est parfaitement en accord 
avec notre resultat. Cependant, nous avons mis en evidence le fait que les modes de basse 
frequence n'agissent pas de concert pour former un condensat coherent, c'est le fait d'avoir 
moyenne les contributions d'un grand nombre de modes dans chaque bin qui est responsable 
de la suppression des fluctuations de la fraction neutre. 



'^Les auteurs en question etudient differents regimes de couplage. La conclusion decrite ici est celle 
obtenue dans le regime de couplage fort, pertinent a la physique des collisions d'ions lourds qui nous 
interesse ici. 
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Finalement, J. Randrup |68| calcule la distribution de la fraction neutre de pions dont 
les impulsions sont inferieures a differentes coupures. II n'observe une loi large (en 
que dans le cas ou le seul mode zero contribue. L'accord entre les resultat du calcul clas- 
sique exact et de I'approximation de champ moyen le conduit a dire : "each k contributes 
pions having an independant orientation in isospin space." C'est precisemment ce que nous 
avons obtenu ici. Cependant Randrup travaille dans un modele avec expansion, c'est a dire 
qu'il modelise I'ensemble de la collision. Par consequent, la plupart des evenements qu'il 
considere n'ont pas ete significativement amplifies (cf. Chap. |2|). Dans ce sens, son analyse 
est incomplete et ne permet pas de conclure quant a la question que nous nous sommes 
pose ici. 

Certains indices concernant la structure d'isospin generee dans le scenario du trempage 
ont ete releves dans la litterature, en particulier I'absence de fluctuations significatives de 
la fraction de pions neutres. Cependant, I'analyse detaillee de la configuration du champ 
issu du trempage n'avait encore jamais ete realisee, en fait la question de la correspondance 
entre cette configuration et le DCC n'avait jamais ete soulevee. Le travail presente dans ce 
chapitre permet de clarifier la situation. 



3.3.3 Speculations 

Le scenario du trempage sous sa forme la plus simple ne permet pas de generer une 
configuration de type DCC. L'inclusion d'effets quantiques dans I'approximation de champ 
moyen, ou encore de I'expansion, n'alterent probablement pas cet etat de faits, les premiers 
par construction, et la derniere n'etant qu'une fagon sophistiquee de modeliser le trempage 
lui-meme ainsi que le decouplage entre les modes. Notre etude indique que c'est I'absence 
de correlations dans I'etat initial qui est a I'origine du probleme. 

En fait, si le phenomene du trempage des fluctuations initiales semble assez naturel dans 
le contexte des collisions d'ions lourds a haute energie, I'hypothese selon laquelle le systeme 
est completement thermalise a I'instant initial est loin d'etre justifiee a I'heure actuelle. 
En particulier, dans ces systemes de petite taille, les modes de grande longueur d'onde 
peuvent ne pas avoir eu suffisamment de temps pour thermaliser, avant que le phenomene 
d'instabilite spinodale ne se produise. On pent imaginer que des correlations, presentes 
dans I'etat initial, soient amplifiees lors de revolution hors d'equilibre ulterieure. En effet, 
on s'attend, d'apres les arguments de champ moyen, a ce que le phenomene d'amplification 
par instabilite spinodale se produise pour une grande classe de conditions initialesl^. Dans 
ce cas, le probleme se ramene a celui, hautement non-trivial, de la construction d'un etat 
initial realiste, ce qui necessite la description des premiers instants de la collison]^. 

A I'inverse, si des correlations etaient presentes dans I'etat initial, elles pourraient ne pas 
survivre a la dynamique non-lineaire. A nouveau, I'etat final serait constitue d'une super- 
position incoherente de modes desorientes. Dans cette situation, pour mesurer de grandes 
fluctuations!^ fraction de pions neutres dans un bin donne, il est necessaire d'isoler 
aussi proprement que possible les modes individuels. Le volume du systeme doit alors etre 
suffisamment grand pour que I'on ait une statistique satisfaisante, et suffisamment petit 
pour que deux modes voisins soient bien separes, la coupure infrarouge etant inversement 
proportionnelle a la taille du systeme. 



^'^La principale restriction etant que les fluctuations initiales ne soient pas trop importantes, de fagon a 
ce que le systeme puisse entrer dans la region d'instabilite. 
^*Le modele de la Ref. Jt^ ] est une possibilite interessante. 
^^Ces fluctuations sont le signal d'un etat semi-classique. 
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3.3.4 Conclusion 

Nous avons montre que I'idce selon laquelle le mccanisme du trempage permet d'expli- 
quer la formation d'un condensat chiral desoriente lors d'une collision d'ions lourds, idee 
tr6s largement r6pandue, est incorrecte. L'hypothetique structure d'isospin du DCC tel 
qu'il a ete proposee a I'origine, n'est pas generee par la dynamique. II n'existe a ce jour 
aucun modele microscopique capable de decrire la formation d'un DCC. C'est la le resultat 
essentiel de ce chapitre. 

D'un autre cote, il est possible que I'image originale du DCC soit par trop idealisee. 
Son extreme oppose est celle d'une superposition incoherente des modes du champ de pion 
classique. Nous avons vu que la manifestation la plus frappante de ces "ondes classiques 
d'isospin" reste la largeur inhabituelle de la distribution de leurs fractions neutres individu- 
elles. 

Pour aller plus loin, des developpements sont necessaires, notamment en ce qui concerne 
la description de I'etat initial, a moins que des donnees experiment ales ne viennent trancher 
la situation. 
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Annexe A 

Creation de particules dans une 
geometrie en expansion 



Dans un espace de geometrie (riemannienne) quelconque, certaines difRcultes appa- 
raissent, comme par exemple le probleme de la definition d'une energie conservee, qui 
empechent la generalisation de concepts usuels en theorie quantique des champs. Dans le 
probleme qui nous concerne au Chap. |2[ nous travaillons dans I'espace plat de Minkowski 
et le probleme de I'energie est simplement une consequence triviale du fait que nous avons 
choisi un systeme de coordonnees curviligne : I'energie dans un co-volume decroit du fait 
de I'expansion (provoquant le trempage des fluctuations). Bien entendu, I'energie totale du 
systeme est conservee. Cependant tout n'est pas aussi facile. En particulier, nous allons 
voir qu'il n'est point besoin de courbure de I'espace pour que le concept de particule de- 
vienne ambigii (ce qui est embetant car nous cherchons precisement a calculer le nombre 
de particules). Ce probleme a une interpretation physique tres simple : nous sommes dans 
un referentiel accelere. D'apres le principe d'equivalence d'Einstein, cela revient a dire que 
nous sommes dans un certain champ de gravitation, lequel est couple au champ de matiere 
(par exemple par I'intermediaire du terme de masse \/—g et agit a la maniere d'une 
source classique en creant des quanta de ce dernier. Definir le concept de particule revient 
a definir le concept de vide. Choisissons une definition du vide a un instant tq, ce "vide" 
est rempli de "particules" a un instant ulterieur r. Nous allons discuter ce probleme sur 
I'exemple d'un champ scalaire massif libre. L'action classique s'ecrit 



-9 (5M.9^'/'9> + mV) 



Le tenseur d'energie impulsion est |t73|, p8|| 
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et satisfait a I'equation 



0. 



En utilisant la propriete V^g^u = 0, on en deduit I'equation du mouvement pour le champ 
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On remarque que la trace T^* 7^ dans la limite m — > 0. Cette forme du tenseur energie- 
impulsion n'est pas invariante sous les dilatations d'echelle dans la limite oii la theorie Test. 
Pour remedier a ce probleme, on utilise le fait que le tenseur energie-impulsion est defini a 
un terme de divergence nulle pres. On remplace done T^u — > T^u + 5T^u, avecQ 

En utilisant la propriete V"^g^v = 0, on obtient 

Le terme entre parentheses est antisymetrique en a, z/, tandis que, dans le cas d'un es- 
pace plat ou les derivees covariantes commutent^, V^V°^ est symetrique. On a done bien 
V^5T^y = 0. En utilisant I'equation du mouvement, on obtient = m'^cfP', qui s'annule 
bien dans la limite conforme. 

Le Hamiltonien du systeme s'ecrit H{t) = f d?x VhTrr, ou \fh est le determinant 
de la metrique decrivant la pseudo-sphere de rayon r 

H{t) = t^ Jd^xVh Q(a,0)2_-L^A(3)</>+!^</>2^^(<^9.0 + a,0<A)) . (A.l) 



Quantification 

Le moment conjugue du champ (p est 

dOr(p 

et les relations de commutation canoniques s'ecrivent 

[0(r, x) ■ n(r, X ')] = i S^^^ {x - x') , (A.2) 

les autres commutateurs a temps propres egaux etant nuls. Dans le systeme de coordonnees 
de Minkowski, la theorie libre s'exprime comme un ensemble d'oscillateurs harmoniques 
independants. Dans notre systeme de coordonnees curvilignes (et, plus generalement dans 
une geometrie courbe), le champ "libre", c'est a dire sans auto-interaction, est couple a un 
champ classique externe : la metrique. Nous allons voir cependant que le probleme pent 
etre ramene a un ensemble d'oscillateurs independants dont les frequences dependent du 
temps. Introduisons tout d'abord la variable de temps adimensionnee u = In(mr) [u est 

^Ceci revient a remplacer + dans I'action, ou R est le scalaire de Ricci (voir On 

obtient alors directement le tenseur T^i/ + ST^i, en utilisant les formules suivantes : 

lip V(7 c 

-9 9 5gp„ 

Dans notre cas R^" = et f? = 0, les equations du mouvement ne sont done pas modifiees. On peut voir 
que, dans le cas d'une geometrie de courbure non-nuUe, le terme £,R(I)^ assure I'invariance d'une theorie de 
masse nulle sous les transformations conformes de I'espace, la valeur du coefficient ^ depend de la dimension 
de I'espace, ici d = 4 et ^ = 1/6 j58|. 

^Le tenseur de courbure R'^" est proportionnel au commutateur [V ; V"]. 
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aussi appelle le "temps conforme"), dans la suite toutes les fonctions exprimees comme 
fonction de u portent un tilde : /(r) = f{u). Definissons les quantites adimensionnees 

(p{u,x) = T(j){T,x) , 7r{u,x) = Tdr^ir^x) = du'fiiu^x) , H{u) = tH{t), 

ainsi que les projections sur les fonctions propres du laplacien tri-dimensionnel A^^) (cf. 



Eq. (IJ) du Chap. D 



'fsiu) = J cfixVhyg{x)i^{u,x) , 
ng{u) = [ d^xVhy*g{x)H 



[u,x 



(A.3) 
(A.4) 



qui sont telles que (ft = {—l)"^Lp_g et vrt = (— l)™7r_^. En termes de ces quantites, 
les modes s et —s apparaissent couples. Pour rendre manifeste le decouplage des modes, 
introduisons les variables normales ||7^, definies ci-dessous pour m > 



En utilisant les Eqs. (U), et (|Aj)), on obtient 

les autres commutateurs etant nuls. Le hamiltonien s'ecrit 

'1 ,0 wK^) 



H{u) 



ou ijJsW) = s + T m . On a done, comme annonce, une superposition d'oscillateurs in- 
dependants dont les frequences dependent du temps^. Les equations de Heisenberg s'ecrivent 

q'(u) = -i H{u) ; q{u) = p{u) , 

p'{u) = —i H{u) ; p{u) = —Cj{u) q{u) , 

ou Ton a note /'(n) = df /du. La frequence physique de nos oscillateurs est 



t \ / " , 2 

a;,(T) = — — = y-2+m2 



(A.5) 



L'effet du (pseudo) champ gravitationnel, c'est a dire de I'expansion, est de decaler les 
frequences vers le rouge au fur et a mesure que le temps avance. On voit apparaitre le 
probleme de la definition de la notion de particule. Nous allons maintenant preciser ce 
point. 



^Dans le Chap, g, la dependance t emp orelle des frequences a deux origines : la geometrie, qui donne 
lieuy au terme s/r dans I'expression (A.5) des frequences physiques, et I'auto-interaction du champ, qui 



donne lieu a la dependance temporelle de la masse effective des quasi-particules dans rapproximation de 
champ moyen. 
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Creation de particules, approximation adiabatique 



Introduisons les operateurs de creation et d'anihilation ot et a, 



qui satisfont les relations de commutation 



(s- s) 



(A.6) 
(A.7) 



(A.8) 



les autres commutateurs etant nuls. La fonctions complexe ips est solution de I'equation 

f" + Cosiu)f = 0. (A.9) 
Les relations ( [A.2|) et( A.8 ) sont equivalentes pourvu que^ 

Ws = rJiu) - rsiu) ^'siu) = i . (A. 10) 

Soit uq (ou To), I'instant de reference, auquel les representations de Heisenberg et de 
Schrodinger coincident, on a 

^g = V^^(no)a,-. + 4*'K)(-iral,-, 

et denotons par U{u,uq) I'operateur unitaire qui connecte ces deux representations^ 

ifs{u) = U{u,uo)ipgU^^{u,uo) 

= Muo)as{u) + rs{uo){-iy^~al,{u), (A.ll) 

7rg{u) = U{u,Uo)TTgU''^{u,Uo) 

= i^',{uo)dg{u) + r;{uo){-iraU^), (A.12) 



ou dg{u) = U{u,uq) dgU ^{u,uo). D'apres 



( [A7| ) et ( |A.11| )-( |AT2|) , et en utilisant 



dg{u) = jls{u)dg+ i>s{u) {-l)"" aig, 



avec 



ii>*{u) = ^ps{uo)i)'^{u) -7p'^{uo)'>ps{u) . 



(A.13) 



(A.14) 
(A.15) 



La transformation de Bogoliubov ( |A.13| ) est unitaire 



1. 



Supposons que le systeme soit initialement (a I'instant uq) dans I'etat |0, mq) defini par 
la relation a^|0,uo) = 0,Vs. Le nombre de particules inital dans le mode s est done 



^Le Wronskien Ws est conserve : W' = 0. 

®Cet operateur peut etre construit explicitement [fea 
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(0, Mol o^«s |0,uo) = 0. Le nombre de particules dans le mode s a I'mstant u est facilement 
obtenu a I'aide de ( |A.13D , ( |A.14D et ( |A.15| ) : 



risiu) = {0,uo\ di{u) dg{u) \0,uo) = \i>s{u)y 



On voit explicitement le probleme mentionne au debut de cette partie : differents choix de 
la fonction ips correspondent a differents operateurs ag et done differents vides |0,no). De 
fagon generale, aucun choix de ag ne correspond a I'operateur d'anihilation de particules 
physiques, le nombre de particules n'etant pas constant dans une metrique dependant du 
temps. Le probleme de la dependance en temps du vide est assez serieux : par exemple, il est 
habituel, en theorie quantique des champs, de soustraire au Hamiltonien sa valeur moyenne 
dans le vide, infinie, les seules quantites physiquement observables etant les differences 
d'energie entre etats. Si la definition du vide depend du temps, cette procedure n'a plus 
de sens car on doit alors soustraire des quantites differentes a des instants differents, on ne 
peux pas comparer deux mesures faites a des instants differents. 

Notons de plus que la dependance temporelle du nombre moyen de particules rend 
la mesure de cette quantite tout a fait incertaine, autrement dit, le nombre moyen de 



particules n'a pas une signification physique bien definie ||63|. En effet, supposons que 
notre observateur en co-mouvement cherche a faire une mesure de cette quantite, mesure 
necessitant une duree At. Le resultat de la mesure est entache de I'incertitude 

AiV > (mAr)'^ + \A\At , 

ou le premier terme vient du principe d'incertitude de Heisenberg et ou A est le taux moyen 
de production de particule par unite de temps dans le co- volume considere. L 'incertitude 
minimale possible est ANmin = 2(|A|/m)^/^, avec At = {m\A\) — 1/2. II n'existe pas de 
valeur de At pour laquelle ANmin s'annule. 



L. Parker a etudie ce probleme vers la fin des annees 1960 [33| et a propose de definir les 
particules physiques, ou, de maniere equivalente, les fonctions modes ipsiu), comme etant 
celles qui minimisent ANmim c'est a dire qui minimisent le taux de production A. Cela 
revient a faire un developpement dans le taux d'expansion, c'est a dire dans un parametre 
decrivant la "vitesse" a laquelle la metrique change avec le temps. Intuitivement, il est clair 
que plus I'expansion est lente, plus on est proche du cas statique (Minkowski), et plus le taux 
de production de particules est faible : c'est le developpement adiabatique. Considerons, 
avec L. Parker et S. A. Fulling |Q, les solutions "de frequence positive" generalisees des 
Eqs. (|AJ|) et ( |A.10D : 



gs{u) = I exp ( -i 



I dvQ.s{v) \ , (A.16) 



ou Os(ii) est une fonction reelle satisfaisant I'equation (cf. (|A.9D ) 

L 'approximation adiabatique d'ordre est obtenue en identifiant, a un instant de reference 

V;s(no)=#(no) , ^'Auo) = -gi^Huo) , 
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ou gi^^ est donnee par I'Eq. ( A.16| ) en remplagant fj^ ^ 



(0) 



Wo 



#(no) 



2a;s(no) 



(p) 

Les approximations successives, gl , sont obtenues par recurence 



2 .2, ^ l^^i''^" 3/1^. 
n) I = a;^(u) - --777^ + 



ri(p)" 



2 Ol^) 4 \ 



Mentionnons enfin le fait que cette procedure de quantification permet de renormaliser la 



theorie de fagon non-ambigiie, c'est a dire de maniere independente du temps |34, 60, 58|. 
On voit en effet que I'approximation adiabatique est d'autant meilleure que s est grand : 
dans la limite s —>■ +oo on a lis ~ s, qui ne depend pas du temps. On comprend done qu'il 
soit possible de traiter les divergences ultraviolettes les plus severes une fois pour toutes 

a I'instant de reference uq. En fait on pent montrer |5|] qu'en identifiant les fonctions 

(p) - 
modes i/jg avec gl a I'instant uo, le taux de creation de particules dans le mode s a un 

instant ulterieur u decroit comme l/s^"^^. 



Annexe B 



Equilibre thermique et equilibre 
thermique local 

Dans cette Annexe, nous presentons le modele cj-lineaire dans I'approximation grand 
N, dans le cas ou le systeme est a I'equilibre thermodynamique a la temperature T. Nous 
calculous ensuite les variances des distributions gaussiennes utilisees pour echantillonner les 
conditions initiales au Chap. Nous utilisons des notations conventionnelles : x = {x, y, z) 
designe le vecteur position dans le systeme de coordonnees cartliesiennes et d^x = dxdydz 
I'element de volume infinitesimal. 



B.l Equilibre thermodynamique 

Considerons un volume infini ou le champ chiral est a I'equilibre thermodynamique, la 
matrice densite du systeme s'ecrit 

P=|e-^/^ (B.l) 

ou H est le Hamiltonien ei Z = Tr[p\ est la fonction de partition. Dans cet etat, le 
parametre d'ordre est uniforme et constant 

= (t>"r, (B.2) 

= 0, (B.3) 

avec la notation {O)^, = Tr[pO]. On definit le champ de fluctuation autour de cette valeur 
moyenne]^ 

5<l)a{t,x) = $a(t,f)-</.^, (B.4) 

5(t)a{t,x) = i^a{t,x). (B.5) 

Dans I'approximation de champ moyen grand A^, le systeme de particules en interaction 
est remplace par un ensemble de "quasi-particules" independantes de masse effective ^/x^. 
En d'autres termes, le hamiltonien est diagonnal dans la base "nombre" associee a ces 
quasi-particules. En denotant respectivement par a} g et ^ les operateurs de creation et 

d'anihilation d'une quasi-particule "thermique" d'impulsion k et de composante chirale a, 



^La fluctuation 8<p deflnie ici, est differente de celle deflnie au Chap. 
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on a, dans la representation de Heisenberg, 

d^k 1 



ou EkiXr) = + Xj, et kx 
relations de commutation sont 



(2vr)3 
d^k 



a t e"*^"" + a"^ ^ e' 

a,k 



kx 



(27r)2 



EkiXi 



a,k 



-ikx 



a,k 



(B.6) 
(B.7) 



EkiXr)^ — k ■ X. Avec notre choix de normalisation, les 



".fc ' b,k' 



{2TTf5a,5^^\k-k'). 

Le hamiltonien decrivant les fluctuations autour du parametre d'ordre s'ecrit 

d^k 



On obtient aisement 



{a\.a^^^,)^ = {2T,f5ab5^''\k-k')Nk{xT) , {%,k%S'^ 



b,k' IT 



OU N^^Xt) 1^ distribution de Bose-Einstein pour les quasi-particules : 



NkiXj 



1 



(B.8) 



(B.9) 



(B.IO) 



(B.ll) 



La masse effective ^/x^ des quasi-particules satisfait la relation d'auto-coherence (equation 
de gap, cf. (|2.13 )) 



A 



A h2 



k^dk 2Nk{xT) + 1 



27r2 



(B.12) 



ou 



4>^ ■ 4>T- La valeur du parametre d'ordre 4 
minimum (cf. Eq. ( 2. 11] )) et depend elle-meme de x 



2Ek{Xr) 

>j, est telle que le potentiel effectif est 

r 

(B.13) 



Les divergences apparaissant dans I'Eq. ( [B.12 ) sont absorbees dans les parametres nus 
u et A de la fagon exposee au Chap. Dans la limite chirale = 0, on voit, d'apres 
I'eq. ([B.13| ), que dans la phase brise^ (T < Tc), ou (j)^ ^ 0, la masse effective des pions est 
nulle, tandis que dans la phase symetrique oii (/>j, = 0, celle-ci est a priori non-nulle. Dans 
la cas ou H ^ 0, cette transition de phase disparait, mais la discussion precedente reste 
approximativement valable : on parle d'une "phase" de haute temperature, ou (p^ f-w, et 
d'une "phase" de basse temperature, ou (f)^ ~ /,r- La variation de Xt avec la temperature 
est montree sur la Fig. 60 pour differentes valeurs de la coupure A. Comme attendu, a 
mesure que la valeur de la temperature se rapproche de celle de la coupure, la presence de 
celle-ci se fait de plus en plus ressentii^. L'etude presentee au Chap. || concerne la valeur 
T = 200 MeV, ou la coupure n'a que peu d'influence, le point de comparaison T = 400 
MeV est a la limite de validite du modele (voir la discussion au Chap. |2|). 

^Dans raproximation grand A'', la transition de phase chirale est du second ordre et la temperature 
critique est donnee par |34 



_ 55|| Tc = vTW^/t 160 MeV pour N ■ 
^Notre modele n'a de sens physique que pour T ^ A. 
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Fig. B.l: La masse effective au carre Xt fonction de T, pour differentes valeurs de la coupure 
ultra-violette A. 

B.2 Equilibre thermique local 

Considerons un tout autre systeme : un petit volume Vq, splierique, ou le champ est 
a I'equilibre thermique local a la temperature T. Cette notion signifie que les fluctuations 
statistiques dans le volume Vq sont les memes que si celui-ci etait une sous-partie du systeme 
considere dans la partie precedente (a la meme temperature). II est bien evident que ceci 
n'a de sens que si ce qui se passe a I'interieur du volume Vq ne depend pas, ou peu, de ce 
qui se passe a I'exterieur, ce qui veut dire que la taille de notre systeme doit etre de I'ordre 
de grandeur de la longueur de correlation au moins. Dans ce qui suit, nous nous interessons 
a la distribution de probability des valeurs possibles des moyennes spatiales du champ et 
de sa dcrivce sur le volume Vq. Ces quantites s'interpretent, pour un observateur vivant 
a I'interieur de Vq, comme les valeurs du parametre d'ordre et de sa derivee temporelle. 
Utilisant la definition de I'equilibre local, nous delimitons une petite boule Vq dans le 
systeme en equilibre thermodynamique decrit plus haut, et calculous la distribution des 
fluctuations statistiques dans ce sous- volume. Pour cviter d'ecrir inutilement la variable de 
temps, nous travaillons dans la representation de Scrodinger. 

On cherche done la distribution de probabilite des resultats possibles et de la 
mesure des observables associees aux operateurs 




(B.14) 



(B.15) 




(B.16) 



Le Hamiltonien etant quadratique, la distribution cherchce, qui est en fait I'exact analogue 
de la distribution de Wigner, bien connue en mecanique quantique, est gaussienne. Les 
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moments du logarithme de Z(p, q), ou cumulants, d'ordre superieur a 2 sont nuls. Les 
valeurs moyennes, ou cumulants d'ordre 1, sont aisement obtenues a I'aide des Eq. (|B.2| ), 
[1) et (iH) : 



E 



0. 



Les variances et covariances (cumulants d'ordre 2) s'ecrivent (cf. ( |B.4D et (|B.5| )) 



Var [(pa 
Var 



^0 JVo 



d^xcPy{5<l)a{x)5<l)a{y)), 



V? 



Gov 



JVq 

1 



(B.17) 
(B.18) 

(B.19) 
(B.20) 
fB.21) 



les termes non-diagonnaux dans les indices chiraux sont nuls. Le hamiltonien etant quadra- 
tique dans le champ Scf) et sa derivee Scj) et ne contenant pas de termes croises du type > 



il est clair que la fonction de correlation apparaissant dans I'Eq. ( B.21 ) est nulle. Calcu- 
lous la variance de la moyenne spatiale du champ, Eq. ( [B.19| ). La fonction de correlation 
apparaissant dans cette formule est 



{6<t)a{x)5(t>a{y)\ 



d^k 2iVfc(x^) + l 



(B.22) 



Introduisons les coordonnees spheriques x = {x,6x,ipx) et y = {y,6y,ipy), oil les angles 
azimutaux 9x et By sont reperes par rapport a la direction du vecteur k. Les integrations 
sur les orientations possibles des vecteurs k, x et y sont triviales. Apres les changements 
d'echelle x — > x/R et y ^ y/R, il vient {k = \k\) 

'•+~ dk 



coth 



2T 



avec 



dx 



dyxy sm{qx) sm{qy) 



r{kRo) , 
(sing — q cosg)^ 



(B.23) 



(B.24) 



Ci-dessus on a denote a\ la variance cherchee (Eq. (|B.19| )). On obtient de la meme fagon 
la variance ( B.20| ), que Ton denote par a"^ ■ 



Vial = ARt 



dkEk coth( ^ ) J^{kRo) 



Reecrivons les formules ( |B.23| ) et ( [B.25| ) sous la forme generique {J^{q) = J^{—q)) 

"+00 



V^a^ = 2R:^r> 



dyg 



(B.25) 



(B.26) 



ou g{k) est une fonction de Ek = y/k'^ + Xt- Dans la limite ou Rq^/x^ ^ 1, on obtient le 
premier terme du developpement asymptotique de I'expression ci-dessus : 

f + OO 



dyT{y) . 



(B.27) 
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En utilisant 



J — ( 



vr 



on obtient (cf. (|2^ et (|]52|)) 



(7^ ~ 



(B.28) 



(B.29) 



Cette formule est valable quand le rayon i?o de la bulle est sufRsament grand devant la 
longueur de correlation At = la dependance en I/Vq traduisant le fait que les 

fluctuations a I'interieur du volume Vq ne dependent pas de son environnement : diverses 
cellules de taille i2o ^ sont statistiquement independantes. 

Dans le Chap. ||, nous sommes interesses a des valeurs de aussi petites qu'il est 
permis physiquement, c'est a dire -Rq ^ -^T- A strictement parler, notre modele est defini 



en presence d'une coupure ultraviolette A, et les integrales ( B.23 ) et ( B.25| ) doivent etre 
calculees avec cette coupure (notons que la coupure intervient aussi dans le calcul de 
Xt)- Cela n'influe que tres pen sur la valeur de la dispersion a\ (Eq. ( |B.23|) ), tant que 
Ai2o S> 1. En effet, la moyenne spatiale ( |B.14| ) supprime les contributions des modes de 
petite longueur d'onde, seuls les modes tels que A;i?o ^ 1 survivent. La Fig. montre 
la variation de a\ en fonction du rapport Rq/Xt = Ro^/X^ pour T = 200 MeV et pour 
diverses valeurs de la coupure. La presence de celle-ci se fait plus ressentir dans le calcul 



de la dispersion a2 qui, a cause du facteur Ek au numerateur dans I'integrand de ( B.25 ), 
diverge logarithmiquement dans 1 'ultra-violet. Cependant, pour des valeurs physiques de la 
coupure (A ~ 1 GeV), cette dependance logarithmique n'a que peu d'influence sur la valeur 



de (72. La Fig. [B.3| est I'analogue de la Fig. [B.2| , pour la dispersion o"2 des fluctuations de la 
derivee temporelle du parametre d'ordre. Dans les deux cas, les courbes en pointilles sont 
obtenues a partir de la formule asymptotique ( [B.2g| ). Les valeurs des dispersions ai et 02 
utilisees dans le calcul du Chap. |^ sont resumees dans le tableau B.l . 



T 




Ro/\t 


0-1 


0-2 


(MeV) 


(fm-2) 




(fm-i) 


(fm-2) 


200 


0.73 


1 


0.21 


0.40 


200 


0.73 


2 


0.11 


0.14 


400 


2.2 


1 


0.37 


1.07 


400 


2.2 


2 


0.19 


0.42 



Tab. B.l: Les differentes valeurs des dispersions utilisees dans le calcul presente au Chap. § 
(A = 800 MeV). 
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Fig. B.2: La dispersion cri (en fm ^) en fonction du rapport Rq/Xt pour T = 200 MeV, et pour 
differentes valeurs de la coupure A. La courbe en pointille represente la formule asymptotique 




Annexe C 

Calcul d'integrales 



Dans cet annexe, nous calculons les integrales utilisees pour renormaliser I'equation de 
gap au Chap. | (cf. Ref. H). 
Definissons les fonctions 



II est facile de voir que 

1 5/„(A,/i) 



(fc2+^2)n/2 



-n/„+2(A,/i). (C.2) 



Commengons par le calcul de Ii, necessaire dans la renormalisation de la masse. En ecrivant 
et en integrant par parties, on obtient 



, /-A 

/i(A,/i) = Av'A2 + /x2- / dk + /i2 

Jo 



En ecrivant 

k^ + ^2 



on obtient aisement 



2Ii (A, /i) = A^ \/l + x2 - ^? sinh"^ (^^^ , 



ou X = /x/A. En ecrivant 

ln(sinh a + cosh a) = In e'* = a = sinh~^(sinh a) , 

d'ou Ton tire 

sinh^^ a = ln(a + yl + cfl) , 
on reecrit I'expression precedente sous la forme 

. . X A2 / u2 / 1 + Vl + 2;2\ ^ ^ 

Ii(A,/x) = — yln . (C.3) 
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Dans la limite ou x <C 1, on a 

/i(A,^)c.i(^A2-^2ij^A^^... , (C.4) 

ou ... designe la partie non- divergent e dans la limite A — > +00. Les resultats precedents 
ont ete obtenus en supposant jj.'^ > Dans le Chap. § on a ausi besoin, pour isoler les 
divergences a chaque instant, du comportement asymptotique de Ii dans la limite a; <C 1 
pour //^ = X < 0. Dans ce cas, on ecrit jj. = i\iJ,\ et, en separant (A > ^) 

+ / ' 



il est facile de verifier que la formule ( |C.3 ) est toujours valable avec la prescritpion Inz = 
ln\z\ +z7r/2. Le comportement ( p.4|) est toujours valable en remplagant /i — > |^|. 

Pour la renormalisation de la constante de couplage, nous avons besoin de I'integrale 
Is. En utilisant (|C.2|) on obtient immediatement 



Dans la limite ou x ^ 1 



l3(A,^) = -== + ln ^ . (C.5) 



/3(A,^) ~ln- + ... . (C.6) 



Annexe D 

Scenario du trempage : I'etat initial 



Dans cette annexe, nous detaillons la discretisation de la transformee de Fourier avec 
les conditions au bords de Neumann, puis nous calculous la distribution de probability des 
valeurs des composantes de Fourier discretes dans I'etat initial du modele de Rajagopal et 
Wilczek (cf. Chap. |). 



Conditions de Neumann 

Considerons, pour simplifier les notations, le cas a une dimension. Nous donnerons plus 
loin la generalisation (triviale) a trois dimensions. Soit H{x) une fonction de la variable x 
sur le segment [0, L], avec des conditions aux bords de Neumann : H' {x = 0) = H'{x = L) = 
0, ou le prime designe la derivation par rapport a x. II est clair que, de la decomposition 
de cette fonction en onde plane exp(iA;x), ne susbsistent que les cos(A:x) qui, contrairement 
aux sin(A:x), sont compatibles avec ces conditions aux bords. 

Pour implementer ces conditions aux bords et la decomposition de Fourier correspon- 
dante, definissons la fonction G'(x), vivant sur I'axe reel — oo < x < +oo, periodique de 
periode 2L, paire, et telle que G{x) = H{x) sur le segment [0, L]. Les composantes de 
Fourier de G{x) sont 



Qn = / dx G{x) exp (2mn-—\ = \ — f dx H{x) cos{kx) 

\/2L J-L ^ 2L/ \l L Jq 



h{k), 



ou k = rnr/L. Les fonctions reelles h{k) sont les composantes de Fourier de H{x), elles 
sont paires {h{—k) = h{k)). 

Discretisation 

La discretisation du probleme ci-dessus pent etre faite de plusieur fagons, la plus simple 
consistant a echantillonner les valeurs de nos fonctions G et H k intervalles reguliers x = ja 
ou a est le pas du reseau de taille L = Na (N est le nombre de sites). Cependant, nous 
aimerions repeter la construction precedente pour les composantes de Fourier de H. Pour 
cela, il est judicieux d'echantillonner les valeurs des fonctions aux points x = a/2 + ja : 
Hj = H{x = {j + l/2)a). L'intervalle —L < x < +L correspond alors a — < i < — 1, 
la periodicite de la fonction H s'ecrit G^+j = Gj^n, tandis que la propriete de parite 
prend la forme Gj = On a, pour les composantes de Fourier, 



— > Go- exp n(? + -) =\ — > Hi cos [ — —] 



N-1 r^- n nr- N-1 

I Trn I 1 -I- I / "/ I \ 
\ N I 

j=0 



78 



Annexe D. Scenario du trempage : I'etat initial 



Les composantes /i„ satisfont a la relation "d'anti-periodicite" hn-N = —hn+N, ce qui 
limite I'espace de Fourier (ou espace reciproque) utile a la premiere zone de BrillouinQ : 
— (A^ — l)<n<A^ — 1. La relation de parite = /i„ ramene le nombre de degres de 
liberte independants aA^:0<n<A^— 1, comme il se doit (dans I'espace direct, on a 
N nombres reels (pj independants). La derniere egalite ci-dessus pent etre reecrite sous la 
forme matricielle (dans la suite nous prendrons a = 1) 

N~l 

hn=^WnjHj ^h = W-H, (D.l) 

i=o 



avec 



2 /7m(i + l/2) 
iv^°n N 



^nj = \l- COS ( ' ' ) . (D.2) 



Certaines proprietes de la transformation ( p.lj ) nous serons utiles pour le calcul de la 
distribution des valeurs /„ dans un modele gaussien. 

Transformation inverse 

On definit la matrice W"^ , transposee de la matrice W par la relation habituelle 

yy^- = Wni ■ 

Calculous les elements de la matrice WW^. Les elements non-diagonnaux sont nuls : 
2 ^"^ 1 

OU 



i=o 



(0<p = nibm< 2N =^ sin(7rp/2A^) 7^ 0). Dans le calcul des elements diagonnaux, il 
faut distinguer les cas n = m = Oetn = m/0. Ceci est relie au fait que la relation de 
symetrie /„ = /_„ est triviale pour n = 0. Au final on obtient 

{WW^)nm = 5nm{l+6n0). (D.3) 

En utilisant les identites 2cos(a) cos(6) = cos(a + 6) + cos(a — 6), et 

E»(^) = 5(i-(-in. 

n=0 

on obtient 

(WtW),fc = 5,, + l. (D.4) 



^La premiere zone de Brillouin s'etend de — A'^ a A^, autrement dit : — 7r/a < k < 7v/a. Dans le cas 
present, on a hjv = 0. 
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A partir de cette relation on deduit 

N-l 



N-l 



N-l 



n=0 



k=0 



j=0 



En ecrivant la somme dans le membre de gauche sous la forme (cf. Eqs. (|D.l|) - ([D.2|) ) 



N-l 



j=0 



on obtient la formule d'inversion |64| 



\ n=0 / 



(D.5) 



ou le prime signifie que le terme n = ne contribue a la somme que par la moitie de sa 
valeur. En utilisant I'une des relations precedentes, il est facile de voir que 



N-l 



■ N-l 



j=0 



(D.6) 



n=0 



Enfin, on a (Det[>V^] = Det[>V]) 

Det[WW^] = (Det[W])2 = 2. 



(D.7) 



Distribution gaussienne 

Dans le modele du Chap. ^, I'etat initial est caracterise par un ensemble statistique 
gaussien : les valeurs de chacune des composantes chirales du champ en chaque noeud 
du reseau direct sont des nombres aleatoires gaussiens independants. II en est de meme 
pour les valeurs des composantes de la derivee temporelle du champ. Imaginons que la 
fonction H{x) represente une des composantes chirales du champ ou de sa derivee. Les Hj, 
j = 0,...,N — 1 sont done des nombres gaussiens independants de variance {H'^). Qu'en 
est-il des composantes de Fourier /i„ ? 

La relation ( p.3| ) implique que les composantes hn sont des variables aleatoires gaussi- 
ennes independantes. En effet, on a {{HiHj) = 6ij {H"^)) 

{hn hm) = {H^) (W W%m OC 6nm , (D.8) 

ainsi que 

{hn km hp kg) = YWni^^mj^^pk^ql{HiHjHkHi) 

ijkl 

= {H^f [(W {W + {W (W W^)mq + {w (W W^), 

= {hn hm) {hp hq) + (/l„ hp) {hm hg) + {hn hg) {hm hp) , 

et ainsi de suite. 
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Annexe D. Scenario du trempage : I'etat initial 



De fagon plus precise, et plus concise, on ecrit la probabilite 

N-l N-l 



Ph{H) W dHj = Ph{h) n dK, 

j=0 n=0 



d'oii Ton tire la densite de probabilite (cf. Eq. ( D.7 )) 

Pf^{h) = ^PH{W-^ -h). (D.9) 

Le facteur l/\/2 est du au fait que le mode zero a une variance egale au double de celle des 
autres modes. Comme on I'a remarque plus haut, cela vient de ce que la relation de parite 
hn = h-n est trivialement satisfaite pour le mode zero : ce degre de liberte est en quelque 
sorte moins contraint et fluctue plus. Dans le cas qui nous interesse, on a 

1 / H-H\ 

et done 

n=0 

ayecal = {H^)il + 5no)■ 



PUh) = n exp , 



Le cas tri-dimensionnel 

La generalisation des considerations precedentes a trois dimensions est immediate. La 
fonction H{x), sur le cube < x,y, z < L = Na de volume V = L^, est echantillonnee aux 
points X = a(l/2 + i, 1/2 + j, 1/2 + /c), ou sa valeur est notee = Hi j^k. On implemente 
les conditions aux bords de Neumann en introduisant une fonction G{x) ayant les meme 
proprietes de periodicite et parite que dans le cas uni-dimensionnel, par rapport a chacune 
des trois variables x, y, z. Considerant le cas d'une des composantes chirales du champ, les 
composantes de Fourier s'ecrivent (a = 1) 

N-1N~IN~1 
i=0 j=0 k=0 



OU W est donnee par (p.2D . Comme precedemment, les degres de libertes independants. 



au nombre de A^^ sont :0</,m,n<A^ — 1. Les relations de parite sont 

On voit que la quantite d'information non-triviale contenue dans ces relations varie selon 
qu'aucun, un, deux ou trois des indices I, m, n sont nuls. Comme dans le cas uni-dimensionnel, 
on s'attend a ce que cela se repercute sur les largeurs des gaussiennes correspondantes. En 



effet, on voit que (cf. Eq. ( p.8|) ) 
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On verifie aisementjj que la generalisation de I'Eq. ( |D.6|) est 

N-l / N-l \ ' 

i,j,k=0 \ l,m,n=0 j 

ou le prime signifie que, dans la somme, les termes dont un des trois indices m, n est nul 
sont comptes avec un facteur 1/2 supplement aire, ceux dont deux des trois indices sont 
nuls avec un facteur 1/4 et celui dont les trois indices sont nuls, avec un facteur 1/8. II est 
alors facile de voir que la distribution des est, dans le cas gaussien. 



— 1 / \ 

pm= n ^=^-p 



avec cT^ = af = {H') (1 + ^^o) (1 + Smo) (1 + S^. 



•nO)- 



Calcul numerique des coefficient /i^ 

Dans le Chap. |3|, nous devons calculer les coefficients de Fourier (ffi et (Eq. ( |D.10| )) 
des composantes d'isospin du champ et de sa derivee temporelle (f), a different instants. En 
pratique, nous nous limitons aux composantes (n, 0, 0), (0, n, 0) et (0, 0, n), qui se calculent 
comme des transformees de Fourier uni-dimensionnelles (cf. Eq. (|D.l| )). Par exemple 

2 ~ 



N 

i=0 



ou 



4>i = Y "^^'J^ • 
j,k=0 

Le calcul de la transformee de Fourier uni-dimensionnelle ci-dessus se fait aisement a I'aide 
d'un programme de type "Fast Fourier Transform" 



^Faisons le calcul dans le cas bi-dimensionnel : on ecrit 

iV-l N-l N-l 

m,n— i,j—0 fc,Z— 

= hIj + j^YY ^''^ ^''■^ + j:fJlJl -^"j + ^J1J1 ^^'^ ^f':' ■ 

i,j i,j k 'i'j I i.j k,l 

N — 1 

^ 2 ^ ^ 2 ^ ^ 1 ^ ^ 1 

On a, de meme, ^ /i™,o = YY ^^'^ + YY ^^'^ ^''■'^ YJl ^''^ ^ 4 ^0.0, 

d'ou Ton tire 

1 1 2 1 2 1 '^"^ 2 1 2 

i,j I m— m — 1 

1 1 2 1 2 

De la meme fagon on obtient — Hij Hk,j = — feo,n + ^ ^0,0 • En combinant le tout 

i,j k n — 1 

JV-1 , JV-1 , JV-1 



i,j m,n—l 

CQFD. 



2 ^ ■ 2 

m— 1 n—1 
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Annexe E 



Modele gaussien 



II est instructif de calculer la distribution de la fraction / de pions neutres dans le cas 
ou les cp^ = (p{k,tf) et (p^ = (p{k,tf) sont des variables aleatoires gaussiennes, indepen- 
dantes, de moyennes nulles et de dispersions respectives <ti et ct2 (celles-ci sont les memes 
pour toutes les direction d'isospin j = 1,2,3 et pour tous les modes k). On ecrit alors la 
probabilite (c = j, k) 

Proba {Wj:}, {^'^j) = JJ ^A^Pc) P-fii^c) dipc dipc , (E.l) 

c 



Les modes etant independants les uns des autres, nous focalisons notre attention sur I'un 
d'entre eux et omettons I'indice k correspondant. D 'apres I'Eq. Q du Chap. |, on a, 
pour cliaque composante d'isospin, 

[2 

ifj = Y — Re Oj = Aj cos 7j , 

(fj = \/2ujlmaj = LvAj siwyj , (E.3) 



ou 7j est defini par la relation Uj = y/nje^^^ et Aj = y^2nj/a;. Les differentes directions 
d'isospin etant independantes, la distribution de probabilite pour les amplitudes Aj et les 
phases 7j est le produit des distributions individuelles PA,7(^j,7j), ou 

Pa,7(^,7) = ioAP^{Acos-f) P^iujAsin-y) . (E.4) 

On en deduit la distribution de probabilite pour les valeurs de I'amplitude dans une direc- 
tion j donnee : 



Pa{A) = / djPA,^{A,j) 
Jo 

Ae-^'^+ /''"rf^g-A2A_cos(27) 

Jo 



27r(7icr2 

Ae-"^^^+/o(^'A_) (E.5) 



0-1(72 
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Annexe E. Modele gaussien 



ou luix) est la fonction de Bessel modifiee de premiere espece, et 

1 



2 2 



(E.6) 



Calcul de la distribution de la fraction de pions neutres 



La distribution de probabilite de la fraction neutre / = ^3^/(^1^ + + ^3^) est 
donnee par 

Pfif) = dxdydz Pa{x) PA{y) Pa{z) S (f - ^ j\ ^ ^ (E.7) 

En introduisant les coordonnees spheriques {r,9,(j)), ou < r < +00, < < 7r/2, et 
< (j) < 7r/2, et en posant u = cosO, on peut integrer la fonction delta {d{f — v?) = 
^{u - V7)/2\/7 pour < u < 1). On obtient 

^/(/) = o (1-/)/ drr5/o(r2/A_)e-^'^+<I>(r2(l-/)A_) , (E.8) 

2 \a1a2J Jo 

avec 

^(•2) = / sin (/) COS </> /o(-z cos^ (/)) Io(z sin^ 0) 

Jo 
1 

= — dxIo{z - x) Io{x) . 

Jo 

Cette derniere integrale est un produit de convolution au sens de la transformee de Laplace. 
En utilisant le theoreme de convolution]^, et la transformee de Laplace de la fonction de 
Bessel Iq 

£{/o(t)}(s) = / dte-''lo{t) = ^== , s>l , (E.9) 
Jo Vs'^ - 1 

on obtient 

^ / X 1 ^ 1 f 1 1 / X sinh z 



2z [s^-lj 2z 
On ecrit alors la distribution de probabilite chercliee sous la forme 

Pjif) = Fi(A+, /A_, (1 - /)A_) , (E.IO) 

ou Ton a defini les fonctions {x = r^) 

"+00 



Fp{a,b, c) 



/ dx x^ e~"-^ loibx) sinh(ca;) 
Jo 



qui satisfont de maniere evidente a la relation de recurrence 

d 

Fp{a,b,c) = -— Fp_i(a, c) . 



^La transformee de Laplace du produit de convolution de deux fonctions F et G est egal au produit des 
transformees de Laplace de ces deux fonctions : 



: dzF{t - z) (.) = C{F{t)}{s) C{G{ms) 
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Le calcul de Fo(a,6, c) est aise : en explicitant la fonction smh(x) sous I'integrale, on 
voit que 

Fo{a, b, c) = ^£{/o(6x)}(a - c) - ^C{Io{bx)}{a + c) , 



d'ou on deduit, en utilisant (E^ 



Fo(a,6,c) = - 



(a-c)'-fe2 ^(a + c)'-62 



a lb c > b . 



(E.ll) 



II vient 



Fi{a,b, c) 



a — c 



a + c 



\ 



(a - c)^ - 62 



3/2 



(a + c)^ - 62 



3/2 



o lb c > 6 . 



(E.12) 



Finalement (cf. Eq. ( [E.1Q| ), notez que pour a = A+, 6 = /A_ et c = (1 — /)A_ on a bien 
a lb c > 6), 



avec 



et 



Pf{f) = -{Fnif) + F^nif)) 



Fnif) = in-il-f)) 



n + 1 



3/2 



n 



1 9 9 

(Jo — w^ur 



(E.13) 



(E.14) 



(E.15) 



Calcul de la distribution du nombre total de pions 

De la meme fagon, on calcule la distribution du nombre total de pions dans le mode k : 



Pn{n) = - Pa^ (-IT' 

LO \UJ 



(E.16) 



avec 



f + OO 



P^,{A^)= [ dxdydzPA{x)PA{y)PA{z)5{A^ -x^ + + z^) (E.17) 
Jo 

Comme precedemment, on integre la fonction delta en passant en coordonnees sheriques. 



En utilisant (E^) il vient (n = 0056* = ^2) 
oil 



8A_ \a1a2) 



^2 e-^^+ ^(^"A_) , 



(E.18) 



z Jo 



dx Iq{zx) sinh(2;(l — x)) 



dyh{y) sinh(z - y) 
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En utilisant le theoreme de convolution de la transformee de Laplace, ainsi que les trans- 
formees de Laplace de Iq (Eq. ([E.9|)) et 



£{smht}(s) ^ 



on obtient 

Finalement, la distribution du nombre total de pions dans le mode k est donnee par 
I'Eq. ( ^l6| ), avec 

PMA') = A'e-^'^-h{A'^-). (E.19) 



Cas particuliers 

II est instructif d'etudier la formule ( [E.13D pour quelques exemples particuliers de la 
classe d'ensemble statistiques gaussiens ( [E.l| )-( p^ , et en particulier les deux cas limites 
$7 ^ 1 et 17 ^ +00. 

Le cas f2 = 1, ou encore A+ = A_, correspond aux deux situations ai = ou bien 
£72 = 0. Cela signifie qu'un des deux vecteurs, tp ou est fixe, et done nul. Le vecteur ol 
est alors proportionnel a un vecteur oriente aleatoirement dans I'espace d'isospin et on doit 
retrouver la loi en 1/ \ff pour la distribution de la fraction neutre /. En effet, considerons 
par exemple la limite 02 0, c'est a dire = 0. Alors a oc et la distribution de la 
fraction neutre tends ver^ (17 — > 1) 



Pfif) 



1 



2^7 



(E.20) 



Un deuxieme cas interessant est la sous-classe d'ensembles tels que o"2 = cjcti ($7 ^ 
+00). En termes de la distribution de probabilite Pa{oLj) pour les nombres complexes aj 
(cf. Eq. (PI) du Chap. |) 



on a, d'apres les Eqs. (|E.3|) 



Pa{a) (fa = PA,-yiA,'y) dAdj , 

(voir aussi ( [E.5| )) 
1 



Pa{a) 



(E.21) 



ou 



^Dans cette limite, Aa 



(^1 o"2 = (n) = J d^ a Pa{ct) \a\ 



comme il se doit, vers la distribution (E.2) 



Pa{A) 



/4cr2 +00, mais A+ — A_ = l/2ai. La distribution (E.5) tend. 



exp 



'2a 



ou I'on a utilise le developpement asymptotique Ii/(z) 
vient du fait que A = \ip\. 



I2^z 



Ci-dessus, le facteur 2 supplementaire 
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(|ap = n est le nombre moyen de quanta dans I'etat coherent \a) et (n) est la moyenne 
statistique de ce nombre dans I'ensemble (|E.2l|) ). La fonction Pa{ot) est appellee la represen- 
tation de Glauber Q de la matrice densite^. Pour la classe d'ensembles ( [E.21 ) la distribu- 
tion de la fraction neutre / est 



Py(/) = 2(1-/) 



(E.22) 



Cette distribution ne depend pas de la largeur a. Un exemple particulierement interessant 
de la distribution ( E.21J ) est celui d'un ensemble de quanta de frequences u) en equilibre 
thermique [^0. 

Enfin, un troisieme cas particulier interessant est celui de I'ensemble initial du le modele 



de Rajagopal et Wilzcek [jT^], utilise au Chap. y. Les variables 4^j{x) et 4'j{x) sont des 
nombres gaussiens de variances respectives {(f)^) et distribues independamment sur 

les noeuds du reseau. Cela implique (cf. Annexe |D|) que les composantes de Fourier 'fj{k) 
et ^j{k) sont des nombres gaussiens independants de variances respectives (en unites du 
pas du reseau a) 



0-1 



ou M est une constante de normalisation dependant de k. La distribution de la fraction 
neutre dans chacun des modes k est alors donnee par les Eqs. ( [E.13D -( ^l5|) avec 



(E.23) 



La Fig. [E.l] represente la distribution de la fraction neutre fk dans les modes k = {k = 
nAk, 0, 0), dans I'ensemble initial du scenario du trempage, pour les deux valeurs extremes 
de la fenetre d'impulsion etudiee au Chap. |3| : < A: < 150 MeV (0 < n < 15). 



P(f) 



1 - 



n=0 



P(f) 



0.2 0.4 0.6 0.8 1 
f 



n=15 




0.2 0.4 0.6 0.8 1 
f 



Fig. E.l: Distribution de la fraction fk de pions neutres dans le mode k ~ {k ^ nAk , 0, ) 
(Ak~ n/Na « 10 MeV), pour n — {a droite) et n = 15 (a gauche), dans I'etat initial (Eqs. (E.l)- 



^On peut toujours representer la matrice densite sur la base des etats coherents. Une tres grande classe 
d'operateur densite peuvent s'ecrire p ~ J d^a Pa{a)\a){a\. 

*0n a alors = (exp(a;/r) — ~ T/lo, la derniere egalite etant valable pour les modes de basse 
frequence {to <C T) pour lesquels rapproximation de champ classique est justifiee. 
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Remarque a propos des conditions aux bords periodiques 

Nous reprenons ici la remarque du bas de la page concernant rinconvenient des 
conditions aux bords periodiques, et en donnons une illustration. 

Avec des conditions periodiques, les sont des vecteurs d'isospin complexes : = a+ib. 
II est facile de montrer que dans I'etat initial, les composantes aj et bj des parties reelle et 
imaginaire sont des variables aleatoires gaussiennes independantes et de meme variance^. 
En d'autres termes, la distribution de probability du vecteur if est une loi du type (|E.21| ). 



Ceci implique que la probability pour que le carre du cosinus de I'angle 9 entre le vecteur 
reel \ip\ et I'axe vrs soit compris entre 003^6* = f et f + df est egaleQ a 2(1 — f)df (cf. 
Eq. ( [E.22| )). Ceci est un artefact des conditions aux bords periodiques : les parties reelles 
et imaginaires sont deux degres de liberte independants decrivant le meme mode k. Dans 
le cas ou les orientations dans I'espace d'isospin sont toutes equiprobables, on attend la 
distribution caracteristique en l/y/f- Les conditions de Neumann sont done mieux adaptees 
a I'etude de I'orientation dans I'espace d'isospin. 



'Excepte bien stir, pour le mode zero qui, lui, est reel : (p*: — 

®La loi en 1/^7, obtenue avec des conditions periodiques par I'auteur de la Ref. |^], n'est valable que 
pour le mode zero, reel. Attention, le / dont il est question ici ne represente pas la fraction de pions neutres, 
mais I'orientation du vecteur \<fi\. 



Deuxieme partie 



VERS LA FORMATION D'UN 
PLASMA DE MATIERE 
DECONFINEE 



Chapitre 4 

Equilibration thermique des gluons 



Un des enjeux principaux de I'etude des collisions nucleaires a haute energie est la 
possibilite de former un plasma de matiere liadronique deconfinee, c'est a dire un systeme 
de quarks et de gluons en equilibre thermodynamique (local). Les donnees accumulees 
durant les dix dernieres annees a TAGS (BNL) et au SPS (CERN) fournissent diverses 
indications ||l| en faveur de cette hypothese, mais il est encore difficile de parler de preuves. 
L'avenement des collisions a tres haute energie, actuellement realisees a RHIC (BNL), et, 
a partir de 2005, au LHC (CERN), ouvre une nouvelle ere dans ce domaine de recherche. 
En particulier les collaborations PHENIX a RHIC et ALICE a LHC, prevoient d'etudier 
les photons et les paires de leptons de grande energie transverse (jusqu'apt ~ 10 GeV), qui 



sont des signatures directes de la matiere deconfinee 76 1. 

Dans la region centrale d'une collision d'ions lourds a suffisament haute energie, on 
s'attend a ce que la majeure partie de I'energie transverse soit produite sous forme d'un 
grand nombre de partons de relativement grande impulsion transverse (pf ~ 1 — 2 GeV). 
Cependant qu'il est rapidement dilue par la forte expansion longitudinale, le systeme ainsi 
produit tend a s'equilibrer localement du fait des interactions mutuelles entre ses constitu- 
ants. Lorsque I'intervalle de temps moyen entre deux collisions successives devient superieur 
au temps typique de formation d'un hadron, le systeme de quarks et de gluons cesse d'ex- 
ister en tant que tel et laisse place a un gaz de hadrons en interaction. C'est une question 
d'importance considerable pour I'interpretation prochaine des donnees que de savoir si le 
systeme de partons initialement produits a le temps d'atteindre un etat d'equilibre local 
avant la phase hadronique. En effet, jusqu'a present les calculs theoriques concernant les 
signatures de la matiere deconfinee reposent sur I'hypothese selon laquelle celle-ci est en 
equilibre thermique local. 



Vers I'equilibre local 

Une premiere etude du processus d'equilibration thermique dans les premiers instants 



de la collision a ete faite par G. Baym en 1984 |77|. Dans cet article, I'auteur suppose 
que revolution du systeme de partons est gouvernee par une equation de Boltzmann qu'il 
modelise par une approximation de temps de relaxation oil le temps de relaxation 6 est 
constant. Baym montre alors de fagon analytique que malgre I'effet de I'expansion, le 



systeme atteint effectivement le regime hydro dynamique de Bjorken |10|, sur une echelle 



de temps donnee par 9. Plus tard, differents travaux consacres a I'etude des proprietes de 



transport d'un plasma de quarks et de gluons (voir par exemple |^, voir aussi |£7|)) o^^t 
permis d'estimer I'echelle de temps caracterisant la relaxation cinetique d'un tel systeme. 
Le resultat obtenu est typiquement de I'ordre de 1 fm, ce qui donne a penser que I'equilibre 



92 



Chapitre 4. Equilibration thermique des gluons 



thermique est rapidement atteint. 

Avec cette idee en tete, differents auteurs (voir en particulier |80|, et plus recem- 
ment ||8l||) etudient alors le processus d'equilibration chimique dans les collisions d'ions 
lourds. Le systeme, en expansion longitudinale rapide est suppose localement isotrope a 
chaque instant (equilibre thermique local), et des equations de transport sont derivees, qui 
decrivent revolution temporelle des variables macroscopiques comme les densites moyennes 
d'energie et de particules (gluons, quarks et antiquarks) par unite de volume. Ces equa- 
tions sont resolues numeriquement avec des conditions initiales directement tirees de codes 
Monte Carlo simulant la production de partons (minijets, voir plus loin) dans les tout pre- 
miers instants de la collision. Les conclusions de ces etudes sont que le systeme n'atteint 
pas I'equilibre chimique, les densites de quarks et de gluons sont inferieures a leurs valeurs 
a I'equilibre. Ceci a des implications importantes pour la phenomenologie, ainsi que pour 
la theorie : d'une part il faut tenir compte de I'aspect "hors d'equilibre chimique" dans 
les calculs d'observables ||8^, |8^, d'autre part il est necessaire de bien decrire revolution 
du systeme pour faire des predictions viables, ce qui implique une bonne connaissance des 
processus d'equilibration (ici chimique), ainsi qu'une bonne description de I'etat initial. 

Au milieu des annees 1990, certains auteurs reexaminent le probleme de la thermali- 
sation dans le meme cadre que I'etude de Baym, mais en faisant I'hypothese plus realiste 
d'un temps de relaxation dependant du temps |^, |86|. En effet, du fait de la forte 
expansion, le systeme est rapidement dilue, et les collisions sont de plus en plus rares, 
il est done naturel que I'echelle caracterisant la relaxation vers I'equilibre croisse avec le 
temp^. Dans ces conditions, il n'est pas clair que I'equilibre puisse etre atteint, le point 
clef reside dans la competition entre les effets respectifs de I'expansion et des collisions. Par 
exemple, dans un modele ou ^ cx t^, le systeme n'atteint jamais le regime hydrodynamique 
si p > 1 |85]. Dans un modele plus realiste, le temps de relaxation a I'instant t depend de 
I'etat du systeme a I'instant t qui depend lui-meme de la valeur du temps de relaxation a 
I'instant antecedent, et ainsi de suite. Contrairement au cas etudie par Baym, ou le temps 
d'equilibration est essentiellement donne par le temps de relaxation (constant), lui-meme 
determine par la nature intrinseque du systeme etudie, dans le cas oii le temps de relax- 
ation depend de I'histoire du systeme, le temps ej^ecii/ d'equilibration depend en plus d'un 
aspect extrinseque : la condition initiale. Pour estimer la duree du regime transitoire il est 
done necessaire de suivre revolution du systeme a partir d'une condition initiale donnee. 
Dans le cas des collisions d'ions lourds, qui nous interesse ici, cela signifie que Ton doit 
avoir une description realiste de I'etat initial du systeme de partons ainsi que des processus 
qui gouvernent son evolution ulterieure. 



Production de gluons dans les premiers instants de la collision : 
le probleme de la condition initiale. 

Le probleme de la caracterisation de I'etat initial du systeme de gluons produits lors 
d'une collision d'ions lourds ultra-relativistes est tres delicat et constitue un domaine de 
recherche en soi. Deux scenarios ont ete proposes dans la litterature : le scenario des 
minijets et le scenario de saturation. Dans ce qui suit nous les decrivons de maniere tres 
schematique, notre but etant de donner une image intuitive de ces scenarios. 

Plagons-nous dans le referentiel du centre de masse de la collision. Avant I'impact, 

^Par exemple, au voisinnage de I'equilibre local, la seule quantite dimensionnee est la temperature 
locale T, et le temps de relaxation est done inversement proportionnel a T. Cependant, pour un systeme 
en expansion longitudinale dans le regime hydrodynamique, T oc t~^^^ jlo[, et done 6 oc t^^^. 
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chacun des noyaux incidents pent etre vu comme un disque tres aplati dans la direction 
de son mouvement, a I'interieur duquel sont confines les quarks de valence des nucleons, et 
avec lequel se deplace un "nuage" de partons virtuels : les quarks de la mer et les gluons, 
les plus nombreux etant ceux dont I'impulsion longitudinale est tres faible comparee a 
celle des noyaux incidents (la fraction x d'impulsion longitudinale emportee par le parton 
est <C 1). Au moment de la collision, les deux noyaux se traversent I'un I'autre, et les 
partons de leurs nuages partoniques respectifs interagissent les uns avec les autres. Du fait 
de ces interactions, certains de ces partons virtuels (mais quasi-reels) sont "projetes" sur 
leur couche de masse, on dit qu'ils sont "liberes" lors de la collision : ils ne font plus partie 
integrante des noyaux qui les "transportaient" initialement. Un grand nombre de partons 
reels d'impulsion transverse pt sont ainsi produits pendant un intervalle de temps ~ l/pt, 
qui est aussi I'echelle de temps typique pour que les partons d'impulsions longitudinales 
differentes se separent pliysiquement les uns des autres [37, 88|. 



Le scenario des minijets 

Si I'impulsion typique du parton produit est sufRsamment grande devant I'echelle typ- 
ique des interactions fortes {pt ^ ^qcd ~ 200 MeV), le processus de production releve 
du regime perturbatif {as{pt) <C 1), ou le taux de production peut etre calcule de fagon 
controlee (voir Fig. |4.1| ). La section efficace de production decroissant comme 1/pj (voir 
par exemple [^]), la contribution dominante a I'energie transverse initialement produite 
vient des gluons d'energie intermediaire (sufRsamment faible pour que la probability de 
production soit considerable, et suffisamment grande pour que le calcul perturbatif ait un 



sens), ptr^l-2 GeV : les minijets |90, 87, 39|. 



MINUET 




Fig. 4.1: Diffusion elastique entre deux partons des noyaux incidents lors d'une collision nucleaire 
a haute energie. Le taux de production perturbatif est domine par les partons d'impulsion pt ~ 
1 — 2 GeV, les minijets. 
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Chapitre 4. Equilibration thermique des gluons 



Le scenario de saturation 



Le scenario de saturation propose une origine differente pour la production initiale 
d'energie transverse. Dans les noyaux incidents, la densite de gluons de petit x par unite 
d'espace des phases transverse ^i^^p^ , ou h mesure la position dans le plan transverse, 
augmente quand pt diminue, et atteint un maximum pour des impulsions transverses pt < 
Qs, oii Qs est I'impulsion de saturation |^, Les gluons de faible impulsion transverse 

etant done en tres grand nombre, on s'attend a ce que leur contribution a la production 
d'energie transverse soit non-negligeable. Cependant, le calcul perturbatif n'a pas de sens 
dans ce cas, et il n'existe a ce jour aucun calcul solide du taux de production de ces gluons. 
On peut cependant se faire un image du processus de production : les gluons consideres 
ici etant tres nombreux, on peut les voir comme formant un champ classiqueQ genere par 
les quarks de valence des nucleons, en mouvement rapide dans le referentiel considere : 



c'est I'equivalent non-abelien du champ de Weizsacker- Williams ||9J]. Dans cette image, la 



collision nucleaire est vue comme une "collision" entre les champs associes a chacun des 



noyaux incidents. Sur la Fig. 4.2, on a represente tres schematiquement cette collision : 
un gluon virtuel d'un des noyaux est libere (projete sur couche de masse) par interaction 
avec le champ classique de I'autre noyau. En pratique, on peut decrire cette collision en 



rfesolvant (numeriquement) les equation de champ classique ||95| , 96|. On peut se convaincre 
avec A. H. Mueller |Q, que la majorite des gluons liberes sont ceux dont I'impulsion 
Pt ^ Qs, et ce durant un intervalle de temps ~ La valeur de I'impulsion de saturation 

Qs depend de la fraction x d'impulsion longitudinale des gluons consideres, et done de 
I'energie de la collision. Aux energies de RHIC (a/s = 200 AGeV), Qs ^ I GeV, et a LHC 
(Vs = 5.5 ATeV), 2 - 3 GeV. 



B 



Fig. 4.2: Dans le scenario de saturation, les gluons virtuels de petite impulsion transverse des 
noyaux incidents peuvent etre vus comme un champ classique de couleur associe aux quarks de 
valence en mouvement rapide. Dans la representation schematique de la collision, ci-dessus, on 
suit un gluon du champ du noyau provenant de la gauche. Celui-ci est libere lors de la collision, 
en interagissant avec I'ensemble des gluons de I'autre noyau, c'est a dire avec le champ classique 
associe. 



^Exactement comme nous I'avons fait pour le champ de pions dans le chapitre precedent. 
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Dans ces deux scenarios la contribution dominante a la production initiale d'energie 
transverse provient des partons dont I'impulsion transverse est de I'ordre de 1 — 2 GeV. 
Cependant, I'origine de ces partons est tres differente dans les deux cas, ce qui conduit a 
des predictions tres differentes pour I'etat initial du systeme. La Fig. ou on a represente 
schematiquement la densite de gluons produits par unite d'espace des phases transverse, 
resume ces deux scenarios. 

SATURATION : 
(Champ classique) 




1-2 GeV Pt 

Fig. 4.3: Les densites de gluons produits par unite d'espace des phases transverse dans les scenarios 
de saturation et des minijets. Dans les deux cas la region qui domine la production d'energie 
transverse est la meme : ce sint les gluons d'impulsion transverse pt ^ I ~ 2 GeV. L'origine de ces 
gluons est cependant tres differente dans les deux cas, ce qui conduit a des predictions differentes 
pour I'etat initial du systeme. 



Recemment, les auteurs des Refs. I^S], 97, 98| d'une part et d'autre part ont etudie, 
independamment et dans des approches differentes, la question de la thermalisation du gaz 
de gluons respectivement dans les scenarios de saturation et des minijets. Dans |^, |97|| , 
A. H. Mueller considere I'effet des collisions elastiques dans I'approximation des petites 
deviations. En supposant une forme idealisee pour la distribution initiale des gluons, il est 
capable de suivre analytiquement revolution du systeme durant les premiers instants de 
la collision. Bien que ses approximations ne lui permettent pas de suivre le gaz de gluons 
jusqu'a I'equilibre, il semble que le processus d'equilibration commence tres tot. Dans la 
Ref. [ p^ , J. Bjoraker et R. Venugopalan reprennent ce modele et resolvent numeriquement 
I'equation de Boltzmann proposee par Mueller. II apparait clairement que le systeme tend 
vers I'equilibre malgre I'expansion. En comparant la dependance en temps de differentes 
quantites macroscopiques (temperature et densite d'entropie par unite de volume) avec 
celle attendue dans le regime hydrodynamique, ces auteurs mesurent des temps d'equili- 
bration de I'ordre de 3 — 4 fm a RHIC et 2 — 3 fm a LHC. Parallelement, G. C. Nayak, 
A. Dumitru, L. Mc Lerran et W. Greiner ||9^ etudient la thermalisation du gaz de mini- 
jets dans une approximation de temps de relaxation et concluent a I'equilibration avec des 
temps caracteristiques de I'ordre de 4 — 5 fm aussi bien aux energies de RHIC que de LHC. 

Les approches utilisees dans |Q d'une part et dans d'autre part sont cependant tres 
differentes, et il est difficile de comparer ces resultats les uns avec les autres. En particulier 
les auteurs de |99| identifient le temps de relaxation avec le libre parcours moyen, et il 
n'est pas clair que la dynamique qui resulte de cette hypothese arbitraire soit directement 
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reliee a I'equation de Boltzmann resolue dans |^].Dans ce chapitre, nous reprenons I'etude 
de la thermalisation du gaz de gluons et comparons les scenarios de saturation et des 
minijets dans une seule et meme approche. Suivant Mueller (9^], nous considerons I'effet des 
collisions elastiques dans la limite des petites deviations (voir 1' Annexe 0). Nous traitons 
le terme de collision de I'equation de Boltzmann correspondante dans I'approximation 
du temps de relaxation et calculous le temps de relaxation de fagon auto-coherente en 
faisant directement intervenir la dynamique de I'equation de Boltzmann sous-jacente (voir 
aussi |84, 85, IM]). Dans le cas de I'integrale de collision de Landau-Mueller, nous sommes 



capables de calculer analytiquement toutes les integrales sur I'espace des impulsions, ce qui 
diminue considerablement I'effort numerique a fournir. En comparant nos resultats avec la 
solution exacte de [^], on peut tester notre approche, qui s'avere etre remarquablement 
bonne compte tenu de son degre de simplicite. 

Nous arguons que le critere utilise dans |98| pour caracteriser I'equilibre n'est pas satis- 
faisant et proposons plutot de tester le degre d'isotropie de differentes observables, critere 
plus intuitif d'une part et plus directement relie a la phenomenologie d'autre part. En con- 
sequence, nos conclusions sont qualitativement differentes de celles de |Q, en particulier 
nous montrons qu'aux energies de RHIC les collisions elastiques ne sont pas sufRsantes 
pour permettre au systeme d'atteindre I'equilibre thermique avant la phase hadronique. 
En ce qui concerne le scenario des minijets, nous arguons que les resultats de |99| ne sont 
pas credibles, les auteurs n'ayant pas correctement implemente les lois de conservation. Ici 
encore, nous montrons que I'equilibre n'est pas atteint a RHIC. 

De meme que dans les chapitres precedents, nous tenterons de mettre en avant les 
idees physiques, les details techniques etant relegues en Annexe. Le travail presente dans 
ce chapitre a ete realise en collaboration avec Dominique Schiff ||1Q0||. 



4.1. Le modele 
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4.1 Le modele 



Considerons une collision entre deux noyaux lourds a tres haute energie. Pour fixer les 
idees, choisissons deux noyaux identiques, spheriques de rayon R. Dans le referentiel du 
centre de masse, du fait de la contraction de Lorentz, chacun de ces noyaux a I'aspect 
d'une galette de rayon R et d'epaisseur / = R\/l — dans la direction z de I'axe de la 
collision, ou /5 < 1 est la vitesse du noyau en unite de la vitesse de la lumiere. Dans une 
collision a tres haute energie, on peut modeliser les deux noyaux incidents par deux disques 
infiniment plats (/? ~ 1, / ~ 0) et de rayon infini]^ Apres la collision, la matiere baryonique 
etant tres peu ralentie les deux disques, qui s'eloignent maintenant I'un de I'autre 
avec des vitesses proches de celle de la lumiere, laissent derriere eux une region de matiere 
non-baryonique hautement excitee. Nous adopterons I'idealisation de Bjorken [||, [l^, ou 
les deux disques s'eloignent a la vitesse de la lumiere. Le probleme est alors invariant sous 
les transformations de Lorentz suivant I'axe de la collision. 

Nous considererons ici des collisions a tres haute energie (RHIC,LHC), et nous sup- 
poserons que la constante de couplage des interactions fortes 05 ^ 1. Durant les tons 
premiers instants apres I'impact, un grand nombre de partons, essentiellement des gluons, 
sont produits. L'expansion dilue tres rapidement le systeme, et les nombres d'occupation 
(les densites de partons par unite d'espace des phases : / oc d^N/SpSx) deviennent rapi- 
dement inferieurs a I'unit^. On peut alors traiter le systeme comme un gaz de particules 
classiques0, et decrire le taux de variation par unite de temps de la distribution / par une 
equation de Boltzmann : 



ou le membre de gauche designe la variation due au mouvement libre des particules du 
gaz : en I'absence de champ externe, ^ = dtf{p,x,t) + v ■ Vxf{p,x,t), avec v = p/p. Le 
membre de droite decrit I'effet des collisions entre les particules du gaz. 



La distribution f{p, x, t) decrit la densite de gluons reels (sur couche de masse : = 0) 
par unite d'espace des phases a I'instant t, et depend en general de I'etat de couleur 
des gluons consideres. Ici, nous cherchons a decrire la relaxation vers I'equilibre cinetique 
(local), c'est a dire de relaxation vers une distribution (localement) isotrope. L'echelle de 
temps caracteristique de ce processus etant grande devant celle du processus de relaxation 
de couleui]^, on peut negliger ce dernier sur les echelles de temps qui nous interessent. 
Bien que la derivation du terme de collision que nous allons utiliser repose sur I'hypothese 
d'une distribution individuelle independante de la couleur, I'equation de Boltzmann obtenue 
doit etre vue comme une equation effective pour les echelles de temps caracteristiques 
de la relaxation cinetique. Dans cette logique la distribution / s'interprete comme une 
distribution moyennee sur les etats de couleur (et de polarisation). 



I'instant de I'impact {t = 0), toute I'energie de la collision est localisee dans une region de taille 
longitudinale nuUe. Cette energie etant finie, la region en question doit avoir une extension transverse 
infinie Cette idealisation est raisonnable dans la region centrale de la collision {z ~ 0, \\xt\\ <C R, ou Xt 
repere la position dans le plan transverse a I'axe de la collision) et pour des temps au moins inferieurs a la 
taille transverse du systeme : t < R. 

*0n peut se convaincre que la description en termes d'une equation de Boltzmann est valable des que 
/ £ I/qs- Quand / < 1, on peut negliger les effets de statistique quantique. 

^Par opposition a I'approximation de champ classique, valide quand les nombres d'occupation sont 
grands devant I'unite. 

®Au voisinage de I'equilibre, les echelles de temps caracteristiques des processus de relaxation de couleur 
et de relaxation cinetique sont bien separees : t~J^i^^^ ~ T~i.-^„- / ag, la premiere, etant beaucoup plus 



sensible a la physique de grande longueur d'onde (voir par exemple |101|) 
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Mentionnons de plus que la description en terme d'une equation de Boltzmann "stan- 
dard" n'a de sens que si le temps moyen entre deux collisions (le libre parcours moyen /) 
successives est sufRsamment grand devant la longueur d'onde des particules considerees 
£ ~ L'echelle typique d'impulsion des gluons decrits iciQ est done p > . 

Du fait de la geometrie du probleme, la physique est invariante sous les transformations 
de Lorentz le long de I'axe z de la collision |77t|^ 



f{pt,Pz,z,t) = f{pt,Pz,r) , 

ou r = y/fi — z"^ et Pz = "fiPz — up), avec p^ = pi +p1, u = z/t et 7 = (1 — u^)~^^'^ = t/r. 
II suffit done de se limiter a la tranche z = 0. II est facile de voir que 



dt 



dtf{p,t) -PldpJ{p,t) = dtfm\r..t=cte 



-0 t 



En ce qui concerne le terme de collision de I'equation de Boltzmann, nous nous limiterons, 
comme nous I'avons mentionne dans I'introduction, aux collisions elastiques entre gluons. 
La section efhcace differentielle de diffusion gg gg est fortement piquee vers I'avant, et 
diverge pour les diffusion a tres petit angle de deviation, exactement comme dans le cas 
de la diffusion entre particules chargees : c'est la divergence de Rutherford due a la portee 
infinie de I'interaction (I'echange d'un quantum de masse nulle). Les collisions a petite 
deviation etant les plus frequentes, elles dominent la dynamique et nous ne retiendrons que 
leur contribution a I'integrale de collision. En developpant I'integrand autour des petites 
valeurs de Tangle de deviation |97, |102| (voir aussi I'Annexe la premiere contribution 



non nulle est egale a la divergence du flux de particules dans I'espace des impulsions : 
I'effet des collisions elastiques a petite deviation est equivalent a un processus diffusif dans 
I'espace des impulsions. En explicitant ce terme, on obtient, dans la region centrale {z = 0) 
et dans 1' approximation de particules classiques (/ <C 1) (cf. Eq. ( F.14|) ), 



dtflp^, = BNoVlf + 2BN_iVpivf) , (4.1) 
ou Ton a defini les moments^] {f^ = 2{N^ - 1) / d^p/{27rf ) 



Nsit) 



{f)= !yf{p,t) (4.2) 



{No = n est la densite moyenne de particules par unite de volume, A^i = e est la densite 
moyenne d'energie par unite de volume). La divergence de Rutherford de la section efficace 



'^Les excitations de grande longueur d'onde (i" 3> ou e est I'energie moyenne par particule) sont des 
modes coUectifs (ils "contiennent" un grand nombre d'excitations dures ~ et agissent comme un 

champ externe classique (qui n'est autre que le champ moy en cree par le mouvement des particules dures) 



sur la dynamique des modes durs (voir par exemple |101|). Dans ce chapitre, nous negligeons ce champ 
externe. 

^L'echelle caracterisant les inhomogeneites spatiales dans la direction transverse est ~ R. Dans la region 
centrale on pent done negliger la dependance en xt de la distribution /. De plus, par symetrie, / ne depend 
que de pt = \\pt\\ et est paire en pz. 

^Nous introduisons ici une notation differente de celle habituellement utilisee dans la litterature : (m) = 
J d^pm{p) f{p,t) designe la moyenne par unite de volume et (m) = {m)/n la moyenne par particule 
(n = (1) est la densite moyenne de particules par unite de volume). 



4.1. Le modele 
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de diffusion gluon-gluon se manifeste|^ par la presence du "grand logarithme" L = 2 f 

Cette divergence logarithmique est regularisee par le phenomene d'ecrantage de Debye : 

(4.3) 




ou m£) est la masse d'ecran de Debye, et oil p designe I'ordre de grandeur de I'impulsion 
typique des particules du milieu {m'j^ <C p) (voir les Eqs.(^) et ( pUOD de 1' Annexe 0). 



4.1.1 Lois de conservation 

Dans la region centrale, le tenseur energie-impulsion s'ecrit 

T^- = {^)= f P!^f{p,t) ^dmg{e, Pt,Pt,Pl), 
P Jp P 

et est diagonal, du fait de la symetrie du probleme {f{p,z = 0,t) = f{—p,z = 0,t)). Ses 
composantes representent respectivement les densites moyennes d'energie et de pression 
(transverse et longitudinale) dans la region centrale a I'instant t : 

e{t) = ip) = m{t) , PrAt) = ipljp) , 

ou = = {Px +Py)/'^- Les particules du milieu etant de masse nulle (T^^ = 0), on 
&e = PL + 2PT. 

Les lois de conservation de I'impulsion et du moment angulaire sont trivialement sat- 
isfaites par symetrie. L'energie etant conservee dans les collisions individuelles, le taux de 
variation par unite de temps de la densite d'energie induit par les collisions est nul : 

pC[f]{p,t) = 0. 

p 

On en deduit I'equation d'evolution de la densite d'energie (e = dte) : 

e + i±^ = 0. (4.4) 

De plus, il est clair que les collisions elastiques conservent le nombre total de particules, ce 
qui implique que le taux de variation collisionnel de la densite moyenne de particules est 
nul. On a done I'equation (n(t) = (1) = No{t)) 



n+'^ = jc[f]{p,t) = d, 



dont la solution est 



n(t)=n(to)y. (4.5) 

Remarquons que ces deux equations, ne faisant pas intervenir le terme de collision, 
gardent la meme forme quelles que soient les approximation faites sur ce dernier. 

^"Les tout premiers termes du developpement dans I'angle de deviation 6 donnent, en principe, des 
contributions plus singulieres. Cependant, ces termes sont nuls dans le cas d'une distribution independante 
des degres de liberte de couleur. C'est la I'origine de la hierarchie des echelles de temps caracterisant 
les processus de relaxation de couleur et cinetique. Cette derniere est en fait sensible a la physique de 
I'ecrantage fchromo-) electrique, c'est a dire aux echelles de distance de I'ordre de la longueur de Debye 
(voir Eq. @). 
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4.1.2 Expansion .vs. Collisions 

Partant d'une condition initiale donnee, c'est a dire d'une distribution initiale foip), 
revolution du systeme est gouvernee par la competition entre deux phenomenes : d'un cote 
(a gauche dans I'Eq. I'expansion, de I'autre les collisions. 

L'expansion longitudinale rapide du systeme tend a rendre la distribution fortement 
anisotrope : a impulsion {pt,Pz) donnee, le nombre moyen de particules qui quittent la 
region centrale librement (c'est a dire par leur mouvement propre, sans I'aide des collisions) 
est superieur au nombre de particules qui y penetrent de la meme maniere, et ceci d'autant 
plus que I'impulsion longitudinale consideree est grande. L'expansion a done tendance a 
"vider" la region centrale de son impulsion longitudinale. 

L'effet des collisions est, schematiquement, de creer des particules d'impulsions diverses, 
d'orientations aleatoires et tend done, a I'inverse de l'expansion, a rendre la distribution 
isotrope (dans la region centrale). Dans 1' approximation des petites deviations, il s'agit d'un 
processus diffusif dans I'espace des impulsions [|102| (voir la discussion plus haut, ainsi que 
I'Annexe ^ : de la meme maniere qu'un gaz tend, par diffusion, a remplir I'espace (des 
configurations) de maniere isotrope par diffusion, on a ici un gaz qui se diffuse dans I'espace 
des impulsions, et tend a occuper celui-ci de maniere isotrope. 

II est clair que les influences respectives de chacun de ces effets sur revolution du systeme 
depend de I'etat initial de celui-ci. L'objet de ce chapitre est precisement I'etude de cette 
competition pour differentes conditions initiales. A cette fin, il est instructif de caracteriser 
les deux regimes limites correspondant chacun a la victoire d'un effet sur I'autre. 

Le cas libre 

Dans le cas ou I'influence des collisions est completement insignifiante devant celle de 
I'expansion, on pent negliger le membre de droite de I'equation de Boltzmann, ce qui est 
formellement equivalent au cas ou les particules du gaz sont libres {as = 0) : 

dtflibre{P,t)\p^t = Cl- 
ou encore 

flibre{P,t) = fo{Pt,PzP) ■ (4.6) 

to 

La distribution devient rapidement ~ 6{pz), on voit par exemple que {pD / {Pt) ~ 



Le regime hydrodynamique 

La situation a I'extreme oppose est celle ou l'effet des collisions est suffisament fort pour 
maintenir I'isotropie a chaque instant. Dans ce cas, c'est le terme de gauche de I'Eq. (4.1) 
qui est negligeable et la forme de la distribution est determinee par I'equation 

C[/e,]=0. (4.7) 

La solution isotrope la plus generale est de la forme 

/e,(p,t) = A(t)e-P/^W, (4.8) 

ou les fonctions du temps X{t) et T{t) (fugacite et temperature locales) sont determinees 
par les equations de conservation qui, ne faisant pas intervenir le terme de collision, gardent 
leurs formes respectives ( |4.4| ) et ( ^.5] ) dans la limite (^.Tj). On a {Pl = Pt = ^/S) 



4.2. L'approximation du temps de relaxation 
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d'ou on deduit, pour la fugacite et la temperature dans le regime hydrodynamique, 

Kydro ~ Cte , Thydro ~ i^^^'^ ■ (4.9) 

Sur une echelle de temps At <^ O^eiax-, ou ^re/ax est Techelle de temps caracterisant 
le processus de relaxation vers I'equilibre local, la dynamique est celle du regime libre. A 
I'inverse, sur des echelles de temps At » Oreiax, la dynamique est decrite par le regime 
hydrodynamique. Enfin, pour des echelles de temps intermediaires At ~ Oreiax, le systeme 
est dans un regime transitoire, la dynamique est decrite par I'Eq. (O). Dans la suite, 



I'echelle de temps a laquelle on s'interesse est donnee par la duree de vie du systeme de 
partons, lequel cesse d'exister en tant que tel quand la densite d'energie moyenne par unite 
de volume devient inferieure a la densite critique ec ~ IGeV/fm^. La question est de savoir 
si le systeme atteint le regime hydrodynamique avant la fin du temps qui lui est imparti. 
Autrement dit, nous voulons savoir comment se compare la duree du regime transitoire 
avec la duree de vie. 

Dans une situation statique, I'ordre de grandeur de la duree du regime transitoire est 
simplement donne par Oreiax- Dans le cas present cependant, la question posee est moins 
triviale qu'elle n'en a I'air. En effet, du fait de I'expansion, le nombre de particules par unite 
de volume diminue et les collisions sont done de moins en moins frequentes. Autrement dit, 
I'echelle de temps Oreiax augmente avec le temps : le systeme approche de I'equilibre local 
de plus en plus lentement a mesure que le temps passe. La duree du regime transitoire 
"effectif" depend done de I'histoire du systeme (et done aussi de son etat initial) et ne 
peut etre obtenue qu'en suivant son evolution au cours du temps, c'est a dire en resolvant 
I'equation de Boltzmann (|4] 



4.2 L'approximation du temps de relaxation 



La resolution numerique de I'Eq. (^Jj) [£8| s'avere relativement lourde en regard de la 
question posee, essentiellement qualitative (le modele tres simple, decrit precedemment, ne 
peut etre credible au niveau quantitatif). Pour cette raison, nous allons encore simplifier 
la dynamique en remplagant le terme de collision par un terme lineaire dans la difference 
/ — /eg. Autrement dit, nous modelisons I'effet des collisions sur le taux de variation de la 
distribution / par unite de temps par un terme de relaxation exponentielle sur une echelle 
de temps : le temps de relaxation (voir par exemple ||7^) : 



0{t) 

ou (cf. Eq. 



dtnp,%^t = -^^^'\^^, (4.10) 



feg{p,t) = X{t)e-P/^^'\ 
les fonctions du temps A et T devant etre determinee a I'aide des lois de conservation^ 

(il)-(ii). 

L'approximation de temps de relaxation consiste a parametriser la distribution /, so- 
lution de I'equation de Boltzmann, a I'aide des trois quantite A(t), T(t) et 0(t), et ce a 
chaque instant t. On ne doit pas s'attendre a ce que la distribution soit bien decrite par 
cette parametrisation. Par exemple, dans cet ansatz on voit que toutes les impulsions re- 
laxent vers I'equilibre avec le meme taux, ce qui n'a a priori aucune raison d'etre le cas. En 

^^Comme nous I'avons vu plus haut, dans le cas oil seules les collisions elastiques sont prises en compte, 
le nombre de particules est conserve. II est alors necessaire d'introduire le parametre A pour prendre en 
compte cette loi de conservation, ce qui a ete omis par les auteurs de la Ref. [B9|. 
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ce sens, I'approximation de relaxation est une sorte d'approximation de champ moyen : les 
differentes impulsions sont effectivement decouplees les unes des autres et relaxent toutes 
avec un meme taux moyen que nous determinerons dans la suite de fagon auto-coherente. 



4.2.1 La methode des moments 

On pent calculer 9 en identifiant le taux de variation collisionnel de certains moments 
de la distribution / avec son analogue dans I'ansatz ([4.10| ) (voir par exemple |85, ^] ). 



Le temps de relaxation ainsi obtenu contient I'information concernant I'effet des collisions 
sur la dynamique. Nous appliquons ici cette methode de maniere systematique pour dif- 
ferents moments caracterisant differentes echelles d'impulsion. Le choix du moment le plus 
adequat depend de la physique a laquelle on s'interesse. Dans le cas du terme de collision 
de I'Eq. ( |4.lD , nous pourrons calculer analytiquement toutes les integrates sur I'espace des 
impulsions, ce qui simplifie consider ablement la resolution numerique. 

Prise au pied de la lettre, I'approximation du temps de relaxation consiste a identifier 
le terme de collision de I'equation de Boltzmann avec le membre de gauche de I'Eq ( [1.1CI| ) : 



/ /eg 



Cette equation, prise au sens fort, c'est a dire telle qu'elle est ecrite ci-dessus, pour tout p, 
est trop contraignante et n'a pas de solution. On doit en fait la prendre dans un sens plus 
faible, en en prenant la valeur moyenne sur les impulsions avec un certain poids : 



(m) — (m) 



eg 



[ m{p)C[f]ip,t), (4.11) 
Jp 



ou 



(m)(eg)(t) = / m(p)/(,g)(p,f) (4.12) 
Jp 

la fonction m(p) etant a priori quelconque, et oil 9m est le temps de relaxation corre- 
spondant. Nous voyons encore une fois que I'approximation du temps de relaxation ( 4.10| ) 



n'est pas, a strictement parler, une approximation pour la dynamique microscopique, mais 
plutot une modelisation de la relaxation du systeme dans sa globalite. 
Les equations de conservation conduisent aux relations]^ 

_ 1 

n{t) = neq{t) = 2^^\{t)T^{t), (4.13) 

e{t) = eeg{t) = 6^^^^X{t)T\t). (4.14) 

A ce point, il est important de faire la remarque suivante : les parametres A et T n'ont pas, 
en general, les significations physiques de la fugacite]^ et de la temperature du systeme, 
au sens thermodynamique (ou plutot hydro dynamique dans notre cas) de ces termes. En 
effet, les notions de fugacite et de temperature (locales) ne sont bien definies que pour 
un systeme a I'equilibre (local), auquel cas elles s'identifient avec les parametres A et T : 

^hydro at T Tfiy^i^Q. 
i2ivr(t)= j p=/,^(p-,i) = (s + 2)!^^i^A(t)r=+3(t). 

J p ^ 

^^On peut ecrire A = exp(— ///T). A I'equilibre (local), = ^hydro est le potentiel chimique (local). 
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4.2.2 Procedure de resolution 



Les parametres inconnus 9, X et T sont obtenus, a chaque instant, a partir des Eqs. 
( 4.111 ), ( [1.13|) et ([4.14 ). Illustrons la procedure de resolution sur un exemple simple : 



ou p\ = Pi /2. On a done (m)eq = 0, et, apres quelques integrations par parties, on obtient, 
pour I'Eq. (pl|). 



{pi - p\) 



4BN^i{Pl-Pt). 



(4.15) 



II nous faut done calculer les quantites (nous laissons de cote le calcul de B pour le moment) 
{p1 ~ Pt/'^)^ -^-ii Pl — Pt, et e a chaque instant t. Pour ce faire, nous ecrivons la solution 
formelle de I'Eq. (p^) || 



fip,t) 



avec 



/o(pl,P.f)e-^W + 

^0 



xit) 



to 



gx(t')~a-(t) ^ 

feq (Pt .PzJ,,t') 



dt' 



(4.16) 



no 0{t') 

Pour une distribution initiale donnee, on en deduit, pour chacune des quantites pre-citees, 
une equation du type 



M{t) = Mito) 4?(to/t) + / dt'^ :F^M'\t'/t) M,,{t') , 



to 



(4.17) 



ou X = x{t), x' = x{t'), 9' = 9{t') et M = {pD, {p]), iV_i, Pl, Pt, ou e, selon le cas. Les 
fonctions T^^j^'^^ ne dependent que de la geometric des distributions /o et feq dans I'espace 
des impulsions, elles sont donnees dans 1' Annexe |G| pour les differents M dont nous aurons 
besoin dans ce chapitre. Les moments M^q, calcules avec la distribution /eg, ne dependent 
que de A et T. 

On resoud numeriquement le systeme d'equations ( 4.13 ), ( 4.14 ), (|4.15D , ( 4.17 ) par itera- 
tions successives, selon la methode proposee dans la Ref. Dans le cas present, toutes 
les integrales sur les impulsions peuvent etre calculees analytiquement. Les integrales tem- 
porelles (cf. Eq. ( [4.17| )) sont calculees par la methode de Simpson (corrigee par un trapeze 
dans les cas ou le nombre de points sous I'integrale discrete est pair) |64|. A chaque pas 
de temps, la procedure d'iteration (decrite dans |Q) est initialisee par une premiere es- 
timation des quantites M, obtenue en calculant les integrales temporelles par la methode 
des rectangles, pour laquelle il n'est pas necessaire de connaitre la valeur de I'integrand 
pour la borne superieure de I'integrale. Cette procedure d'initialisation n'est pas stricte- 
ment necessaire, mais ameliore la convergence de la procedure d'iteration. On verifie la 
bonne precision de la methode en resolvant, en parallele avec notre systeme d'equations. 



I'Eq. (|4.17|) dans des cas solubles analytiquement, par exemple en calculant n(t), non pas 
a I'aide de I'Eq. (^^), mais a partir de I'Eq. ( 4.17 ) avec Tn{a) 



a. 



4.2.3 Le choix du moment 

Comme nous I'avons remarque plus haut, dans l'approximation du temps de relaxation 
la distribution / relaxe vers la distribution feq d'une maniere independante de p : toutes 
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les echelles d 'impulsions relaxent avec le meme taux moyen 9 (cf Eq. ([4.1QD ). D'un autre 
cote, nous avons vu que le calcul de ce temps de relaxation necessite de choisir un moment 
particulier de la distribution, c'est a dire une fonction m{p) (cf. Eq. ( 4. 111 )). II est clair 
que les differents moments (m) sont sensibles a differentes parties de la distribution /, 
c'est a dire a differentes echelles d'impulsion, et il en est par consequent de meme pour les 
temps de relaxation 0m correspondant. Ceci est illustre sur la Fig. 4^, ou on a represente 
revolution temporelle du rapport Pl/Pt des pressions longitudinale et transverse pour 
differents choix de m, pour une condition initiale donnee. 

e(fm) 



3 
2 

1 






5 10 

temps (fm) 



15 




5 10 
temps (fm) 



15 



Fig. 4.4: Le rapport des pressions longitudinale et transverse (a droite), ainsi que le temps de 
relaxation en fm (a gauche) en fonction du temps (en fm), pour differents choix du moment (m) : 
(mi) ^ Pl — Pt (courbes en trait-plein) et (7712) = {pi — p\) (courbes en trait-tiret). Le premier, 
etant sensible a une echelle d'impulsion plus faible que le second, correspond a une relaxation plus 
rapide vers le regime hydrodynamique. Ceci est du au fait que les impulsions les plus grandes sont 
les plus sensibles a I'expansion (voir la discussion page |100[ ). Ces courbes ont ete obtenues dans le 
scenario de saturation a RHIC (voir la section Resultats). 

On pent done preciser la signification de I'approximation du temps de relaxation : on 
remplace la dynamique microscopique par une relaxation globale vers I'equilibre, avec un 
taux caracteristique d'une certaine echelle d'impulsion. Pour un probleme donne, il est 
done judicieux de choisir le moment m selon I'echelle typique d'impulsion qui caracterise 
la physique a laquelle on s'interesse. En ce qui concerne la question de la relaxation vers 
le regime hydrodynamique, I'echelle d'impulsion typique qui nous interesse est celle qui 
caracterise les variables macroscopiques telles que la pression ou la densite d'energie. Le 
choix (m) = e (c'est a dire m{p) = p) correspondant deja a I'equation de conservation de 
I'energie, nous opterons pouij^ (m) = Pl — Pt, c'est a dire 

pI-pI 

mip) = . 

P 

Apres quelques integrations par parties, on obtient, a partir de I'Eq. ( [4.11|) , 

^^-^^ =4BNo{^i^) + 2BN^i{^^). (4.18) 



^*0n obtient des resultats identiques a ceux presentes dans la suite en c hoisi ssant (m) — Pl ou (m) — Pt. 
La difference des pressions n'est autre que la viscosite de cisaillement C, |l03|] : Pl — Pt = 2C,/'it. 
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Comme nous I'avons explique dans I'introduction, notre but est ici d'etudier de fagon 
systematique le processus de relaxation vers I'equilibre local (le regime hydrodynamique) 
pour differentes conditions initiales proposees dans la litterature pour decrire I'etat du 
systeme immediatement apres I'impact dans les collisions d'ions lourds a RHIC et LHC : le 
scenario de saturation et le scenario des minijets. Bien que cliacun de ces deux cas ait ete 
etudie recemment ||8^, 97, |9^, ils n'ont jamais ete compares I'un a I'autre dans une seule 
et meme approche. En comparant nos resultats avec ceux de la Ref. |98], ou I'Eq. (4J_) est 
resolue numeriquement dans le scenario de saturation, nous pouvons tester notre approche, 
qui s'avere etre remarquablement bonne, compte tenu de son caractere tres qualitatif. 
Nous verrons que le critere utilise dans |98| pour caracteriser I'ecart a I'equilibre n'est pas 
satisfaisant. Dans cette etude, nous proposons de mesurer le degre d'isotropie de differentes 
observables, en consequence de quoi nos conclusions sont qualitativement differentes de 
celles de |98]. Nous etudions ensuite le scenario des minijets, recemment etudie dans la 
Ref. |99|, mais ou les auteurs n'ont pas implemente correctement les lois de conservation. 
Enfin, nous examinons la robustesse des resultats par rapport aux details du calcul, en 
particulier en ce qui concerne les approximations de couplage faible et logarithmique. 



4.3.1 "Temps de vie" 

Dans ce qui suit, nous examinons la vitesse avec laquelle differentes configurations 
initiales relaxent vers le regime hydrodynamique. Du fait de I'expansion, le systeme est 
rapidement dilue et, apres un certain temps, la description en termes des degre de liberte 
partoniques n'a plus de sens. Le point important est done de comparer le temps carac- 
teristique d'approche de I'equilibre local avec le temps de vie du systeme. Autrement dit, 
nous voulons estimer le degre d'isotropie de la distribution / a I'instant tmax ou le systeme 
de partons cesse d'exister en tant que tel. Nous donnons ci-dessous des arguments tres 
qualitatifs]^ qui n'ont de sens qu'a fortiori, les ordres de grandeur obtenus etant en accord 
avec les resultats des calculs sur reseaux (voir par exemple [|105| ). 

Une description en termes des degres de liberte partoniques n'a de sens que si la distance 
typique entre les excitations du systeme est inferieure a la taille typique d'un hadron ~ 1 fm, 
autrement dit, si la densite moyenne n d'excitations par unite de volume est superieure 
a la densite critique ric ~ 1/fm^. Dans le cas ou le systeme de gluons est a I'equilibre 
thermodynamique]^, la densite de particules est reliee a la temperature T^q par la relation 
rieg — 2Tg^g. On en deduit la temperature critique Tc ~ 1 fm~^ ~ 200 MeV, et la densite 
d'energie correspondante ec ^ 6T^ ~ 1 GeV/fm^. 

Plus la densite de particules est elevee initialement, plus il faudra de temps pour la 
diluer et done plus la duree de vie sera longue. Dans le cas present, ou le nombre total de 
particules est conserve, on pent deduire I'ordre de grandeur du temps de vie du systeme a 



^^En particulier, les estimations ci-dessous reposent sur I'hypothese que le plasma de gluons pent etre 
vu en premiere ap prox imation comme une collection de particules de masse nulle, ce qui n'est pas le cas 
(voir par exemple [|lO^ ]). 

^^Pour un ensemble de gluons a I'equilibre thermodynamique a la temperature Teq, la distribution est 
donnee par la distribution de Bose-Einstein, pour laquelle on a : Ueq — 2{N^ — 1) ("(3) Te^/Tr^, teg ~ 
6{N^ - l)C(4)r^V^^ ou C(3) ^ 1.202 et C(4) = 7r*/90 ~ 1.1. Pour iV^ = 3, on a (A^^^ - 1)/tt^ ~ 1. 
Attention, les relations precedentes sont obte nues avec une distribution de Bose-Einstein et ne doivent pas 

et ( |4.1^ ). Cependant, on note que la confusion n'aurait rien change a 
nos estimations. 



etre confondues avec les Eqs. (4.13 
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partir des seules conditions initiales : en utilisant I'Eq. (|4.5D , il vient 

tmax ^(^o) 

to ric 



(4.19) 



ou on a defini tmax 

par la relation '^(irnax^ — ^c- 

Dans chacun des deux scenarios consideres dans la suite, I'estimation de cette duree 
de vie a partir du critere e > ec donne grossierement le meme ordre de grandeur que celui 
obtenu a partir du critere ci-dessus. En revanche, le critere qui consiste a comparer le 
parametre T (Eq. (|4.8|) ) avec la temperature critique Tc n'a de sens que si le systeme est 
a I'equilibre. En effet, en general, le parametre T mesure I'energie moyenne par particule 
(T oc e/n) et n'a la signification physique d'une temperature, au sens thermodynamique 
(hydrodynamique) du terme, que si le systeme est a I'equilibre (local). 

II est clair que ces considerations tres qualitatives ne servent qu'a donner un ordre d'idee. 
Cependant les criteres n > Uc et e > ec delimitent le domaine de coherence physique de 
notre description. 

4.3.2 Le scenario de saturation 

Nous suivons ici les auteurs de la Ref. |^] pour caracteriser I'etat initial des gluons 
produits dans le scenario de saturation. La distribution initiale idealisee ( [4.20| ) ci-dessous a 
ete proposee par A. H. Mueller |^, ^ : les gluons de la fonction d'onde du noyau incident 
dont I'impulsion transverse est ~ Qs sont simplement liberes dans la collision. En resolvant 
numeriquement les equations de champ classique pour I'interaction entre les champs de 
Weizsacker- Williams non-abelien des noyau incidents, A. Krasnitz et R. Venugopalan |^ 
calculent le temps de production des gluons reels. lis obtiennent ti « 0.3 fm a RHIC 
{Qs = 1 GeV) et U « 0.13 fm a LHC (Q, = 2 - 3 GeV). 

On voit que le temps de formation est typiquement ~ 1/Qs- Get ordre de grandeur 
correspond aussi au temps necessaire pour que les gluons de differentes impulsions longi- 
tudinales (plus precisement dans differentes unites de rapidite) se separent physiquement 



les uns des autres fS^, ^3]- De meme, au bout d'un temps ~ 1/Qs, les gluons de grande 
impulsion longitudinale ont quitte la region centrale, les gluons restant ont une impulsion 
longitudinale Pz ^ Qs- De plus, du fait du phenomene de saturation, la densite de gluons 
liberes par unite d'espace des phases est ~ 1/asNc. Le nombre d'occupation dans la region 
centrale a I'instant ti peut s'ecrire 

fsatip, ti) = ^-r 5{pzti) Q{Ql - pI) , 

asNc 



ou c est une constante estimee ~ 1 par Mueller |88|, et calculee numeriquement pour 
un groupe de jauge SU{2) dans la Ref. |@ : c = 1.3. La fonction delta modelise une 
distribution tres etroite dans la direction longitudinale : Pz/Qs ~ Pzti <S 1. On choisi 



I'instant initial to comme dans la Ref. |g8[ : on laisse le systeme evoluer librement jusqu'a 



ce que le nombre d'occupation, initialement ~ l/agNc, devienne ~ 1 : 

c 

*0 = ti ■ 

L'instant to est le temps a partir duquel I'approximation de particules classiques commence 



a etre raisonnable, c'est I'instant initial pour notre etude. On a alors ||98| 

Mp) = fsatip, to) = 6{pzti) Q{Ql - pI) , (4.20) 
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Ici nous prendrons as = 0.3, Nc = 3. Les valeurs des differents parametres correspondant 



aux energies de RHIC et LHC sont resumees dans le tableau ^ 



On peut estimer le temps au bout duquel la densite de particules par unite de volume 
devient ~ nc ~ 1 fm~^ directement a partir des conditions initiales a I'aide de I'Eq. ( 4.19 ). 
On obtient tmax — 7 fm a RHIC et tmax — 30 fm a LHC. Cette derniere valeur depasse la 
limite au dela de laquelle Thypothese d'expansion longitudinale n'est plus raisonnable]^. 
En suivant revolution du systeme au cours du temps (voir ci-dessous), on voit que le 
critere e > ec donne des estimations du meme ordre de grandeur (~ 4 fm a RHIC, et 
~ 15 fm a LHC). Dans la suite, nous prendrons ces valeurs comme limites superieures (voir 
le tableau O). 



SATURATION 


Qs (GeV) 


to (fm) 


n(fo) (fm-^) 


e(to) (GeV/fm^) 


tmax (fm) 


RHIC 


1. 


0.4 


18.1 


12.0 


~ 10 


LHC 


2. 


0.18 


163.4 


217.9 


~ 30 



Tab. 4.1: Valeurs des parametres correspondant aux energie de RHIC et LHC dans le 
scenario de saturation (voir la Ref. ||9^). Le temps tmax, est la limite au dela de laquelle 
notre description n'a certainement plus aucun sens. 



Pour aller plus loin, nous devons preciser la forme du logarithme L (Eq. (^)). Dans 
les premiers instant apres t^, les particules de la region centrale ont des impulsions es- 
sentiellement transverses et interagissent entre elles en echangeant des gluons d'impulsion 
transverse. Dans ce cas, la masse de Debye (cf. Eq. ( [4. 3D ) se reduit a la masse d'ecran 
transverse [||] (voir I'Annexe |^) : 



ml'{t) 



asNc f d?p 



P 



f{p,t) 



4TTasNc 



N-i{t), 



(4.21) 



et la quantite p dans I'Eq. (|4.3|) est I'echelle d'impulsion transverse typique. Dans la suite, 
nous prendrons 



L = In 



(pi) 



nin 



(4.22) 



ou (...) = {■■■)/n designe la moyenne par particule. Bien que ce choix soit motive par 
la forme de la distribution initiale ([l.20D , hautement anisotrope, on peut voir qu'il reste 
raisonnable dans le cas d'une distribution isotrope, oii il est plus approprie de choisir 
Lisa = ln{{p'^) / rnjj) . En effet, pour une distribution isotrope m|, = 2m^ (voir Annexe 
{pf) = 2 (p2)/3, et on a L - Lisa = In 4/3 ~ 0.3. 

II est important de remarquer que la validite de I'approximation logarithmique est 
rendue tres marginale par le fait que la valeur as = 0.3 est a la limite du regime de 
couplage faible. En effet, de fagon generale, si p est I'echelle typique d'energie du probleme 
(ici p ~ Qs), on aura m'^ ~ asp'^, et L ~ ln(c/a5'), ou c ~ 1. Ici, as n'est pas reellement 
petite devant I'unite (en particulier asNc ~ 1) et L ^ l|^. C'est la valeur de la constante c 
qui controle la valeur de L. Autrement dit, la valeur de L depend des details du probleme. 



Cette hypothese est certainement valable pour t < R, on R est le rayon des noyaux incidents (pour 
un noyau de Plomb {A = 208), on estime R ~ 1.2 A^^^ fm ~ 6 fm). Pour notre propos, a caractere 
essentiellement qualitatif, la limite superieure f < 15 — 20 fm est optimiste mais pas deraisonnable. 

^^Certains choix pour L peuvent meme etre tels que L < a I'instant initial. C'est le cas par exemple 



de I'Eq. (27) de la Ref. 



pour les valeurs des parametres correspondant aux energies de LHC. 
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ce qui traduit le fait qu'a strictement parler, I'approximation logarithmique ne s'applique 
pas. Avec notre choix ( 4.22| ), il est facile de voir qu'a I'instant initial, pour la distribution 

(S) 

-^(^o) =ln( TT^Qsti 



ou ti est le temps de formation des partons de I'etat initial (Qsti ~ 1), le facteur n/iasNc 
0.8 < 1 est donne par la geometrie de la distribution initiale dans I'espace des impulsions. 
On voit que la valeur precise de Qsti est cruciale pour que L > 0. 

Pour la distribution ( 4.20| ), L ^ 1 a I'instant initial, quelque soit le choix que Ton 
fait pour I'echelle p et la masse d'ecran. Bien que la situation s'ameliore par la suite, 
comme le montre la Fig. on reste constamment a la limite du domaine de validite de 
I'approximation logarithmique. 




LHC 
RHIC 



5 10 

temps (fm) 



15 



Fig. 4.5: Le "grand" logarithme ( 4.22 ) en fonction du temps (en fm), pour les valeurs des 
parametres correspondant aux energies de RHIC et LHC dans le scenario de saturation. 

Pour tester notre approche, nous comparons nos resultats avec la solution exacte de la 



Ref. ||98|. Le premier graphe de la Fig. ^^represente revolution des quantites e, 3Pl et 3Pt 
pour les valeurs des parametres correspondant aux energies de RHIC (voir tableau |0|), les 
deux suivant representen t les moyennes par particule des carres des impulsions longitudinale 
et transverse, {p1){t) et {p\){t), pour RHIC et LHC respectivement. Ces courbes sont les 
analogues de celles des Figs. 7 et 9 de la Ref. fo^ , on a choisi les echelles et les unites de 
maniere a faciliter la comparaison : non seulement I'approche a I'equilibre est correctement 
reproduite, mais on voit que dans I'ensemble, nos resultats sont en accord semi-quantitatif 
avec la solution exacte, ce qui est remarquable compte tenu du degre de simplicite de notre 
approche. 

Venons-en maintenant aux resultats de notre etude du scenario de saturation. Pour 
caracteriser I'ecart a I'equilibre cinetique local, c'est a dire le degre d'anisotropie de la 
distribution, nous mesurons les rapports 



Rk 



(pi// 



(4.23) 



en fonction du temps. Pour une distribution completement isotrope, K"^" 



moments du type {pz. 



1 , V/c. Les 

Ip^) sont sensibles a des echelles d'impulsion de plus en plus faible 
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Fig. 4.6: Reproduction de differentes courbes de la Ref. ||9^. Les differentes quantites presentees ici 
sont obtenues dans rapproximation du temps de relaxation. Les echelles et les unites ont ete choisies 
de maniere a faciliter la comparaison. En se reportant a ||9^ on trouve un accord semi-quantitatif. 
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a mesure que k augmente, les rapports Rk (plus precisement I'ecart 1 — Rk) mesurent done 
le degre d'anisotropie des differentes parties de la distribution, c'est a dire des differentes 



echelles d'impulsions. La Fig. 4/7 represente revolution temporelle des rapports Rk{t) pour 
/c = 0, 3, a RHIC et LHC. On observe le comportement qualitatif attendu : les moments 
sensibles aux echelles d'impulsions les plus grandes sont ceux qui relaxent le moins vite 
vers I'equilibre. Cette hierarchie de vitesses de relaxation, deja observee dans [^], est un 
simple effet de la geometric du systeme0, laquelle est prise en compte exactement dans 
I'Eq. ( CT ). 
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Fig. 4.7: Evolution temporelle des rapports Rk definis par I'Eq. ( |I2^ pour le scenario de satu- 
ration a RHIC (a gauche) et LHC (a droite). Les ecarts 1 — Rk mesurent le degre d'anisotropie de 
la distribution microscopique f{p,t), pour differentes echelles d'impulsions. 



Le processus de relaxation est bien plus rapide a LHC et, la duree de vie y etant plus 
grande, le degre d'isotropie atteint est bien meilleur. On voit clairement que pour t ~ t^ax, 
le systeme est encore tres anisotrope a RHIC, tandis qu'a LHC le degre d'isotropie est 
tres bon. Dans ce dernier cas, le degre d'isotropie est deja satisfaisant pour t > 15 fm oii 
^fc=o,...,3 > 0.8. 

Ces resultats semblent en complete contradiction avec ceux de la Ref. ou les auteurs 
mesurent les "temps d'equilibration" teq = 3.24 fm et teq = 2.36 fm pour les valeurs des 
parametres correspondant a RHIC {as = 0.3, Qs = 1 GeV) et LHC {as = 0.3, Qs = 2 GeV) 
respectivement. Cette contradiction apparente est due a ce que les criteres utilises pour 
caracteriser I'ecart a I'equilibre ne sont pas equivalents : dans |Q, teq est defini comme le 
temps auquel les quantites Tt^^^, et s/n, I'entropie moyenne par particule, atteignent 90% 
de leurs valeurs respectives a I'equilibre (ces deux quantites sont constantes dans le regime 
hydrodynamique). Par ailleurs, on pent voir sur les Figs. 7 et 9 de |Q (voir aussi notre 

{pI)/{pI 



Fig. |4.6|) , que les valeurs des rapports Ri = Pl/ Pt et Rq = (pj,) / {p\) = {pi 
du temps sont en accord avec celles de notre Fig. [4.7|. En particulier, pour t 



"eq 



au cours 
: 3.24 fm, 

on & Pl/Pt — 1/2. Le critere que nous avons utilise pour definir I'equilibre cinetique, 
c'est a dire I'isotropie^, est d'une part plus intuitif et d'autre part plus precis. En effet les 



^^Les particules de grande impulsion longitudinale desertent rapidement la region centrale du fait de la 
seule expansion. II s'ensuit que les particules "dures" (de grande energie) de cette region ont une impulsion 
essentiellement transverse. II faut done plusieurs collisions a petite deviation pour creer une particule de 
grande impulsion longitudinale, cela prend beaucoup de temps, plus qu'il n'en faut a cette derniere pour 
quitter la region centrale. 



II est facile de montrer que I'Eq. (4.1) ne pent avoir de solution isotrope que dans le cas oil le membre 
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valeurs asymptotiques des rapports Rk sont bien connues : Rj^ ^° = I, tandis que celles des 
quantites utilisees dans [98| dependent de I'histoire du systeme. De plus le degre d'isotropie 
du systeme est I'observable importante du point de vue phenomenologique. 



4.3.3 Le scenario des minijets 

Nous suivons ici le choix des auteurs de la Ref. |Q pour caracteriser la distribution ini- 
tiale dans le scenario des minijets. Dans la region centrale, celle-ci est prise de la forme ( |4.8| ) : 

Mp) = fjetip, to = l/po) = Xjet e-P/^^^' , (4.24) 



A(fo) et Tjet = T(to) sont determines a partir les densites de 
'eq{to) et d'energie ejet = e(to) = ^eqito) du systeme de partons 



ou les parametres Xjet 
particules njet = n{to) ■ 
initialement produit, a I'aide des Eqs. (|4.13| ) et ( 4.14 ). La contribution dominante vient des 
partons d'energie po, dont le temps de formation est to = 1/Po- Les valeurs des parametres 
de ce scenario pour des collisions frontales Au-Au a RHIC et Pb-Pb a LHC, sont donnees 
dans le tableau 4.2. Dans les deux cas, les limites superieures tmax pour les temps auxquels 



n ~ nc et e ~ ec sont du meme ordre de grandeur que celles obtenues dans le scenario de 
saturation. 



MINUET 


to (fm) 


rijet (fm"^) 


ejet (GeV/fm^) 


Xjet 


Tjet (GeV) 


tmax (fm) 


RHIC 


0.18 


34.3 


56.0 


1.0 


0.535 


~ 10 


LHC 


0.09 


321.6 


1110.0 


1.0 


1.13 


~ 30 



Tab. 4.2: Valeurs des parametres correspondant aux energie de RHIC et LHC dans le 
scenario des minijets (voir la Ref. p9| ). 



De meme que dans I'etude du scenario de saturation, nous prendrons as = 0.3 et ferons 



le choix ( 4.22 ) pour le logarithme L. En ce qui concerne la validite de I'approximation 
logarithmique, bien que toujours a la limite, la situation est plus confortable ici que dans 
le cas du scenario de saturation : dans I'intervalle de temps etudie on a 1.5 < L < 5. En 
particulier, il est facile de voir qu'a I'instant initial 



L{to 



In 



7tt 



1 



AasNc Xjet 



Encore une fois le profil de la distribution joue un role important (c'est lui qui determine 
le facteur Jn/^asNc). 

On voit sur la Fig. ^]8| que I'isotropie initiale est rapidement detruite dans les premiers 
instants de revolution. C'est la phase initiale de regime libre dont nous avons parle plus 
haut : a I'instant initial, le temps de relaxation est non-nu|^ et il faut attendre un temps 
At = t — to < 9{tQ) ~ 1 fm) avant que I'effet des collisions ne se fasse sentir, I'isotropie est 
alors lentement reconstruite. 

Apres que I'isotropie initiale ait ete detruite par I'expansion, la distribution est tres 
piquee autour de pz = 0, et revolution ulterieure est tres semblable a celle du scenario 
de saturation. On voit cependant que le processus de relaxation est plus lent dans le cas 



de gauche est negligeable, ce qui correspond, par definition, au regime hydrodynamique (cf. Eq. (1.7)) 
^^Mentionnons a ce propos un detail technique concernan t le c alcul du tem ps de relaxation a I'instant 



initial dans le scenario des minijets : pour la distribution ( |4.24| ), I'Eq. ( 4.11 ) est triviale (0 = 0) et ne 
permet done pas de calculer le temps de relaxation 9{to). On contourne la difficulte en exploitant le fait 
que le systeme est initalement dans un regime libre (voir l'Annexe|^). 
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Fig. 4.8: Evolution temporelle des rapports Rk (Eq. ( 4.23 ) pour le scenario des minijets a RHIC 
et LHC. 



present et les degres d'isotropie atteints sent par consequent nettement inferieurs, aussi 
bien a RHIC qu'a LHC (les temps de vie sent du meme ordre de grandeur dans les deux 
scenarios (voir Tabs. ([4.1D et ([4.2] ) ). Dans ce dernier cas, bien que le degre d'isotropie finisse 
par etre relativement bon, il faut attendre t ~ 30 fm pour avoir -Rfc=o,...,3 ^ 0.8, soit deux 
fois plus longtemps que dans le scenario de saturation. Sur la base de ce critere, les temps 
d'equilibration obtenus dans le scenario des minijets sont » tmax a RHIC, ~ tmax a LHC. 

Ces conclusions sont en disaccord qualitatif avec celles des auteurs de la Ref. qui 
mesurent des temps d'equilibration similaires ~ 4 — 5 fm a RHIC et LHC. Cependant 
I'approche utilisee dans cet article n'est pas satisfaisante, en particulier a cause du fait que 
bien que seuls les processus elastiques y soient consideres, le nombre total de particules n'est 
pas conserve. En fait il augmente avec le temps (la densite moyenne par unite de volume 
n decroit moins vite que ce qui diminue les temps d'equilibration (les collisions etant 
plus frequentes). Le fait que ceux-ci soient du meme ordre de grandeur a RHIC et a LHC 
vient de ce que dans |99|, la constante de couplage n'est pas constante, mais augmente 
avec le temps, sa valeur etant calculee a chaque instant en fonction de I'echelle typique 
d'impulsion des particules du milieu. Celle-ci etant plus petite a RHIC qu'a LHC, I'intensite 
de I'interaction, et done la frequence des collisions y est augmentee. Dans la partie suivante, 
ou nous considerons le cas d'une constante de couplage variable, nous verrons cependant 
que cet effet ne modifie pas nos conclusions presentes. 



4.3.4 Robustesse des resultats 

Pour obtenir une information la plus liable possible dans I'approche utilisee ici, il est 
important de tester la sensibilite des resultats precedent par rapport aux details de la 
description, comme par exemple la valeur de as, ou encore le choix de I'echelle d'impulsion 
p sous le logarithme L (cf. Eqs. ( f4.3D ). En effet, etant a la limite du domaine de validite de 
I'approximation logarithmique, on s'attend a ce que de petits changement sous le logarithme 
aient des effets non-negligeables. 



Constante de couplage variable 



II est interessant de considerer le cas ou la valeur de la constante de couplage as 
augmente avec le temps, a mesure que decroit I'energie typique des particules du milieu. 



4.3. Conditions initiales et resultats 



113 



Comme nous I'avons mentionne plus haut, on s'attend alors a ce que les particules inter- 
agissent plus fortement a RHIC qu'a LHC, I'energie moyenne par particule y etant plus 
faible. 

On calcule la constante de couplage a chaque instant en fonction de I'energie moyenne 
par particule : as{t) = asi/J- = e = 3T), avec la parametrisation 

as{,)=asiM,) , (4-25) 

ou Mz — 90 GeV est la masse du boson Z, as{Mz) = 0.1, et ou on a pris Aqcd — 200 MeV. 
Dans le scenario des minijets, a LHC, on a 0.2 < as{t) < 0.4 sur I'intervalle de temps 



considere et on obtient un resultat tres similaire a celui de la partie ^.3.3| , obtenu avec 
as = 0.3. Bien que, par rapport au cas ou as est constante, la situation s'ameliore a RHIC, 
ou 0.3 < as{t) < 0.5, I'equilibre n'est pas atteint. Les courbes representant revolution de 
i?i sont presentees dans la partie suivante. 

Dans le scenario de saturation, I'energie moyenne par particule est plus faible et la 
constante de couplage devient rapidement > 0.5. Les resultats obtenus n'ont pas de sens 
physique. Ceci doit etre pris comme une manifestation de la fragilite du calcul perturbatif. 

Estimation des incertitudes 

Pour nous faire une idee plus generale de I'ordre de grandeur des incertitudes dues 
aux details de notre approche, nous remplagons simplement L ^ 2L, puis L L/2, 
dans le terme de collision^ (on travaille ici avec as = cte = 0.3). Les courbes des parties 
precedentes deviennent des "bandes" dont la largeur represente le degre de confiance avec 
lequel interpreter nos resultats. Les bandes correspondant a I'observable Ri = Pl/Pt (cf. 
Eq. ([4.23D) pour les scenarios de saturation et des minijets sont representees sur les Fig. 



et [4.10| respectivement. Sur la Fig. fl.lOj , on a aussi represente les courbes correspondant 
au cas ou la constante de couplage depend du temps (voir la partie precedente), on voit 
que celles-ci sont a I'interieur du 'domaine d'incertitude" decrit ci-dessus. Les observations 
qualitatives des parties precedentes restent valables. Cependant, a cause des incertitudes 
liees a notre description, on ne pent pas dire si I'equilibre est atteint a LHC dans le scenario 
des minijets. Dans le scenario de saturation, ou I'equilibre est certainement atteint, il est 
est impossible d'estimer precisemment la duree du regime transitoire. 

4.3.5 Conclusions et perspectives 

En resume, nous avons etudie revolution du systeme de gluons, produit dans les tout 
premiers instant d'une collision nucleaire a tres haute energie, et en particulier la contribu- 
tion des collisions elastiques de petite deviation au processus d'equilibration cinetique. Nous 
avons etudie et compare deux scenarios differents pour I'etat initial du systeme, en resolvant 
les equations cinetiques dans une approximation de temps de relaxation "auto-coherente". 



laquelle donne des resultats en bon accord avec la solution exacte obtenue dans |98| pour 



le scenario de saturation. Nos conclusions sont resumees dans le tableau [4.3| . Dans les deux 
scenarios etudies, les processus elastiques ne sont pas sufRsamment efficaces pour rendre 
la distribution isotrope avant la phase hadronique aux energies de RHIC {teq ^ tmax)- 
Dans le scenario de saturation, pour LHC, le regime hydrodynamique est atteint. La deter- 
mination exacte du temps d'equilibration n'est cependant pas possible etant donnees les 



^^Le changement L ^ aL revient a remplacer le temps de relaxation 6 par 9 /a dans I'Eq. ( [l.l6| ). Quand 
= 2 (a = 1/2), le systeme relaxe plus rapidement (lentement). 
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Fig. 4.9: Dependance des resultats par rapport aux details de la description pour le scenario 
de saturation. Les courbes delimitent la marge de variation de revolution temporelle du rapport 
i?i — Pl/Pt quand on change L 2L (pointilles) et L L/2 (tirets). 
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Fig. 4.10: Dependance des resultats par rapport aux details de la description pour le scenario 
des minijets. Les courbes representent la marge de variation de revolution temporelle du rapport 
i?i — Pl/Pt quand on change L ^ 2L (pointilles) et L L/2 (tirets). Les courbes en trait 
continu sont obtenues avec une constante de couplage variable, dependant de I'energie moyenne 
e = 3T des particules du milieu (voir Eq. ( 4.25| )). 
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incertitudes de la description utilisee. Dans le scenario minijet, il apparait meme difficile 
de dire si I'equilibre local est atteint. 





SATURATION 


MINIJET 




RHIC 


LHC 


RHIC 


LHC 


tmax if'^^ 


~ 10 


~ 30 


~ 10 


~ 30 


Ri 


<0.8 


>0.8 


<0.8 


>0.6 


Regime liydrodynamique ? 


non 


oui 


non 


? 



Tab. 4.3: Le systeme a-t-il atteint le regime hydrodynamique ? 

Les auteurs de la Ref. [|106| ont recemment argue que les processus inelastiques de 
branchement (5+milieu gf^f+milieu), contribuant au meme meme ordre que les processus 
elastiques dans I'approximation logarithmique, jouent un role important dans le processus 
d'equilibration cinetique, et accelerent considerablement celui-ci. II est fort probable que 
ces processus inelastiques permettent de rendre le temps d'equilibration negligeable a LHC 
(quel que soit le scenario le plus adapte). La situation est moins claire a RHIC. II est 
tres important, pour I'interpretation des donnees experiment ales actuellement accumulees 
a RHIC, de savoir si le regime transitoire pent etre neglige, ou si I'aspect hors equilibre 
du systeme doit etre pris en compte dans les calculs concernant les signatures directes du 
plasma. 
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Annexe F 

L'integrale de collision de Landau 



Dans le Chap. § nous nous interessons a la matiere produite immediatement apres 
I'impact, dans la region centrale d'une collision d'ions lourds ultrarelativistes. Nous sup- 
posons que la densite d'energie est suffisament elevee pour que cette matiere puisse etre 
decrite par un ensemble de partons (essentiellement des gluons) interagissant faiblement 
les uns avec les autres. Autrement dit, nous supposons que la constante de couplage forte 
est tres petite : 05 <C 1. Nous supposons de plus qu'apres quelques instants, le systeme est 
suffisament dilue pour pouvoir etre decrit par une equation de Boltzmann pour la densite 
de gluons par unite d'espace des phases, que nous supposons independante de I'etat de 
polarisation et de couleur des gluonsQ 

fiP,x,t)- ^ ' 



Construisant le terme de collision de I'equation de Boltzmann correspondante, nous nous 
limiterons a la contribution des collisions elastiques gg gg, a I'ordre dominant en as- 
Nous devons considerer tous les processus du type (e = +, — designe la polarisation et 
a = 1, — 1 la couleur) : 

ip,e,a) + {p',e,a') (pi,ei,ai) + {pi',e[,a'i) , 

ainsi que le processus inverse. Dans les deux cas il nous faut sommer sur tous les etats 
d'impulsion, de polarisation et de couleur des gluons autres que celui auquel on s'interesse 
(le premier ci-dessus), et moyenner sur les etats de polarisation et de couleur de ce dernier, 
schematiquement : 



2{N^ - 1) 



t',a' ei,ai e[,a[ 



II est commode de retablir une certaine symetrie entre les deux gluons {p, e, a) et {p' , e' , a') 
en exprimant la probabilite de transition apparaissant sous I'integrale de collision en termes 
du carre de I'amplitude de diffusion \ Aigg^gg\'^ , somme sur les etats de polarisation et de 



strictement parler, la construction du terme de collision qui suit n'est valable que dans le cas 
d'une distribution independante de la couleur et de la polarisation. Cependant, on pent voir I'equation 
de Boltzmann obtenue comme decrivant revolution du systeme sur des echelles de temps grandes devant 
Techelle caracterisant la relaxation de couleur. En effet, cette derniere est beaucoup plus rapide que I'echelle 
de relaxation vers I'equilibre cinetique a laquelle on s'interesse ici. La distribution / doit etre vue comme 
une moyenne sur les etats de polarisation et de couleur. 
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couleur de I'etat final et moyenne sur ceux des deux gluons de I'etat initial. Ceci revient a 



multiplier et diviser par 2{N^ 



2{Nl 



11 



2(iV2 - 1) ^ 2(iV2 



1) ^ 



EE 



L'integrale de collision s'ecrit alors 



c[fm 



2(iV2 



11 



d^p' d^pi d^p[ 



w{pp' 



Pi Pi ') 



(27r)3 (27r)3 (27r)3 

[/i/{(l + /)(l + /')-//'(l + /i)(l + /0] 



(F.ll 



ou / = f{p), f = f{p'), fi = fipi) et /( = f{pi') (toutes les fonctions aparaissant 
dans le terme de collision sont prises au meme point d'espace-temps, nous avons done omis 
d'indiquer la dependance en x et t) et oil 



w{pp' pi pi ') = w{pi pi ' pp') = (27r)^(5(^) {P + P' -Pi-P[) 



\M 



16pp' Pip[ 



Le fait que les gluons sont des bosons indiscernables est pris en compte dans le calcul 
de I'amplitude de transition \Aigg^gg\'^ (terme d'echange), et resulte dans le fait que la 
probabilite de transition w est symetrique dans I'echange pi <-> pi ' (voir la symetrie sous 
I'echange u ^ t dans I'expression ( |F.3D ). Apres avoir integre sur toutes les valeurs possibles 
de pi et pi on a done compte deux fois I'ensemble des configurations possibles, d'ou le 
facteur 1/2 devant l'integrale. 

Dans le cas qui nous interesse ici, il est possible de calculer la contribution dominante 
a cette integrale en exploitant le fait que I'amplitude de diffusion diverge fortement dans 



la limite des diffusions a petit angle de deviation |102|. 



La section efficace de diffusion gg gg 

La section efficace differentielle de diffusion s'ecrit 

. .<^.,P + P' - P, - P;, ^1=^ ^ ^ , (F.2, 

ou \Aigg^gg\'^ cst Ic carrB de I'amplitude de diffusion gg — > gg somme sur les etats de 
polarisation et de couleur des gluons de I'etat final, et moyenne sur ceux de I'etat initial. 
En termes des variables de Mandelstam 

on a (voir par exemple ||88|| 1 

, . , ,9 A f ut US st\ , , 

avec g'^ = Anas- Dans le referentiel du centre de masse de la collision gg gg '■ 

Pf" = p (1,0, 0,1) , Pf =p(l,0,sin6',cos6l) , 
(P'r = p(l,0,0,-l) , {P{r =p {1,0, -sin 9, -cos 9) , 
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, t = -2p^{l-cose) , u = -2p^ (1 + cos( 



on a 



La section efRcace ( [F.2| ) diverge fortement pour les petites valeurs de t ss 0) et, par 
symetrie, pour les petites valeurs de u {6 ^ it). Ces deux regions de I'espace des phases 
dominent done I'integrale de collision, et leurs contributions respectives sont identiques 
par symetrie. II sufRt done d'en calculer une et de multiplier le resultat par un facteur 2, 
qui vient se simplifier avec le facteur de symetrie 1/2 devant I'integrale (Eq. ( [F.l| )). Tout 
se passe alors comme si les gluons etaient des particules discernables. Dans la suite nous 
omettrons des le depart les facteurs de symetrie et calculerons la contribution de la region 
« (diffusion a petit angle de deviation), oii Ton a 

.to^W,-^^^, (F.4) 

L'integrale de collision dans la limite des petits angles 

Du fait de la conservation de I'impulsion, I'integrale sur pi ' est triviale, et la probabilite 
de transition w ne depend que de trois vecteurs. II est utile de I'exprimer comme une 
fonction des demi-sommes des impulsions initiales et finales des deux particules : {p+pi)/2 
et (p' +pi')/2. Remplagant I'integrale sur pi par une integrale sur I'impulsion ecliangee 
q = Pi — p, on ecrira 

w{pp' pipi ') = w{p+ - ~ + ^ {diW - d[w) , 

ou u; = w{p,p';q 0), et di = d/dp^, le prime designant la derivation par rapport aux 
composantes de p' . On developpe de meme les poids statistiques : 

q^q^ 



h f + q'dif + ^d,,f 



2 



OU dij = d^/dp^dp' . Dans le developpement de la combinaison des poids statistiques ap- 
paraissant sous I'integrale de collision, le premier terme non-nul est {F'dif — Fd[f')w, 
ou F = /(I + /). La probabilite de transisiton w{pp' PiPi') etant symetrique dans 
I'echange des etats initial et final (autrement dit, c'est une fonction paire de q), ce terme 
ne survit pas a I'integration sur d^q. La premiere contribution non-nulle a I'integrale de 
collision vient du terme du second ordre dans le developpement de I'integrand 

C ^ / [{d,B'^) - (dlB^^)] [F' {d,f) - F {d'^f')] 
Jp' 

+ / {F' {d,,f) + F (dy) - [1 + / + /'] [(dj) (d'^f) + {d,f) {d'J')] } B^^ , 

Jn' 



ou 



et ou Ton a note 



/ 



d^p' 
(2^ 
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En utilisant la relation B^^ = B^^, on reecrit le dernier terme sous la seconde integrale : 

[1 + / + /'] (d'^f) + (5,/) id'j')] B'^ = B'^ [{d,F) {d'^n + id,f) idlF')] , 

et, en reorganisant les differents termes, on obtient (les derivees agissent sur tout les termes 
entre accolades) 

C^- [jd,- d[) {B^^ [F {d'^f) - F' (a,/)] } . 

Le second terme est l'integrale d'une divergence totale et est done nul. On obtient alors 
(voir ill) 

C ~ -di , (F.6) 



ou 



s^= f B'^ [F{d'jf')-F'{d,f)] . (F.7) 

II n'est pas surprenant que le terme de collision s'exprime comme une divergence dans 
la limite des petits angle : dans I'espace des impulsions, la particule est reperee par le 
vecteur p. Apres une collision typique, elle se retrouve a la "position" voisine reperee par 
le vecteur p + 6p. On voit done que les particules se deplacent le long de "trajectoires" 
continues dans I'espace des impulsions. II s'agit d'un processus diffusif dans cet espace et, 
par consequent, le terme de collision correspondant doit s'exprimer comme la divergence 
du flux de particule dans I'espace des impulsions. On pent construire explicitement ce 
flux ||102|] , le resultat obtenu comcideQ avec I'Eq. (F.7). 



II nous reste a calculer les composantes du tenseur B'^^ . Par definition de la section 
efficace (Eq. ( |F.2| )), on a 

W{pp -^pipi) = Vrel da , 

avec 

p.p' 



Vrel = 1- = \/ {v -V'Y - {v X v'Y 



On a done (g, = —q^) 

B'^ = Bij = ^ y qiQj Vrel da . 



En suivant la discussion de la Ref. [[102| , il est facile de deduire la structure generale de ce 



tenseur. Pour des particules de masse nulle, on obtient, en utilisant I'Eq. 

= L[{l-v- v')6ij + Viv'^ + vjv';\ , (F.8) 

ou 

, fde^ 

±J - 



Cette integrale sur les angles diverge logarithmiquement. La divergence correspondant aux 
petites valeurs de a pour origine physique la portee infinie de I'interaction entre les deux 
partons (gluons), dans la limite de couplage faible (echange d'un gluon de masse nulle). 



2^ 



Ceci est fait explicitement dans la Ref. |102| pour le cas d'un gaz de particules classiques. II est facile 



de repeter la discussion en prenant en compte la statistique quantique. 
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En realite la portee de cette interaction est ecrantee par la presence des autres partons 
du milieu, c'est le phenomene d'ecrantage de Debye : le gluon echange, de type electrique, 
a une masse effective (il est habille par les interactions avec le milieu) m^), la masse de 
Debye. L'integrale doit done etre limitee a la valeur minimale Omin, correspondant a la 
valeur minimale de la quadri-impulsion echangee, egale a la masse de Debye : 

-tmin = m\,r^j^ el^i^ , (F.9) 

ou p designe I'ordre de grandeur de I'impulsion typique des partons du milieu. La divergence 
de l'integrale L pour les grandes valeurs de 9 est non-physique et reflete le fait que les 
formules precedentes ont ete obtenues dans I'approximation des petits angles, qui cesse 
d'etre valable quand ^ ~ 1. On a done, a I'approximation logarithmique (les termes negliges 
dans le developpement doivent etre petit devant le logarithme L), 

V ^min / 

Calcul du terme de collision 

En developpant des formules ( [F.6| )-( ]FJ| ) il vient, apres quelques integrations par parties, 
C = {d,jf) [ B'^F' + d,F [ f'id'^B'i) 

Jp' Jp' 

+F [ f{d'^d,B'^) + (djf) [ F'im^) , 
Jp' Jp' 

avec diB''^ = 2Bv'^ /p et d'jdiB"^^ = AB/pp' (on a note B = 7ra| j^'L^ L). Dans la region 
centrale de la collision (voir le Chap. la symetrie du probleme est telle que la distribution 
f{p) est une fonction paire de p. Les integrales du type / d^p f v'^ sont done nulles et 
I'expression ci-dessus se simplifie consider ablement. En remarquant que = V -v = 2/p, 
on obtient finalement 

C[/]|,=o = SMo V2/ + 2eiV_i V [vf{l + /)] , (F.ll) 
ou Ton a defini les moments 



d?p 
(2^ 



Ms = 2(iV2-l) / /(! + /), (F.12) 



Ns = 2{N',-l) I ^ff- (F.13) 



L'Eq. (F.ll) ci-dessus se simplifie encore dans le cas d'un gaz non-degenere, c'est a dire 



d'un gaz classique (/ <C 1), ou on a 

C[/]L=o = + 2BN^i V{vf) . (F.14) 

Nous utilisons cette derniere forme du terme de collision dans I'etude presentee au Chap. ^. 
Cette approximation est motivee par le fait que le systeme est rapidement dilue a cause 
de I'expansion. Les partons peuvent alors etre consideres comme des particules classiques 
(par opposition a I'approximation de champ classique du Chap. ||, valable quand f ^ 1). 
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Annexe F. L'integrale de collision de Landau 



Masse de Debye et masse d'ecran transverse 

On peut comprendre simplement I'origine physique du phenomene d'ecrantage de Debye 
de la fagon suivante : supposons que deux particules, plongees dans un milieu, tentent 
d'interagir I'une avec I'autre en echangeant un quanta virtuel. Celui-ci est emis par la 
particule 1 et se propage "en direction" de la particule 2. II est clair que plus la distance 
a parcourir est grande, plus la probabilite pour que celui-ci soit absorbe par une particule 
du milieu avant d'arriver a destination augmente : la portee de I'interaction, infinie dans le 
vide, est limitee du fait de la presence du milieu. Ceci se traduit par le fait que le quanta 
echange correspond a un champ effectivement massif, "habille" par ses interactions avec les 
particules du milieu. La masse correspondant a I'ecrantage de la partie electrique du champ 
correspondant (la seule qui joue un role dans notre cas) est appellee masse de Debye. Dans 
le cas d'un systeme de gluon, les auteurs de |107| obtiennent rexpression]^ : 



ml = lim / d^pil^^Ml) , (F.15) 

ouv = p/p et q = \ \q\\, q etant I'impulsion du gluon echange. On voit que, de fagon generale, 
la masse d'ecran depend de la direction de cette derniere. Dans le cas d'un gluon d'impulsion 
transverse q= {qt,qz = 0) et pour une distribution independante de la direction dans la 
direction transverse f{p,t) = f{pt-,Pz-,t), on obtient, apres integration par parties, 

Dans le Chap. nous aurons a considerer le cas oil la largeur de la distribution dans la 
direction longitudinale est negligeable : f{p,t) oc S{pz)- On obtient alors I'expression de la 
masse d'ecran transverse 

2 asNc f d?p 
TT^ J p 

Nous aurons aussi besoin du cas ou la distribution / est isotrope. On a alors Vp/(p, t) = 



V f'{p,t), oil /' = df /dp, et la masse d'ecran ( F.15 ) est independante de la direction du 
gluon echange, comme il se doit (il est facile de voir que I'expression ( F.16| ), obtenue pour 
une distribution generale, se reduit a la formule ci-dessous dans le cas d'un distribution 
isotrope). On obtient, apres integration par partie. 



2 2asN^ fd^p 
"^0 = —^ 



l^fip,t). (F.18) 



Bien que I'expression ( |F.17| ) de la masse d'ecran transverse ait ete obtenue pour le cas 
d'une distribution oc 6{pz), il nous sera pratique de I'extrapoler au cas oii la distribution 
est isotrope. On a alors m'^ = ml/2. 



''La derivation de cette formule n'est pas presentee dans la reference citee. On peut obtenir une formule 
identique, au facteur de coule ur pre s. dans le cas d'un plasma QED, dans la theorie de la rfeponse lineaire. 
Nous pensons que la formule ( F.15 ) est obtenue avec le meme type d'approximations pour un systeme de 
gluons 



Annexe G 

Le calcul des moments 



Dans cette annexe, nous derivons les formules utilisees au Chap. ^ pour calculer le 
temps de relaxation dans la methode des moments. Dans le cas du terme de collision de 
Landau (voir Annexe |F|), il est possible de calculer explicitement toutes les integrales sur 
I'espace des impulsions pour tous les moments dont nous avons besoin au Ch ap. ^ . 

La solution formelle de I'equation du temps de relaxation s'ecrit (cf. Eq. ([4.16| )) 

f{p,t) = Mpt,P.-)e~^^'^+ dt'—^^f,,{pt,p,-,t'), (G.l) 

oil 

. dt' 

x{t) 



to 



9{t') 



et ou fo{p) est la distribution initiale et feq{p,t) = A(t)e 



'p/m 



De fagon generale, nous voulons calculer les integrales du type 

M{t) = {m)= [ m{p)f{p,t), 



avec = 2{N^ ~ 1) / d^p/(2vr)^. Considerons par exemple le cas M = Ng = (p^) (cf. 
Eq. ([4.2D ), et concentrons nous d'abord sur le second terme du membre de droite de 
I'Eq. (|G.1|). On doit calculer I'integrale 



i!\t,t') = J^yf,,{pt,p,^,t')). 



En faisant le changement de variable Pzt/t' — > pz et en utilisant un systeme de coordonnees 
spheriques pour p, on obtient aisement 



h,{t,t') = jN^'^{t')hs (0 



avec 

Jp 



et ou la fonction 

-1 



hs{a)= I dx[l-{l-a^)x^X''^ . (G.2) 
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Annexe G. Le calcul des moments 



Pour calculer la contribution du premier terme du membre de droite de I'Eq. ( G.l ), on 
doit specifier la condition initiale. Dans le Chap. on considere deux type de conditions 
initiales : le scenario des minijets, ou fo{p) = feg{p,t), et le scenario de saturation, oil 
foip) = ^(Pz) g{p±)- Dans la suite, nous definissons la notation 

A Si fo{p) = feq{p,to) 



{A\\B] 



B si /o(^ = S{p^)g{p2 



On obtient finalement 



te^Nsit) = toNM {hs{to/t)\\l} + I'^t'e^'N^^t') hs{t'/t) , (G.3) 

Jto ^ 

oux = x{t), x' = x{t'), e' = e{t'). 

On calcule de maniere analogue les moments "longitudinaux" = (par 
exemple, la pression longitudinale Pl = {p1/p) = ) : 

j-t If/ 

^.3 „x ATZ/+\ _ 4-3 AT thzu^ /^.\\\n^ i / W\3 „x' ]\Teq W\ hz f4-' /-4\ (G 4) 



t'e^N^{t)=4NM {K{to/t)\\0}+ / {t're^ N!^{t')K{t'/t) 



to 



avec 



a)= [ dxx^ [l - (1 - a^) 
Jo 



(G.5) 



Enfin, les moments "transverse" = {p\p'^ ^) sont aisement obtenus a partir des prece- 
dents : 

A^, = Nl + 2N^ . 

Nous donnons ci-dessous les expressions des fonctions hs et /if pour s = —2, ...,2. Les 
integrales ( G.2| ) et (|G.5D se calculent aisement. On peut aussi utiliser la relation suivante : 

s da? 



(G.6) 



En notant A = ^/\ — a?, 



hs{a) 



[ dx a^) 
Jo 



/if (a) = / dxx^ [1- (l-a2)x2] 
Jo 



2\^21*/2-l 



/12(a) 



2 + 



If arcsinA 
/11(a) = - ( a + J 



/io(a) = 1 
/i_i(a) 
h-2{a) 



arcsin A 
I 

arctan A 
A 



/if (a) 
/if (a) 

hl{a) 



1 

3 



1 / arcsin yl 



2^2 V A 

1 / arctan yl 
1^ A 



, ^ , , 1 /I arcsin A 
'-^'^) = I^ a A- 



1 /I arctan A 



2A2 Va2 



A 
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Dans le Chap. ^, nous avons resume les Eqs. ( |G.3D et (|G.4D ci-dessus sous la forme (cf. 
Eq. §13) 

M(i) = M(io)^S^(to/t)e-^'+ fdt'^J^^^''\t'/t)M,,{t'), 

ou M designe un moment quelconque. On a 

= a{hs{a)\\l} , T^^f{a) = ah,{a) 
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Annexe H 

Calcul de 6{t(y) dans le scenario des 
minijets 



La methode des moments, exposee au Chap. ^, ne permet pas de calculer le temps de 
relaxation dans le cas oil / — feqj comme c'est le cas dans le scenario des minijets, a I'instant 
inital. Pour calculer 0(to); il sufRt cependant de s'ecarter de feg de fagon infinitesimale en 
expoitant le fait que dans les tout premiers instants, le systeme ne ressent pas encore I'effet 
des collisions et est dans un regime quasi-libre. Nous detaillons ce calcul dans cette annexe. 

La distribution initiale est donnee par (Aq = Xjet et Tq = Tj^t) 

/(p,to) = Aoe-P/^». 

Le temps de relaxation etant non-nul a I'instant initial, il est clair que pour les temps t 
tels que t — tQ^ S{to), le systeme est dans un regime essentiellement libre. L'evolution est 
alors donnee par I'Eq. (|4.6D, et on a 



V Tjet P 

On en deduit, pour les moments Ns, et N^, definis dans I'annexe |G|, 

Ns{to + 6t) = Nsito) - Ns+iito) , 

61 



et 



N^ito + 5t) - N^ito + 5t) = ——- Ns+i{to) , 

15To 



oil Ton a utilise le fait que la distribution est isotrope a to. Avec ces notations, I'Eq. ( |4.18D 
s'ecrit : 

ATZ _ AT± / \ / 



oil 



Ar2 



et L est donne par les Eqs. ( [4.2lD -( ^^ ). En se plagant a I'instant t = to(l + et 
en utilisant les expressions ci-dessus, il vient, au premier ordre non-trivial en e (on ecrit 

e{to{i + e)) = e{to) + ee{to)), 

^ = AB{Nof + 2BN^iNu 
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Annexe H. Calcul de 9{to) dans le scenario des minijets 



ou toutes les quantites sont calculees a to [Oq = 0(to))- En utilisant 



Nsito) = is + 2)l^-^XoT^ 



s+3 



On obtient finalement 



Avec les valeurs des parametres as = 0.3, A'^c = 3 et Aq = 1, correspondant au scenario des 
minijets (voir le Chap. on obtient 

Jo 

c'est a dire 6*0 ^ 0.7 fm a RHIC, et Oq ~ 0.3 fm a LHC. 



Conclusion 



Dans la premiere partie de cette these, nous avons etudie la possibilite qu'un condensat 
chiral dcsorientc soit forme lors du passage rapide de la transition de phase chirale dans 
une collision d'ions lourds. Nous avons propose une methode originale d'echantillonnage 
des conditions initiales pour le champ de pion en equilibre thermique local dans une petite 
bulle formee lors de la collision. Dans les travaux antecedents, ces conditions initiales sont 
habituellement choisies arbitrairement. Nous avons ainsi pu calculer une toute premiere 
estimation de la probabilite de formation d'une configuration classique du champ. En fait 
nous en obtenons une limite superieure, qui s'avere relativement faible, typiquement de 
I'ordre de 10~^. Ce resultat a des consequences importantes aussi bien experiment ales 
que theoriques. En particulier une strategic de detection evenement-par-evenement est 
preferable. 

Nous avons ensuite etudie la structure d'isospin de la configuration du champ produite 
par trempage des fluctuations initiales, et avons montre que, contrairement a ce qui a ete cru 
jusqu'a maintenant, celle-ci n'est pas une configuration DCC : les directions d'oscillation 
des differents modes dans I'espace d'isospin sont statistiquement independantes les unes 
des autres. Le modele le plus simple utilise jusqu'a present permet d'expliquer la nature 
classique du DCC, pas sa polarisation collective hypothetique. Nous avons montre que la 
dynamique ne genere pas les correlations cherchees entre les modes, une description plus 
realiste de I'etat initial pour le champ de pion est necessaire pour trancher la question de 
la possibilite de formation d'un DCC dans le scenario du trempage. 

Dans la deuxieme partie de la these, nous avons etudie la question de la thermalisation 
des gluons produits dans les tout premiers instants de la collision. Nous avons considere 
les collisions clastiques de petite deviation et modelise I'equation de Boltzmann correspon- 
dante par une approximation de temps de relaxation ou le temps de relaxation est calcule 
de fagon auto-coherente. Nos resultats sont en accord semi-quantitatif avec la solution 
exacte recemment obtenue pour une condition initiale donnee. Nous avons compare dif- 
ferents scenarios proposes dans la litterature pour decrire I'etat initial. Nous arguons que 
les criteres utilises dans les travaux precedents pour caracteriser I'ecart a I'equilibre local 
n'est pas satisfaisant. En mesurant plutot le degre d'anisotropie de differentes observables, 
nous arrivons a des conclusions qui contredisent celles precedemment obtenues. En partic- 
ulier, nous montrons que les collisions elastiques ne sont pas suffisantes pour thermaliser le 
systeme aux energies de RHIC. A LHC, I'equilibre est atteint dans le scenario de saturation, 
les incertitudes de I'approche, en particulier celles liees a la fragilite du calcul perturbatif, 
ne permettent pas de conclure dans le scenario des minijets. 

II est necessaire d'aller phis loin et d'inclure les contributions des processus inelastiques 
de branchement. Cependant, notre etude indique qu'il faut d'ores et deja s'attendre a ce 
que le regime transitoire ne soit pas completement negligeable a RHIC. Nous prevoyons 



d'inclure les quarks dans la description afin d'estimer I'influence de ces efFets hors d'cquilibre 
sur les spectres en energie transverse des dileptons produits, qui seront etudies par les 
collaborations PHENIX a RHIC et ALICE a LHC. 
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